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SYNOPSIS 

Pure nuclear quadrupole resonance (NQR) spectroscopy of 
solids is generally characterised by broad lines with widths of 
the order of kilohertz, in contrast to high resolution nuclear 
magnetic resonance (NMR) spectroscopy where the line widths are of 
the order of a fraction of a hertz. This is mainly because the 
electric field gradient (efg) inhoraogeneity and dipole-dipole 
interactions between the quadrupolar nucleus and other nuclear 
dipoles are not averaged out as in fluids. Dislocation, strains, 
presence of impurities, torsional motions of units in the crystal 
and fluctuations in the temperature cause a random distribution of 
the efg at the quadrupole nuclear site. Hence, there will be a 
distribution of quadrupole frequencies (aMq) aroundthe central 
frequency (o)q) in NQR spectra. In some cases, the presence of 
chemically inequivalent sites leads to multiple lines in the NQR 
spectrum. 

Uniform excitation of such broadband NQR spectra is not 
possible by means of a single radio frequency (rf) pulse of finite 
width and usually limited power. In addition to this, the presence 
of rf field inhomogeneity (Ao) ^ ) results in poor NQR signal inten- 
sity. Similar problems have been overcome to a great extent in 
NMR spectroscopy using composite pulses . These are a cluster of 
rectangular rf pulses whose amplitude and phase are piecewise 
constant and are applied without any time delay between them. 
Composite pulses consisting of n rf pulses can be completely 
defined by 2n variables, namely, n flip angles and n rf phases. 
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signal intensity from powder specimens. 

The thesis consists of seven chapters and five appendices. 

The first chapter is a brief introduction to NQR spectroscopy, 
where the possible role of composite pulses in NQR has also been 
outlined. The necessary quantum mechanical background for study- 
ing the spin dynamics in NQR spectroscopy is also presented in this 
chapter . 

Chapter II considers various techniques that have been 
successful in deriving NMR composite pulse sequences with a view 
to assessing tneir suitability for NQR spectroscopy. One cannot 
use the geometrical approach in NQR since, unlike in NMR of I = l/2 
case, one cannot easily visualize the spin trajectories for spin 
I 1 cases. A multidimensional search of a n flip angles and n 
rf phases defining an n-pulse sequence using a purely numerical, 
analysis approach involves an enormous amount of computer time and 
hence is not recommended for NQR spectroscopy. Iterative schemes 
(recursive expansion, reversed nutation pulse and fixed-point 
methods) yield, in general, sequences with a large number of rf 
pulses. As pointed out earlier, T 0 values in NQR being shorter 
than in NMR, it is not desirable to construct larger duration 
sequences for NQR. Moreover, the mechanism by which composite 
pulses overcome the limitations of a single rf pulse is not clear 
with these procedures. The Magnus expansion approach for NMR 
composite pulses, proposed by Tycko (Phys. Rev. Lett, 775 

(l983))> has been adapted here for NQR because it has the follow- 
ing advantages: (t) it reduces the problem of finding composite 
pulses to that of solving a set of equation, ( 2 ) it treats the 
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composite pulse sequence entirely rather than a pulse-by-pulse 
analysis, and (3) it can be used for nuclei with any spin I value. 
A general theoretical procedure for constructing composite pulses 
using the Magnus expansion technique in the case of NQR spectro- 
scopy is given in Chapter II. 

Chapter III presents the design of NQR composite pulses for 
physically equivalent spin 1=1 nuclei in single crystals using 
the Magnus expansion approach to compensate for efg and rf field 
inhomogeneities . The fictitious spin- 1 / 2 operator formalism of 
Vega and Pines (J. Chem. Phys. j6 6, 5624 (l977)) has been employed 
to carry out the spin dynamics calculations. A set of zero th- 
order composite pulse sequences consisting of two and three rf 
pulses has been obtained to overcome efg and rf field inhomoge- 
neities, and the performance of these composite pulses is demons- 
trated by computer simulations. If one limits oneself to phase- 
alternating rf pulses in the sequence then the procedure of deriv- 
ing composite pulses using the Magnus expansion approach is simpli 
fied and such sequences are also easy to implement experimentally. 
A group of phase-al ternating composite inversion 71 and 71/2 pulse 
sequences, including the "sj'mrnetric" ones i.e., symmetric with 
respect to the centre rf pulse in the sequence, has been derived. 
Computer simulations show that all the sequences reported here 
perform much better than a single rf pulse to compensate for A£o n 
and/or A(0-j* It has been demonstrated through the derivation of 
symmetric phase-al ternating composite pulses that there exists a 
one-to-one correspondence between the spin I = l/2 NMR case and 
pure NQR of spin 1=1 case. Results presented here are shown 
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to be independent of the asymmetry in efg. 

Chapter IV is concerned with the generalization of the appli- 
cability of composite pulse sequences to an arbitrary spin I 
nucleus in NQR spectroscopy. It has been proved that the compo- 
site / 2 and tt pulses derived using the Magnus expansion technique 
to compensate f or A(Oq and/ or Aco^ for physically equivalent spin I = 

3/2 nuclei in single crystals are identical with those of the 
sequences constructed for spin 1=1 nuclei. Here we have applied 
the irreducible spherical tensor operator formalism proposed by 
Bowden et al . (j. Magn. Reson, 67 . 4 1 5 ( 1986) ) for studying the 
spin dynamics . 

Chapter V develops and exemplifies the procedure for calcu- 
lating rf pulse responses in powders. A numerical approach has 
been employed to construct composite 71/2 pulses for powders. Compo- 
site 71 / 2 pulses thus derived have a great potential in NQR spectro- 
scopy since a large number of samples investigated in practical NQR 
spectroscopy are in the powder form and the use of these sequences 
should result in better NQR signal intensity in comparison to 
intensities obtained using single rf pulses. 

In Chapter - VI, experimental verification of the performance 

of composite pulses derived using the Magnus expansion technique 

is presented. The performance of both composite 71 and 7 t /2 pulses 

35 

has been evaluated by measuring the Cl NQR signal intensity as 
a function of resonance offset (Ao>q) and flip angle errors at room 
temperature in samples of powdered HgCl^ , single crystalline and 
powder specimen of NaClO^. Experimental results confirm that almost 
all the composite 7 T pulses function much better than a single tt 
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pulse in all the systemswe have tested. But, somewhat surpri- 
singly, the performance of composite 71 /z pulses is markedly 
excellent in HgCl,, powder, while not being effective in NaCIO 
Representing each pulse of slip angle 9 and phase 0 as (e^/the 
quence (A5) 1 g 0 ~ ( 135) 0 is recommended as a composite n/Z pulse 
while sequences (*»)„- ( 270 ) 9Q - ( 90) 0 and ( 3 1 5 ) Q - ( ) , SQ - ( 90 ) Q are 

recommended as composite * pulses for routine experimental sort in 

Nqft. 

Ihe last chapter of the thesis summarizes the present work 
end concludes with a brief outline of the scope for future work 
in the area of composite pulse NQR spectroscopy. 

Appendix A presents the definition of fictitious spin-l/2 
operators and their commutation relations, while Appendix B deals 
with irreducible spherical tensors for spin 1 = 3 /2 case and their 
commutation relationships. The importances of interaction frames 
in magnetic resonance have been explained in Appendix C and 
Appendix D discusses the general rf pulse response evaluations in 
NQR spectroscopy. Appendix E contains details of the derivation 
of composite n/Z pulses for quadrupole echo without phase distor- 
tions in the case of MMR „f spin I = 1 nuclei using the Magnus 
expansion method. 
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CHAPTER I 


INTRODUCTION 


I. A Nuclear Quadrupole Resonance (NQR) Spectroscopy 
I.A.i Phenomenon and Detection 

The interaction between the electric quadrupole moment, eQ, 
of a nucleus with spin 1^1 and eq, the electric field gradient 

(efg) surrounding it causes the splitting of energy levels even 

1 * 

in the absence of magnetic field. This interaction, defined by 
the quadrupole Hamiltonian, is given by^ ’ ^ 

A = 4x(^i-i 7 E 3I z 2 -i 2 .d (I x 2 -I y 2 )] ..(1.1) 

in terms of the angular momentum operator I. ^ is the asymmetry 
parameter, defined as 

V -V 

n _ _SE XX 

‘ ~ V 

zz 


.. ( 1 . 2 ) 


are the components of efg tensor defined such that 


zz 


> V S 


XX 


Since the electric field surrounding the 


nucleus is produced only by charges external to the nucleus the 


♦References appear at the end of each chapter 
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Laplace equation, given below, holds. 


V + V 
xx yy 


+ V 


zz 


0 


• • (-^ » 3 ) 


Thus n ranges from 0 to 1 . If T) = 0, then the efg is axially 
symmetric. The term eqQ in Eqn. (i.l), where eq = V , is 

ZZ 

referred to as the nuclear quadrupole coupling constant. 

In a representation where Hq is diagonal and = 0, the 
quadrupolar energy levels are expressed as^ 

E m = 4l(2I-1 ) [ 3 m — !(l + 1 )J » ..(1.4) 

where m is the z-component of the nuclear spin I. The energy 

levels are doubly degenerate in m except for the m = 0 energy 

level in the case of integral spins. The presence of T) mixes 

states corresponding to An* = ±2 and the m's are no longer 'good' 

quantum numbers. In addition, 77 lifts the degeneracy of _+m energy 

levels in the case of integral spins. For various values of the 

nuclear spin I, the secular equations and corresponding eigen- 

13 4 

values are available in the literature . * Transition between 

the energy levels established by ^ can be produced by the appli- 
cation of an oscillating radio frequency (rf) field which inter- 
acts with the magnetic dipole moment of the nucleus. The NQR 
phenomenon has been observed by various techniques including conti- 
nuous wave (CW) and pulsed Fourier Transform (FT) methods,'**^ 
as in Nuclear Magnetic Resonance (NMR) spectroscopy. In this 
thesis, we shall be concerned with only pulsed NQR . 
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If one applies an rf pulse of a frequency go n , to the 
system under study for a short time, t^, (known as rf pulse width), 
with an rf field strength w ^ (rad/sec) (where f = gyromag- 

netic ratio; = rf field amplitude in Gauss) and detects the 
following NQR signal in time domain, the observed response of the 
system decays with a time constant (T^) known as spin-spin relaxa- 
tion time. This response is called Free Induction Decay (FID). 

By Fourier -transforming this FID, one gets the NQR signal in the 
frequency domain.** Generally broad (line width of the order of 
kilo Hertz) and multiple lines are observed in the NQR spectrum 
because of various internal interactions present in the system. ^ 
Hence, using a single rf pulse with finite power, it is not 
possible to excite these broad spectra uniformly. Exactly with 
this problem in view, we have undertaken the task of constructing 
broadband excitation rf pulse sequences for NQR spectroscopy in 
zero dc magnetic field, which forms the bulk of the subject 
matter contained in this Thesis. 

For a better understanding of the problem we present, in 
Section I. A. 2, various line broadening mechanisms which charac- 
terize NQR spectroscopy. The need for broadband excitation in 
NQR will also be illustrated in Section I.B. A powerful theore- 
tical tool for studying the pulsed response of quadrupolar 
systems is the density matrix approach, and we present the back- 
ground for such spin dynamics calculations in Section I.C. 

Finally, this chapter ends with an outline of the present thesis 


in Section I.D 
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I. A. 2 Line Broadening Effects in NQR 

In pure NQR spectroscopy (i.e., zero applied magnetic field), 
line broadening originates from two main sources; these are the 
(i) electric and (ii) magnetic interactions present in the solid 
material. In this section, a brief discussion of these inter- 
actions is given. 

I.A.2(i) Electrical Sources of Line Broadening 

NQR frequency corresponding to the coherence between any 

two energy levels in a spin I quadrupole nucleus is given by 
2 1 

(Oq a © qQ» where q is the efg component. Effects which cause 
a change in efg-s of quadrupole nuclei lead to a distribution of 
NQR frequencies (go q _+ Agoq)* Eor example, strains and disloca- 
tions in the crystal or in powder grains create a random distort- 
ion of intermolecular interactions and a corresponding random 

distribution of efg is evidenced at quadrupole nuclear sites. 

7-9 

Strains which result from the crystal growth process' and the 
presence of impurities in the crystal lattice^ * ^ have been 
found to broaden NQR lines. In NQR literature, the effects of 

1 

molecular motion upon quadrupole frequencies have been studied. 

There are two kinds of molecular motions, namely, torsional and 

the hindered rotational motions of the molecule about some axes 

1 11 

in the crystal. Baeyer has shown that torsional motions of 
units constituting the crystal lattice lead to (a) a shift in the 
resonance frequency, due to fluctuations in efg and (b) a change 
in the relaxation process affecting quadrupole resonance. Both 
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of these are, of course, temperature dependent. As such, NQR 
frequencies are highly sensitive to changes in temperature which 


affects the random motions of units in the crystal lattice. 


12 


Generally, the temperature coefficient of the NQR frequency. 


da Q 

dT 


i is negative. This effect has been explored in the cons- 


truction of NQR thermometers which can measure temperatures 
12 13 

accurately. ' Overall, these abo vementioned effects are 
responsible for experimentally observed line widths (often they 
are of the order of a few kiloHertz to a few megaHertz) . This 
means that the spin-spin relaxation time T and spin-lattice 
relaxation time, T^ (which play very important role? in deter- 
mining NQR line widths), are much shorter J 


Apart from the inhomogeneous broadening^ of NQR lines, 
electrical interactions create multiple spectral lines in an NQR 

3 "| Kmm “| *7 

spectrum. * Often, the presence of chemically inequivalent 

sites^ in the crystal specimen leads to multiple lines in the 
pure quadrupole resonance spectrum itself. Here, two or more 
quadrupole nuclei with spin I>1 are said to be chemically inequi- 
valent when they have different chemical environments or different 
efg-s in the crystal lattice. Two chemically inequivalent sites 
cannot be related by any kind of symmetry operation. Chemical 
inequivalence may arise due to the occupation of atoms with 
quadrupole nuclei, at chemically different environments in the 
free molecule itself or it may also originate even in cases where 
the atoms are in equivalent positions in the free molecule but 
the interactions between neighbouring molecule in the solid destroy 
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their equivalence J The splittings produced by the presence of 
chemically inequivalent atoms in the free molecule are much larger 
than those produced in the latter case, because of the comparative 
weakness of intermole cular forces relative to that of the forces 
that exist inside the molecule. In NQR literature, multiple lines 
for a given quadrupolar nuclear species have been observed in many 

<j ^ -a ry 

systems and has been attributed to chemical inequivalence. * ' 

3 5 

For example, at 77 sixteen lines have been detected in the Cl 

8 

NQR spectrum of solid CCl^ with a spectral width of a> 350 kHz. 

I.A.2(ii) Magnetic Interaction s 

There are two kinds of magnetic interactions between a resonant 

nucleus and other non-resonant nuclei of the same or different 

1 7 

species. They are, namely, (a) direct magnetic dipole-dipole 
interactions, and (b) indirect spin-spin interaction, also called 
J-coupling The spin-spin coupling is usually much smaller 

than the dipole-dipole interaction and is not of interest to us 
in the present work. Hence we consider only the dipole-dipole 
interaction here. 

(a) Dipole-dipole interactions : Dipole-dipole interaction is the 
interaction between the magnetic moment of a resonant nucleus with 
the dipolar field due to the magnetic moments of all other spins 
in a rigid lattice. This depends on the magnitude, orientations 
of the magnetic moments and also the lengths and orientations of 
the vectors describing their relative positions. When the quadru- 
pole nucleus has one very close neighbour, then the dipolar inter- 
action with it predominates over that of distant neighbours and 
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fine splitting may be observed in NQR spectrum, similar to the 

20 

' Pake doublet* in NMR. Slow beats in the spin echo envelope 

22 

have been observed in solid N^, because of the intramolecular 
dipolar interactions. In the chlorine NQR spectrum of solid 
chlorine ( T] £ o),^ due to "^C1-"^C1 and ^^Cl-^^Cl homonuclear 
dipolar interaction a doublet has been detected. He teronuclear 

dipole-dipole interaction has been observed by Maclcowiak and 

22 35 

Brookeman for HC1 as fine splittings in the Cl quadrupole 

line. Other typical examples may also be found in the NQR 

literature 

I.B Necessity for Broadband Excitation in NQR 

As pointed out in the previous section, NQR spectral lines 
are inherently broad and the presence of multiple lines leads to 
a broadband NQR frequency spectrum. Usually, in FT NQR, a single 
rf pulse with a constant frequency and power is used to excite 
the whole wide spread spin populations so as to obtain the entire 
NQR frequency spectrum.^ However, it is often not possible to 
excite the broadband NQR spectrum uniformly by a single rf pulse. 
Spins which are off-resonant with respect to the rf pulse frequency 
are poorly excited and hence their resultant signal intensity will 
be poor. This problem can be partly overcome by applying high rf 
power. But, because of difficulties in the experimental technique, 
the production of high power rf pulses is limited. 

In addition to this, the response of a single rf pulse is 
highly sensitive to rf field inhomogeneity ( Aoij ) across the NQR 
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sample placed in an rf coil. Because of this 0)^ hh distribut- 

ion in the system, spins in different parts of the sample feel 
different flip angles ( 9± A B) which are directly proportional 
to the strength and time duration of an applied rf pulse. Thus, 
the observed signal intensity is poor and, in powder specimens 

particularly, it is difficult to find out the exact tt/ 2 and n 

rf pulse time durations in an experiment. So it is essential to 
find a suitable replacement for single rf pulse in FT NQR experi- 
ments. Somewhat similar problems, such as resonance offset and 

rf field inhomogeneity, have been overcome by the use of ' composite 

2k 

pulses 1 instead of single rf pulses in NMR spectroscopy. A 
’composite pulse' is a sequence of phase shifted, constant ampli- 
tude rectangular rf pulses, without any time delay between them, 
to excite the spin populations over a larger range of some experi- 
mental parameters compared to a single rf pulse. To the best of 

our knowledge, no attempt has so far been made in the realm of NQR 

spectroscopy to design 'composite pulses'. Therefore, a programme 
of designing composite pulses that compensate for efg inhomoge- 
neity or resonance offset ( Ac<^) and rf field inhomogeneity ( Aw^ ) 

in NQR spectroscopy seems worthwhile. To construct such rf pulse 

2k 

sequences different routes have been analysed from the NQR point 
of view. These will be presented in Chapter II of this thesis. 

As a convenient background for understanding the theoretical 
procedures involved in this work, we present the spin dynamics 
studies using the density matrix approach in the next section. 
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I . C Spin Dynamics Studies 

W _ 

I.C.1 A Brief Survey of* Theoretical Procedures ,-Bhpploved in NQR 
Spectroscopy for Spin Dynamics Studies 

In this subsection, we briefly survey the different proce- 
dures used in NQR literature to follow the time evolution of the 

25 

spin system. In 1 955 > Bloom et al . used the time dependent 
perturbation theory to evaluate the state of the spin system in 
the presence and absence of rf pulses. Later, Das et al. 
developed the density matrix mathod for calculating the FID and 
spin echo signal for a general spin I in the case of single 
crystals. Pratt et al . have investigated the response of a 
spin 3 1 2 system in NQR to (i) a long off-resonant rf pulse ' and 

2g 

(ii) two-pulse sequence. In these cases, the response of the 
quadrupolar nucleus has been proved to be similar to that of 
spin 1/2 nuclei in NMR experiments. In 1977» Pratt has intro- 
duced an interaction representation approach to analyse the 
Slichter-Holton and Lee-Goldburg pulse sequences in spin 3/2 NQR 
spectroscopy. A common difficulty with these methods is in obtain- 
ing complete solutions for differential equations at different 
times. Using these methods one cannot explain the experimental 
results of many newer pulse techniques in NQR. Furthermore, it 
is not necessary to consider the entire form of the density 
matrix. This problem can be simplified by defining the density 

matrix of a system containing spin I nuclei in terms of properly 

2 

selected (21+1 ) -1 traceless independent Hermitian operators. 

Such a set of operators, which give rise to the very useful 
'fictitious spin-1 /2* operator formalism , J 


has been proposed by 
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Vega and Pines^ for the use of multiple quantum studies in spin 

1=1, NMR spectroscopy. This has been effectively used later 

in spin 1 NQR spectroscopy to study the Overhauser effect by 
32 

Vega and also for analysing the effect of various pulse 

33 34 

sequences. Provo torov et al. have used this formalism 

along with the procedure of canonical transformation to analyse 

the line narrowing effect of multiple pulse sequences in NQR 

spectroscopy. Rue to the simple commutation relationship of 

these operators, we have adapted this formalism in our procedure 

3' 

for designing composite pulses for application to spin 1 cases. 

A detailed description of the fictitious spin-1 /2 formalism is 

presented in Appendix A of this thesis. This formalism has been 

generalised in recent years for higher spin cases by Vega J and 

37 

tfokaun and Ernst independently. A generalised version of this 


formalism has been applied for studying the behaviour of multiple 


pulse sequences in spin 3 /2 case, 


Projection operators have 


been employed in a theoretical study of multiple pulse sequences 


in NQR spectroscopy. 


Recently, Singh and Armstrong have 


employed SU ( 3 ) group generators for investigation of the zero- 

40 4 1 

time resolution technique and Jeener-Broekart pulse sequence 

in spin 3/2 NQR spectroscopy. Irreducible spherical tensor 
42 

operators have been used by Sanctuary et al. to solve the 

4 

Liouville equation directly for the cases of 1=1, 3/2 and 5/2 

4 ' 

and to study the selective excitation in pulsed NQR spectroscopy. 
Somewhat similar approach has been used in a theoretical argument 


questioning the existence of the phenomenon 1 quadrupole precession' 

/ 44 

in spin 1=1 and 3/2 cases in the absence of Zeeman field. 



A modified version of spherical tensors, proposed by Bowden and 

45 

Hutchison, has been used for multiple quantum NQR studies in 
/ 46 

1=5/2 case. Since spherical tensors form a complete basis 
set and their commutation relations are well established, we 
have employed them in the procedure for constructing composite 
pulses for spin 3/2 case in this thesis. In Appendix B, we have 
presented an introduction to irreducible spherical tensors, their 
commutation relationships and their matrix representations. Since 
the present work employs the density matrix method, it is felt 
that a brief discussion about this method is essential, and 
this has been done in the following section. 


X.C.2 Density Matrix Method 
I.C.2(i) Density Matrix 

We consider a multilevel spin system with a complete set 
of eigen states The state of the system can be written 


as 


j'l|/^>= E Cjj. * A mixture of independently prepared 


states l^n> (n = 1,2,...) with statistical weights can be 
completely described by an operator, jj , as 


T = l v n 


.. ( 1 . 5 ) 


The operator a is generally called the 'density matrix' or 

47 

'statistical operator'. <3 may also be viewed as a vector in 
Liouville space . ^ a is Hermitian, that is, <^^i|d|0j^ = 

. The diagonal element (o r mm ) of the density matrix 
corresponds to the probability of finding the system in the state 
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io- Trace of 0 , Tr (a) , is equal to 1, shoving that it is 
independent of the choice of the representation. • The main 
advantage of using the density operator is that all information 
on the behaviour of a given system can be expressed in terms of 
expectation values of suitably chosen operators. For example , 
the expectation value of any operator Q is given by 


= Tr(aQ ) 


.. ( 1 . 6 ) 


49 

This approach has been introduced by Fano in 1957» * Since then, 

it has been effectively used for the spin dynamics studies in 

, . 50-52 

magnetic resonance. 

Thermal equilibrium density matrix, a (o) , can be given 
according to the Boltzmann distribution of populations among 
the energy levels. In our case, an ensemble of quadrupole nuclei 
with spin I>1 in the absence of a magnetic field can be defined 
as 

}" 

50 

Using the high temperature approximation, we write 

-1 

CT (0)s[Tr( 1 )] (l -^q/kT) .. (1.8) 


(o) = exp 


r 

( ) 4 Tr[ exp(-# Q /kT)] 


where is the quadrupole Hamiltonian given by Eqn. (i.l), k is 
the Boltzmann constant and T is the temperature. Off-diagonal 
elements of (j(o) can be taken to be zero according to the random 
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CO 

phase approximation . ~ >J The first term in Eqn. (E.8) is unaffected 
by an evolution of the system, so we need to follow only the evolu- 
tion of the second term. For simplicity, we write 

a (0) » [Tr (£ )] 1 [i. -£^1] .. ( 1 . 9 ) 

where P(o) is the reduced thermal equilibrium density matrix 
defined as 


P(0) = % 


Q 


.. (i.io) 


I.C.2(ii) Time Evolution of the Spin System 

The time development of the spin system under various inter- 
actions defined by the Hamiltonian, % , can be studied using the 
equation of motion of the density matrix^ written as 


Jj S il ■ [#. M 


.. (i.n) 


This can be derived from the time dependent Schrtfdinger equation 

47 

in the Hilbert space. This differential equation is often 

called the von Neumann equation or Liouville equation (because it 

assumes the same form as the equation of motion for the phase 

\ 48 

space probability distribution in classical mechanics). The 
solution of Liouville -von Neumann equation is given as 


P(t) = U(t) p(o) U(t) 

t„ 


-1 


where U( t) = T exp [ -i J %{t') tit'"] 


.. (1.12) 

.. (1.13) 
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47 

is called the Liouville operator or Dyson expression or the time 
evolution operator corresponding to the interactions given by . 

T is the Dyson time-ordering operator, which orders operators of 
longer time arguments in the expanded exponential to the left. For 
example, 


T {ft(t* ) M{t ") j = 


•jUt* ) ft (t") if t ' > t" 


# ( t ” ) !&(t f ) if t M > t ' 


.* (I. 14) 


During such an evolution of the density matrix: 

(a) the eigen values of the density matrix are time independent; 

(b) the diagonal elements of the density matrix will not change 
with time and (c) its off-diagonal elements, which correspond to 
the coherence between different energy levels, oscillate with 
time * 

In magnetic resonance, the time evolution of spin population 
can be thought of as rotations of spin vectors in a (21 4* 1 ) 
dimensional, space* Commonly, exponential operators are used to 

54 

explain such time evolution and their properties are as follows: 

(a) Evolution operator U(t) unitary, i*e., TJU*= J. ** (I. 15 ) 

2 3 

(b) e A = 1 + A + jj- + jj- + . . . .. (1.16) 

(c) e A = e A+ ^ = e A if and only if [A, If] =0 .. (I. 17 ) 

A B 

where A and B must be anti-.Hermi tion so that e and e will 


be unitary 



15 


(d) 


, . 55 

Baker-Campbell-Hausdorff ^BCH) formula to second order is 
given for [ A,b] 0 case as 


e A e® = exp-£ A + B + l/2 ( A+B) + l/l2 [A, ( A+B) ] + 

1/12 [[A,B], B] + ... J .. (1.18) 


Using the BCH formula, one can write. 


iAd „ -iAa 
e Be 


B, if [ A,B] = 0 


.. ( 1 . 19 ) 


This explains the constant of motion of B under the effect of A 
and 

e iAa B e '*‘ Aa = B cos (a) -C sin(a) , if [ A,B] = iC; cyclic. .. (l.20) 


One may also independently derive the above result using the 

56 

following nested commutation relationship 


s iAa B e" iAa = B + id [ A,B] + [A, [A, B]] + 


i±al 3 


31 


i £ A,£ A, £ A , Bj + ... ., (l.21 ) 


In some cases, evaluating such a time evolution becomes a diffi- 
cult task. For this purpose, computer programs have been develop- 
ed in NMR literature for the automatic generation of the commuta- 
57 

tor algebra. Howevex, in solving time dependent equation of 
motion, one has to make approximations to get time independent 
solutions. In the next subsection, we review the different 
strategies that have been used in magnetic resonance for studying 
time dependent perturbations. 
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I.C.3 Methods for Solving Time depe ndent von Neumann Equation 

Evolution of the spin systems under the effect of a time 
dependent Hamiltonian can be brought out by 

.. ( 1 . 22 ) 

In magnetic resonance, the applied rf radiation is time 
dependent and in some cases we often make the line broadening 
interactions time dependent by transf ering into a suitable inter- 
action frame. For example, in average Hamiltonian theory 
(AHT) ’ which is used to explain the multiple pulse effects, 
one uses the rf interaction frame to induce a time dependence on 
the internal interactions called ^ t ) and it is time— averaged 
out to achieve line narrowed spectrum. Similar is the case with 

C O 

the construction of composite pulses# In general, however, an 
exact solution of Eqn* ( I * 22) is not possible* Often, the inter- 
action *( t) is a small perturbation and Eqn. (i. 22) can be solved 
by approximation. Some of the methods applied in the literature 
are presented below. 

X.C.3(i) Time dependent Perturbation Theory 

47 

In this theory, one gets an iterated solution given by 

u(t) = i- 4 o / dt i ^ (t i> + ( f )2 7 dt i t y dt 2 **(* 2 ) + 

.. (1.23) 

The disadvantage with this solution is that it is not unitary. 


k JL £r 1 =[«(*). p(t>] 
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But in magnetic resonance, unitary operators in the exponential 
form are needed because they do not change tr(P) and they corres- 
pond to rotations in a (21+1 ) dimensional space for spin I case. 
Moreover, this solution is expected to be accurate only for small 

K(t). 


I.C.3(li) 


In 1954, Magnus'* solved the time dependent differential. 


equation with a procedure which is essentially a rearrangement 
of the perturbation series. In this procedure, the major assump- 
tion is that 'H (t) is constant in smaller intervals of time t^ , 

^ 2 » *•* with t = tj + t^+ ... + t n . t) is expanded to an 

infinite series as 


#(t) = 1 ^ + ... + ... .. (1.24) 


where ' is the nth-order term in the series ♦ The zero th- 

47 

and first-order terms are expressed as 


« (o M /*<*■> 

o 


.. (1.25) 


and J dt 2 J 2 dt,[ ^(t 2 ), %(t t )] ..( 1 . 26 ) 


Similarly, an nth-order term is given by (n+l) integrals with n 
commutators of *(t). Here, each term is a Hermitian operator 
and the infinite series can be limited to any order so as to get 
a unitary solution. In this series, is the main term and 

other higher order terms are quantum corrections to the classical 
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solution , Hamiltonians corresponding to different time 

intervals do not commute with each other. As Magnus has pointed 

59 5 *s 

out, the Magnus formula is the extention of BCH formula, J which 

is for ^C(t) with only two different time intervals. The exponen- 
tial solutions for a time dependent differential equation has been 
also derived independently by Robinson in 1963^ and by Pechukas 
and Light in 1 965 • ^ ^ 


Magnus expansion has been effectively used to synthesize 

51 52 58 

multiple pulse^ * and composite pulse^ sequences. Warren and 

_ 62 

Zewail have extended these techniques into the optical region. 
Cady has used it in the theory of spectral line broadening and 

64 

Sheck et al . have applied Magnus expansion in the theory of 

65 

multiphoton absorption. Cross used it in his formulation of 

66 

scattering theory, while Saxon and Light have applied it in the 

treatment of various problems in molecular scattering and energy 

transfer. It has even been used in the formulation of the Born- 

6 7 

Oppenheimer approximation. It suffices to say here that the 

present literature is extensive on the applications of the Magnus 

expansion technique. However, several groups have reported that 

the application of the Magnus expansion in the Schrddinger 

representation to explain the observed experimental results was 

less than satisfactory. ” The reason for this could be due to 

the possible errors in the forms of the Magnus expansion terms 

used. It may be mentioned here that there exists a number of 

derivations of the Magnus expansion terms^’^ and there are 

variations in the resulting definitions of different order terms 

71 

between authors. Later, in 1985 , Salzman developed an 
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alternative to the Magnus expansion and lie notes that it is easier 
to use than the original Magnus expansion itself. His procedure 
leads to a straightforward derivation of unique higher-order 
Magnus expansion terms. Salzman concludes that the problems with 
Magnus expansion more likely arise from a lack of the convergence 
of the series. 

Magnus expansion approach has been successful in predict- 
ing the behaviour of the spin system at relatively short times; 

however, because of the convergence problem which is called as 

68 

•Magnus paradox', it fails to account properly for the long time 
70 72 

behaviour. * This is the case regarding its use in the analy- 

70 

sis of multiple pulse spin-locking experiments on dipole-coupled 
many-body spin systems. In the literature, the convergence of 
Magnus expansion for a two-level case has been studied in detail 

7 

and even the radius of convergence has been extended considerably. 
The criterion proposed as a lower limit to the convergence radius 
of the Magnus series is given by 



dt < 2n 


. . ( 1 . 27 ) 


where + w(t)/2 are the time-dependent eigen values of the 
Hamiltonian #(t). 

Furthermore, the equilibrium properties of the system can 
not be determined in the Magnus expansion procedure by statis- 
tical mechanical techniques analogous to those for time independ- 

51 52 

ent systems. In AHT ’ based on the use of Magnus expansion. 
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the observation of the system are possible only at the end of a 
cycle (U ^.(t c ) = ±X )» where the rf interaction frame and the 
observation frame coincide. These problems have been to some 
extent overcome by the application of Krylov, Bogoli jubov and 
Mitropolski j approach (KBM)^ and Floquet theory*^ which are 
suitable alternatives to the Magnus expansion. KBM averaging 
method has been used to observe theoretically the slow motion 

evolution of the spin density matrix during the applied multiple 

75 76 77 

pulse sequences. Using the Floquet theory, Maricq * 1 analysed 

pulsed spin-locking phenomenon resulting with the application of 
multiple pulse sequences in magnetic resonance. Equilibrium pro- 
perties of a model time dependent Hamiltonian in an interaction 
reference frame have been extensively studied with the help" of 

rj /T 

Floquet theorem. For further details about spin dynamics in 
magnetic resonance, one is referred to the literature.^ *52,56,78 


It should be mentioned here that, as far as the construction 

of composite pulses is concerned, the use of Magnus expansion 

approximation is the suitable one and this has been completely 

24 58 

successful in NMR spectroscopy. * J In this thesis, we have 
applied the Magnus expansion technique , for designing composite 
pulses suitable to NQR spectroscopy, and this has been detailed 
in Chapter II. 

I .D Scope of the Present ¥ork 

Uniform excitation of broad quadrupole frequency spectrum 
in zero dc magnetic field is necessary (i) to improve the NQR 
signal intensity and (ii) while searching for chemically 
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inequivalent sites present in the crystal system by observing the 
multiple lines of NQR spectrum. This would be of greater impor- 
tance in investigating high temperature superconductivity materials 
whose NQR line widths are of the order of a few MHz.^ It will be 
also useful in improving the and measurements with better 
inversion and refocussing of spin populations over wider range* of 
in spin inversion recovery and spin echo methods, respectively. 
Similarly, excitation over a larger range of AoXj would compensate 
for flip angle errors of rf pulses employed in pulsed techniques 
of NQR. For this purpose, in the present work, we propose compo- 
site K and IT/ 2 pulses for the replacement of single rf pulses in 
compensating efg and rf field inhomogeneities in NQR spectroscopy. 

In Chapter II, various routes used in NMR spectroscopy for 
designing composite pulses have been analyzed from the NQR point 
of view and a general theoretical procedure based on the Magnus 
approximation is given for NQR of spin I nuclei. 

Chapter III deals with the design of composite pulses 
using the Magnus approach for spin 1=1 nuclei. Fictitious 
spin- l/2 operator formalism is employed for the spin dynamics 
calculations. In Chapter IV, this procedure is extended for a 
general spin I nucleus by applying spherical tensors. Composite 
7T / 2 pulses that can enhance the NQR signal intensity in poly- 
crystalline samples have been designed using a purely numerical 
approach in Chapter V. Experimental verification of composite 
pulses designed using the Magnus expansion method, on powder 
specimens of sodium chlorate (NaClO^) and mercuric chloride 
(HgCl^) have been presented in Chapter VI. 



In the last chapter of the thesis important conclusions 
and points for further work in this area of research have been 
given. Five Appendices have been included to clarify some aspects 
of the present work. 

Summary 

In this chapter, a brief introduction to the phenomenon 
of NQR and its detection as well as a discussion of various line 
broadening mechanisms in NQR have been given. However, no attempt 
has been made to cover the full literature of NQR spectroscopy. 

The need for designing composite pulses to overcome resonance 
offset and rf field inhomogeneity problems in NQR has been 
emphasized. This chapter also presents an introduction to the 
density matrix approach to the study of spin dynamics in NQR 
spectroscopy. The density matrix approach is of great value in 
understanding magnetic resonance phenomena in general and NQR 
spectroscopy, in particular. In this context, the Magnus expan- 
sion approach to the solution of time dependent problems has been 


critically examined 
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CHAPTER II 


TECHNIQUES • FOR THE DESIGN OF 
COMPOSITE PULSES FOR BROADBAND 
EXCITATION IN NQR SPECTROSCOPY 

There are various ways of designing broadband excitation 
composite pulses and they have been used effectively in NMR spec- 
troscopy. The purpose of this chapter is to discuss these methods 
briefly, especially in relation to NQR spectroscopy and to outline 
a thoretical procedure for constructing composite pulses for NQR 
spectroscopy using the Magnus expansion approach. This chapter 
is organised as follows. In Section II. A, the geometrical approach 
of deriving composite pulses has been outlined. Different iter- 
ative schemes have been considered in Section II. B, while in 
Section II. C, we deal with the 'coherent averaging* approach. 
Section II. D presents the details of the use of Magnus expansion 
in coherent averaging approach to construct composite pulses for 
NQR spectroscopy. 

II. A Geometrical Approach 

It is well known that the state of an isolated spin l/2 
nucleus can be represented by a vector in three dimensional space 
and the dynamics of the spin can be understood in terms of the 
rotation of the vector in an appropriate coordinate space which 
is the Bloch-vector model. * In this model, the thermal 
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equilibrium spin system in NMR can be defined by a vector along 
the magnetic field direction ( z-axis ) . Nov in the rotating co- 
ordinate a 90 ° rf pulse applied along the x-axis ( 90 °) rotates the 
magnetization vector from z- to y-axis. Similarly a 180° pulse 
will leave the vector along -z-axis. The Bloch- vector model has 
been very useful in NMR to visualize the spin dynamics in many 

experiments and has in fact been employed to follow the spin 

2 3 

trajectories under multiple pulse sequences. 5 J 

4 

In 1979, Levitt and Freeman used this Bloch— vector model 
to design composite pulses that compensate for resonance offset 
(A« o ) and rf field inhomogeneity (A 03 ^ ) . They have examined 
Bloch-vector trajectories as a function of A0) o and A03^ • From the 
geometrical pictures, the conclusion emerges that combinations of 
various flip angles and rf phases constituting the composite pulse 
have to be carefully chosen such that it cancels out the defects 
of a single rf pulse. Relaxation of spins is neglected while 
considering the motion of the nuclear magnetization subjected to 
a composite pulse in the Bloch-vector model. In this method the 
overall effect of a composite pulse is obtained by the solution 
of a simple problem in spherical geometry. Instead of seeking an 
analytical solution for this problem which is geometrically compli 
cated, computer simulated trajectories have been used. The geo- 
metrical method has been successful in developing composite pulse 
5-7 

sequences which are widely used in NMR today. 

The disadvantages of the geometrical approach are given 
below. The variety of composite pulses » and their degree of 
compensation achieved, have been restricted by the limits of 



human ability in visualizing a set of three dimensional rotations, 
each subject to errors . Though the extent of compensation can be 
increased by increasing the number of rf pulses in a composite 
pulse, no easy method has been proposed to construct such sequen- 
ces. Increasing the number of pulses in a composite pulse compli- 
cates the geometrical picture also. At a time, one can construct 
a composite pulse to compensate only one parameter and the assump- 
tion that the initial condition is always defined by a vector 
along the z-axis limits its applicability. Moreover, it is not 

J 

valid for coupled spin systems. 

From the NQR point of view, the geometrical method is not 
as attractive as in J NMR for constructing composite pulses. The 
simple vector model is not possible for spin I>l/2 cases, unlike 
in spin l/2 case. Imagination of a vector representing a spin I 
nucleus in a ( 21+1 ) ‘- 1 dimensional space is necessary, if one 
wishes to use the geometrical method for a spin I > 1 /2 quadrupolar 
nucleus. s 

II. B Iterative Schemes 

In this approach, rather than attempting to construct 

a sequence with the desired propagator directly, a set of opera- 
tions is defined th&t will be applied repetitively to an initial 
pulse sequence, thuS- generating a series of iterative sequences. 
The procedure of using iterative schemes is depicted in Fig. II. 1. 

The aim of iterative procedure is to construct the opera- 
tions in such a ways that the propagators for iterative sequences 
converge to a desired form. Constructing higher iterative 
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Fig. II. 1 Schematic illustration of* an iterative scheme . ^ 

Operations Eire applied repetitively to an initial 
pulse sequence S q , and generating a series of* 
iterative schemes , S^ etc. 

sequences becomes increasingly complex and it will be composed 
of increasing number of individual pulses. The higher the iterate 
sequence the better the performance will be. Iterative sequences 
perform independent of the choice of the initial sequence. For 
example, regardless of S q , may be a broadband excitation 
sequence. Therefore, the development and evaluation of an itera- 
tive scheme must focus on the operations that comprise the sequence 
rather than on the specific sequences that the scheme generates. 
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In NMR literature, operations treated include phase shifting of 
rf pulses,** 11 concatenation of sequences** - 1 1 permutation of rf 
pulses, 6 fixed-point theory, 11 ’ 12 formation of inversion 

8 ili it; 

sequences and pure numerical analysis* f J It is beyond the 
scope of this dissertation to explain all these procedures in 
detail but for the sake of completeness, a brief presentation is 
given here . 


II.B.1 Recursive Expansion Procedure 

g -j o 

This is an iterative scheme ’ J for constructing composite 
pulses of arbitrary flip angles, not necessarily rr/2 or tt radians, 
to achieve higher degree of compensation for resonance offset and 
rf field inhomogeneity in NMR. The procedure for generating a 
better composite n/2 pulse is as follows: 

(i) Choose a composite n/2 pulse, P Q ra ; m is the order of 
iteration. 

(ii) Find out the inverse sequence (P^ ) 1 by phase inversion 
and time reversal. 


(iii) Phase shift the inverse sequence in phase by ±n/ 2 to get 

(iv) Concatenate with the original sequence P Q m , in any order, 


to get an improved sequence like p m+1 = ( P m ' T ' m 


,-1 , 
±n/2 ' P o 


(v) Steps (i) to (iv) may be applied in turn to P™ +1 and the 
process repeated to get a desired sequence. 

(vi) At any stage, the composite 7t/2 pulse P^ may be converted 
iqto a composite pulse of flip angle j3 (arbitrary) by forming 
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( p rp-1 p m . for example , if p = 135°, then P™ +1 = 

O o 

is a composite 135 pulse. If P q = (rr/2) y , a single Tr/2 pulse 
applied along the y-axis of the rotating frame. 

p o - ( ,r /2) x ( ,t /2) y • 

This is a composite TT / 2 pulse which compensates Tor the rf 
field inhomogeneity effects* 


Starting from P q , we can construct a composite inversion 


pulse as follows : 


R 1 = (P 1 ) _1 P 1 

O V TT 7 o 


(7T /2) y (7t/2) x (7X/2) x (7T/2) y 


= (7t/2) y 0l) x (n/2) y . 

Sequences constructed by this procedure have been effect- 
ively used in high resolution NMR, NMR imaging and broadband 
decoupling experiments. * * 


II. B. 2 Reversed Nutation Sequences 

Following the recursive expansion procedure, Shaka and 

9 

Freeman constructed a set of composite invers ion TT pulses 
using a method which is similar to the recursive expansion 

9,13 

approach. In this method, called the reversed nutation method, 
starting from a single TT pulse, composite inversion pulses are 
obtained by the following procedure: 
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> ra #. P m (p™)" 1 P m = P m +i 


Let P = Tt 


x 


then P = 7T ^ Tty. TT x , is a composite inversion TT pulse. Similarly, 

starting from the sequence P with (P ) = tt tt TT , we get 

y —x y 

P 2 = P 1 (P 1 )” 1 P 1 as 

p2 “ "x \ "x "y *.x "y "x "y *x 

2 

P is also a composite tt pulse f whose performance is better 
than . 


HeB.3 Phase Shifting and Catenation Operation Scheme 

In this procedure, a class of iterative scheme ^ ° is defined 
by the following operations: 

(i) Consider an rf pulse sequence and form N phase shifted 

versions of S ± with phase shifts ^ , $ 2> . , . (where N 

is an odd integer)* 

(ii) A new sequence ^ is obtained by concatenat ing the N phase 

shifted versions in order from 1 to N as explained in an 
example below. Thus is N times longer than 

An example of constructing a composite tt pulse is given 
here. Ve assume that is acting on an isolated spin to produce 
a rotation R^ f given by 

= exp[ -i (l x cos$ o + I sing!S o )tt] exp(-i6l) .. (il.l) 

where <jj& Q is the rf phase and £ has a z-component of zero. is 
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in fact tiie propagator for S ± i.e., R ± = U ± (t). If the inversion 
is nearly complete then 5 is small. Applying the iterative 
scheme, a new sequence S i+1 is generated to first order in |6|. 
The propagator R i+ is given by 

R i+1 = ® x p[ -i(l x cos^ T + I sin^k] exp(-i6 T I ) .. (II.2) 


N 


where $ = + E (-1 ) n+1 $ 


n=1 


.. (II. 3) 


N 


and 


6 t = r {a x cos0 n -6 y sin^ n , (-1 ) n+1 6 y cosj 3 n +(-l) n+1 6 x sin0 n ,O 

.. (II. 4 ) 

/3 is given by 


n-1 


= 0 n + E (-l) m+1 2^ for odd n 


.. (I 1 . 5 a) 


m=1 


n-1 


j3 a* <j6 - 2 d> + E (-1 ) m 2 4 for even n 

n r n T o „ v ' 'm 

m= 1 


.. (II. 5b) 


If 6 ,p can be made to vanish by properly selecting phase 
shifts, the first-order deviation from complete inversion in S_^ 
will be nullified in Thus, it is possible to generate 

composite n pulse sequences , S^, etc. starting from S q for 

broadband spin inversion over large ranges of Aco^ and Ato • One of 
the composite inversion 7 T pulse thus achieved is 


S 1 * (") o (*) o (l0, 2o (n) 6o (*), 2o 
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IX. B . k Fixed Point Theory 

The propagator that corresponds to any rf pulse sequence 

can be represented by a point in the propagator space, which is 

a subset of Liouville space. For example, the propagator space 

for a quadrupolar spin l»1 nucleus is an eight dimensional space 

but if the quadrupole interaction vanishes then the dimension of 

the propagator space reduces to three. The purpose of this 
1 2 

method is to find a proper function (F) which generates the 

higher iterative sequences (U_^ + ^ = FU^) in euch a way that (i) 

the higher iterates , U^, etc. should be determined by the 

initial one (U q ) and (ii) for a given there should be only 

one U . „ . 

x+ 1 

When there is only one sequence given by the propagator U 
describing the desired response, that is, = U if U q = U and 


U = FU or U = F U 


.. (II. 6) 


where F is the n iteration of the function F, U is then said to 
be a fixed point of F. In general, there may be a set of propa- 
gators which give the desired results, 


{u< k >} = f{u^^J .. (II. 7) 

Then both and are said to be invariant sets of F. 

An iterative sequence with a stable fixed point (shown in 

Figure 11.2(a)) will generate broadband excitation sequences 

whereas F with unstable fixed point (shown in Figure 11.2(b)) 

1 2 

generates narrowband excitation sequences. A range of initial 
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Figure II. 2 (a) Points near a stable fixed point converge 

to the fixed point with successive iterations 
of F. (b) Points near an unstable fixed point 
diverge. The axes are along directions in a 
multi-dimensional space. 

pulse sequence propagators, arising from a range of experi- 
mental parameters, converges to the desired propagators upon 
iteration in the case of stable fixed point. Generally, a 
simple expression for F is not possible. For this, basins 
must be mapped with a computer program. These basin images 
define the dynamics of F, superstable fixed points. Basin 
images will have regions shaded according to the number of 
iterations required for an initial point in a given region 
to converge to a fixed point. Basin images can be used 
to guide the selection of an initial pulse sequence so as to 



get better composite pulses. An initial sequence is selected 
for which the locus of the initial points conforms to the basin 
of the desired fixed point. 

For example, the composite inversion ft pulse 300 120 is 

obtained by a computer search conducted over possible two pulse 

sequences with phases 0 and 180. Rotation in three dimensional 

space can be represented by matrices which are elements of a 

spherical subspace belonging to the group called S0(3)?^ The 

points in S0(3) space have been calculated for each possible 

sequence with the variation Aw from -1.4 go ^ to +1.4 q 0 . 

Similarly, ( 1 65) 0 ~( 1 65 ) 1 Q 5 -( 1 65 ) Q and many other sequences 

11 12 

have been generated by this procedure. * 

11 12 

The advantages of this approach * are the following: 

A unified set of concepts and the accompanying geometrical 
picture allow visualizations, comparisons of scheme and explan- 
ations of their essential features to be made. Dynamical concepts 
allow a clear treatment of schemes that will be normally difficult 
to consider. However, in higher dimensions, or in problems with 
arbitrary dimension, the basins are more difficult to be genera- 
ted by computer programs. 

XI. B . 5 Numerical Approach 

14 

This procedure removes restrictions on the pulse para- 
meters and allows the design of composite pulses with fewer 

components and shorter duration. The efficiency of a composite 

1 4 

pulse sequence can be defined by the function 
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F(x) = S E a pq [l(a p ,b q ) -S(x,a p ,b q )] 2 .. (II.8) 

where (ap,bq) is a point on the Aw /&) ^° -Aco^/o)^° plane at which 

the desired response is I ( ap , bq ) , a^ q determines the relative 

contribution of each point to F(x). Minimising F(x) with respect 

to the parameters x will give an optimum composite pulse sequence. 

This is carried out by a multi-dimensional search of the global 

minimum corresponding to the best set of parameter values. Some 

of the composite pulse sequences have been verified by this pro— 

1 4 

cedure and newer ones have been constructed. This method does 
not provide insight into the mechanism of error compensation. 
Though it generates efficient pulses, each iteration increases 
the number of pulses in the required sequence. 

Shaka^ ^ has applied a numerical non-line air optimization 
procedure arid obtained several phase alternated composite pulses 
for ultra broadband spin inversion in NMR spectroscopy. These 
phase alternated pulses have a phase of either 0° or l 80 ° only. 

In principle, it is possible to consider arbitrary phases also 
for the pulses. 

Generally, iterative schemes and pure numerical approaches 
demand larger computer time for generating composite pulses and 
they yield mostly longer duration pulse sequences, which will 
not be useful in NQR spectroscopy because of the shorter spin- 

spin relaxation time, T^. The composite n / 2 pulse derived from 

8 13 8 
recursive expansion 9 approach, namely, (90 ) o -(90)^q is not 

a better pulse sequence to compensate for efg inhomogeneity in 

NQR, Moreover, the mechanism by which these sequences compensate 
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for different effects so as to give broadband excitation is not 
yet clearly understood. 


II . C Magnus Expansion Approach 

In Chapter I, we have discussed the importance of Magnus 

approximation in magnetic resonance. Pulse sequences have been 

discovered for line narrowing, time-reversal and other such 
2 3 

experiments using the AHT which employs the Magnus approxi- 
1 8 

mation. The basic idea behind this theory is to get a pulse 
sequence, under which a particular line broadening interaction 
is averaged out. In AHT, Magnus expansion is used to solve the 
time dependent von Neumann equation defining the evolution of 
spin population under a desired rf pulse sequence. The solution 
obtained by Magnus expansion is expressed with a time independ- 
ent Hamiltonian called average Hamiltonian or effective Hamil- 

1 9 

tonian which is effective during the desired pulse sequence. 

AHT concept has been used to construct broadband decoupling 

^ ijr *|3 *|6 *] 7 

composite pulse sequences in NMR. » * » » In a similar 

manner, Magnus expansion has been used in the construction of 

11 , 20,21 

composite pulses for broadband excitation in NMR spectroscopy. 

The purpose of broadband excitation problem is to find a pulse 
sequence, whose net propagator will have a desired form and 
independent of a particular effect which causes the broadening 
of the frequency spectrum. 

In this approach, only a certain part of AHT is appli- 


cable. The important and common idea behind AHT and broadband 
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excitation problem is the use of 'toggling frame (or the rf 
interaction frame) to impart a time dependence to the internal 
Hamiltonian of the spin system. Spin evolution under a constant 
time independent effective Hamiltonian derived using the Magnus 
expansion is also common for both the cases. Otherwise, the 
requirement that U rf (t) = JL ^ (for a cyclic rf pulse sequence) 
of AHT is not at all a necessary condition for broadband excita- 
tion problem; rather, U r ^(t) has to be a rotation that creates 
transverse magnetization. Since composite pulses are not applied 
repetitively like multiple pulse sequences, the condition of 
stroboscopic observation (or the so called 'window' in AHT) ^ 
is not needed for the broadband excitation problem. 

Advantages of Magnus expansion approach are as follows. 

20 

This technique of finding out composite pulses is valid for a 
general spin I case. This method can be used for constructing 
sequences to overcome all kinds of interactions and at a time 
more than one effect can be considered. Here, the task of 
achieving broadband excitation sequences boils down to solving 
a set of analytical equations. This method treats composite 
pulses entirely rather than pulse-by-pulse analysis and, there- 
fore, the results are independent of the initial state of the 
system. Pulse sequences obtained by this procedure are equiva- 
lent to constant, pure rotations produced by single rf pulses. 

In addition to this, the mechanism by which different sequences 
are developed is easy to understand and to improve further. Hence, 
this method has a general applicability to many areas. Especially 
for NQR spectroscopy, it would be possible to derive shorter time 
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duration (<<T 2 ) sequences. Therefore, we have chosen this 
particular approach in the present work. 

Before concluding this section on different approaches 

employed in the design of NMR composite pulses, we would like 

to survey the following important recent literature. Counsell 
22 

et al. have proposed an analytical theory using the quaternion 

23 24 

formalism ’ for the analysis of several NMR composite pulses. 

2 5 

In 1987, Sanctuary and Cole showed that the analytical expre- 

26 

ssions generated by the multiple theory agree with other theore- 
tical approaches and experimental results. Alternative approaches 

27 28 2.Q 

like continuous phase modulation, amplitude modulation ’ ~ and 

29 30 

frequency switching * of rf pulses have been used to devise 

composite pulse sequences in the NMR literature. Recently, Zax 
31 

and Vega have used a perturbation treatment with Floquet 
formalism to derive broadband excitation pulses of arbitrary 
flip angle for two-level systems. In what follows, a general 
Magnus expansion technique for designing composite pulses in 
spin I NQR spectroscopy is discussed. 

II . D Magnus Expansion Approach for the Construction of 
Composite Pulses in NQR 

We consider an ensemble of physically equivalent spin I >l/2 
nuclei at the sites of a single crystal in the absence of Zeeman 
field. Defining a common quadrupole principal axis system ( QPAS ) 
for all physically equivalent spin systems and assuming that QPAS 
coincides with the laboratory frame, the total Hamiltonian, , of 
the system in QPAS is given by 
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^ = ‘Kq + # rf (t) + V 


.. (II. 9) 


o o 

where is the quadrupolar Hamiltonian" 5 and it is explicitly 
given in terms of Cartesian coordinate angular momentum opera- 
tors (I , I » I and I) in Chapter I (see Eqn. (i.l)). The 

X. y Z 

Hamiltonian in Eqn. (II.9) defines the interaction 

between the spins and an applied rf radiation of strength ( 0 ^ 
(rad/sec) with phase t) which is piecewise constant. Typically, 
the rf radiation is applied in the form of pulses that are ideally 
square, so that go ^ is piecewise constant having values either 0 
or only. It can be expressed in QPAS as 


^ r f(t) = — 2 ( 0 1 ( t ) c o s ( G)q t + 0(t)) I x 


.. (IX. 10) 


where is the quadrupole frequency. In writing Eqn. (II.IO), 
we have made the assumption that the x-axis of QPAS is aligned 
along the rf coil axis. Unlike in NMR, here both components of 
the linearly polarized rf radiation (i.e., two circularly polar- 
ized rf components with strength and rotating in opposite 

direction at a frequency (jq) interact with the spin magnetiza- 

33 

tion in zero magnetic field. 


In general, the Hamiltonian V in Eqn. (II. 9 ) may include 
any other interactions in the spin assembly. In NQR, these can 
be different effects that lead to the efg inhomogeneity (6q) 
which in turn create a quadrupole frequency distribution ( Ao^) » 
dipole-dipole interactions (( 1 )^)* spin-spin coupling, resonance 
offset and rf field inhomogeneity (Aoo^)* The mechanisms 



47 


32 

responsible for these interactions have been described in 

Chapter X. Here, resonance offset is the frequency difference 

between the applied rf radiation and the quadrupole frequency. 

Since Hamiltonians corresponding to the effects of resonance offset 

and efg inhomogeneity transform in the same way under unitary 

transformations, we denote both of them by Aco a in this thesis. 

The various terms of the Hamiltonians in Eqn. (II. 9 ) may then be 

arranged in an order, namely, #q» ft rf (t) V. Therefore, in QPAS, 

£ 

the time evolution of the spin system is largely brought about by 
the quadrupole Hamiltonian, Ttq which is already well-characterized 
and therefore contributes little new information. But usually in 
NQR one is interested in studying the effect of rf radiation on 
the spin system. So, it is necessary to suppress ftq. Furthermore, 
in QPAS, ft r ^(t) is itself time-dependent and it would be highly 
cumbersome to carry out the theoretical analysis of NQR experi- 
ments. These difficulties can be overcome by a transformation of 

ft to a new frame called 'quadrupole interaction frame 9 (QIF) 

34 

deFined by the unitary operator U * This Frame is somewhat 

3 5 

similar to that of the 'rotating Frame* used in NMR ; however, 

QIF cannot be visualized by a simple geometrical picture* 

Now in QIF, becomes 

* = U Q * ^rf (t) + V •• (11.11) 

QIF is represented by a tilde (*») over the Hamiltonians, 
where Uq = exp(-iT^t) 


• ■ • 


(XI. 12 ) 
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Derivation of Eqn. (XI.1l) from Eqn. (II. 9 ) can be easily under- 
stood with the help of Appendix C. In QIF, the transformed rf 

r*> 

Hamiltonian *rf( t) is time dependent and it possesses rapidly 

oscillating terms such as cos(nu^t) and sin(na)qt). Effects 

corresponding to these high frequency terms are similar to the 

Bloch-Siegert effect in NMR - ^ and in NQR, this has not been 

extensively studied because of its limited chemical applications. 

37 

However, in 1984, Furman reported the first observation of Bloch- 
Siegert effect in NQR. To a first order approximation these high 
frequency terms can be dropped, and the resultant rf Hamiltonian 
is time independent in QIF* Furthermore, the main inter- 
action may also be suppressed, and Eqn* (ll.1l) may be written 
as 

# = * r f + V •• (11.13) 

The procedure of truncation of high frequency terms will be 
explicitly shown in later chapters for spin 1=1 and 3 /2 cases. 

But it can be generally said that truncation means rejection of 
all those terms whose operators depend on time in such a way that 
their time average vanishes. 

The concept of QXF is more than a mere mathematical trick. 

In NQR experiments which employ phase sensitive detection, the 
observed experimental results correspond to the theoretical response 
of the spin system in QIF. We can then say that usually NQR 
experiments are performed in QIF. Now in QIF, the time evolu- 

O O 

tion of the spin system is described by the von Neuman equation -3 
as 
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= -i[*. ?(0)] .. (II. II.) 

where f(o) = U Q -1 P(o) U Q » 7£ Q .. (II. 15 ) 

in 

Derivation of Eqn. (IX. 15 ) for a general case is given^ Appendix C. 
The reduced density matrix P( 0 ) has been derived in Chapter I. 

The solution of Eqn. (II. 1 4) is given by the Dyson expression as 


U(t) = T exp [-i 




o 


X(t') d t '] > 


.. (11.16) 


where T is the Dyson time- ordering operator (see Chapter l). The 
density matrix defining the spins at time ' t' in QIF is given by 


P(t) = U(t) P(o) u(t)- 1 


.. (11.17) 


where U(t) defines the evolution of the spin system under the 

>V A/ /S/ 

simultaneous effects of *H and V. Here, ft V* Therefore, 

r 1 ri 

the main time dependence of P ( t ) comes from Since we would 

like to design composite pulses that reduce the effects of the 

A/ 

interaction V, it is necessary to understand further the evolution 

/s/ ^ 

of the system under V alone. The important parts of V can be 
made apparent by a transformation to a new frame called the 
'toggling frame' ( i . e . , the rf interaction frame). The concept 
of this new frame is similar to the one that has been implemented 

19 

in AHT. The frame is so called because it toggles with respect to 
QIF, according to the rhythm of rf pulses in the applied compo- 
site pulse sequence. Toggling frame and QIF coincide in the 
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absence of rf pulses in the present problem, but in the case of 
AHT the toggling frame and QIF coincide only at windows. In the 
toggling frame, the net Hamiltonian is expressed as 

^ * u r ;V)a u rf (t> = v(t) .. ( 11 . 18 )' 

Throughout this thesis we represent the toggling frame by double 

tilde* (#). It should be noted that the transformed Hamiltonian, 

a 

namely, V(t), is time dependent. Alternatively, we may say that 
the internal Hamiltonian is deliberately made time dependent. 

Now, the equation of motion becomes (Appendix C) 

= -i[ v(t), F(o)] ..( 11 . 19 ) 

where P(o) = P(o) .. (ll.20) 


The solution to this differential equation is derived to be 

%. & 

U v (t) = T exp [ -i / V(t« ) dt' ] . . (11.21 ) 

o 

In QIF, the complete evolution operator, i . e . , Eqn . (II. 1 6 ) , 
becomes 


U(t) = u rf (t) U v (t) . 


.. ( 11 . 22 ) 


One can notice the importance of the toggling frame in separat- 
ing the simultaneous time evolution under 71^ and V* From 
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Eqn. (II.19) » it can be stated that in the absence of internal 

interactions represented by the Hamiltonian V, the spin system 

in the toggling frame denoted by ^(t) remains stationary (or it 

a 

is said to be a constant of the motion). Since V(t) is time 
dependent, U v (t) is time dependent. Therefore the problem now 
is to solve a time dependent differential equation, i.e., Eqn. 
(II. 19 ). As has been discussed in Chapter I, we use the Magnus 
approximation for this purpose in a manner similar to that of AHT, 
to write 


V = 


. v(°> 


(XX. 23) 


In Chapter I, the different terms in the above mentioned Magnus 
expansion have been separately explained (Eqns . (I.25) and 
( 1 . 26 )). These terms can be derived for various spin values and 
they can be expressed in terms of rf pulse phases (0 n ) and flip 
angles ($ n )» Eor example, considering a composite pulse contain- 
ing N-pulses, it can be completely defined by 2N parameters, 
namely, the N-f lip angles and N-phases. Schematically, a N pulse 
composite pulse has been represented in Figure II. 3. Now 
can be written in terms of these 2N variables (or unknowns). At 
this point, if we can make V(t) = 0 by adjusting the 2N unknowns 
in such a way that goes to zero, then the propagator 

U (t) = H . This means that the spin populations would not feel 
the effect of different internal interactions. Therefore, the 
applied composite pulse sequence excites the spin system uniformly 

irrespective of the strength of the interactioi^j¥^^j^:|s 

I I.T., KANPU". 




Figure II. 3 Schematic representation of a composite pulse 

consisting of N rf pulses with individual time 

durations denoted as t, , t_, t„, ... t. T . 

l 2 3 N 


achieved by selecting an mth— order composite pulse which satis— 

39 

fies the following conditions: 

(i) all the Magnus terms should be equal to zero, i . e . , 


V 


(n) 


0 for 0 ^ n ^ m 


.. ( 11 . 24 ) 


As m increases, the degree of compensation increases. So, 
an mth-order composite pulse will produce better broadband 
excitation than an (m-1 ) th- order composite pulse. However, 
commonly higher order sequences will be of longer time 
duration. 

(ii) The composite pulse that satisfies the above condition should 

behave like a single rf pulse. For example, a composite 7T / 2 

pulse is defined by the maximum spin magnetization in the 

39 

transverse plane of QPAS i.e., 

Tr [M , U „(t) P(0) U ^(t)" 1 ] = <(l } 

L x’ rf' ’ ' ' / rf' ' J N x ' maximum 


• ( 11 . 25 ) 
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where M is the observable operator in QIF. The derivation of M 

x 

is given in Appendix D for I = 1 and 3/2 cases. Similarly, a 
composite inversion u pulse is given by the minimum expectation 
value of I z and by zero transverse spin magnetization. It can be 
expressed as 


Tr [M z , U rf (t) P(0) U rf (t) _1 ] 



minimum 



m aximum 


<* x > = (l y > = 0 . 


.. ( 11 . 26 ) 


Therefore, an mth-order composite pulse consisting of N rf pulses 

can be derived by solving (m+l) equations, namely, Eqn . (11.24) 

and Eqn. (IX . 25 ) or Eqn. (11.26) simultaneously so as to get 2 N 

variables which define the complete sequence. The efficiency of 

these composite pulses can be examined by evaluating the response 

of the spin system (refer Appendix D) as a function of a parameter 

which defines the interaction V. For example, for a composite 

7l / 2 pulse, the signal anplitude, , should be evaluated by the 

39 

following expression i 


Tr[M x , P(t)] 


.. (11.27) 


Similarly, the extent of spin inversion by a composite % pulse 
is given by 


¥ = 
4*2 


•Tr[ M z , P(t )] 

Tr[ I 2 ] 


. (11.28) 
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Using this approach composite pulses with any arbitrary net flip 
angle can be designed. But as far as NQR spectroscopy is concerned, 
only composite ft/2 and composite ft pulses are important. 

The main goal of this thesis is to design shorter time 
duration composite pulses (ft / 2 and ft) with broadband 

and 6)^ '+. Aco-j ) excitation performance that would be appli- 
cable in any kind of FT NQR spectroscopy. In the present work, we 
t ourselves to zeroth order composite pulses because usually 
higher order sequences are of longer time duration and therefore 
they may not be of much importance in NQR spectroscopy. In this 
approach, composite pulses compensating simultaneously, for 
different effects can be derived. But for the sake of simplicity 

we have considered different effects separately. This procedure 

39 

applied for spin 1 and spin 3/2 cases have been presented in 
Chapters III and IV, respectively. 

Sufflm ary 

In this chapter, the various procedures for designing 
composite pulse sequences reported in the recent NMR literature 
have been briefly reviewed in the light of their applicability/ 
adaptation to NQR spectroscopy. The merits of using the Magnus 
expansion approach have been pointed out. The purpose of the 
Magnus expansion employed in both AHT and broadband excitation 
problems has been clearly brought out. A general theoretical 
treatment of designing composite pulses for an arbitrary spin 
I ( ^. 1 ) nuclei for use in NQR spectroscopy has been developed. 
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CHAPTER III 


DERIVATION OF NQR COMPOSITE PULSES 
USING THE MAGNUS EXPANSION APPROACH 
FOR SPIN I = 1 NUCLEI IN SINGLE 
CRYSTAL SPECIMENS* 

In the previous chapter, we have developed a general theore- 
tical procedure based on Magnus expansion approach for the cons- 
truction of composite pulses in NQR spectroscopy. This chapter 
deals with the procedure for deriving composite n /2 and composite 
7T pulses that can compensate for efg and rf field inhomogeneities 
in pure NQR spectroscopy of single crystals containing spin I = 1 
nuclei . 

This chapter is organized els follows* Section III. A presents 
the Hamiltonians that represent the various interactions present in 
the model system assumed for our theoretical study and they have 
been expressed in terms of fictitious spin-l/2 operators. Cons- 
truction of composite pulses that overcome the problems of efg 
and rf field inhoraogenei ties through the Magnus expansion technique 
and the examination of composite pulses have been presented in 
Section III.B. Section III.C covers the derivation of phase alter- 
nating composite pulse sequences (PACPS) and their performance with 
respect to offset or efg inhomogeneity effects. 

♦Part of the material in this chapter has been published in 
J. Molec. Struct. 122, 333 ( 1 9S9 ) . 




6i 

XIX. A Hamiltonian of the Spin System 

We consider an assembly of non-interacting, physically equi- 
valent spin 1=1 quadrupole nuclei in the absence of a Zeeman 
field. Only the following interactions are taken into account in 
the model system assumed here, quadrupole interaction,^,^, various 
interactions as defined in Chapter I which lead to efg inhomo- 
geneity, rf interaction, the inhomogeneity in the applied rf field. 
Since only physically equivalent spins are considered, we can 
define a common quadrupole principal axis system (QPAS)^ for the 
total spin system. For simplicity, coincidence of QPAS and labo- 
ratory frames is assumed. Although it would seem to be a crude 
approximation to consider the spins as an ensemble of isolated non- 
interacting single spins, it is a good starting point for cons- 
tructing composite pulses for NQR spectroscopy. In what follows, 
spin Hamiltonians are explained in detail. 

XII. A. 1 Quadrupole Hamiltonian 

A general form of the quadrupole Hamiltonian, for a 

spin I nucleus in terms of cartesian angular momentum operators 

is given in Chapter I (vide Eqn. (i.l)). For spin 1=1 nucleus, 

"fi , , can be written as^ 

Q 

\ ( 3 i z 2 -i 2 +r? (i x 2 -i y 2 )) .. (m.i) 

In the present spin 1 theoretical analysis, we use the fictitious 

2 

spin-l/2 operator formalism of Vega and Pines as a basic mathe- 
matical tool. This formalism is explained in Appendix A. Using 
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fictitious spin-1 /2 operators, is expressed in QPAS by 



2 C\ 

_ £ ( x _ x 

2 v X3 Y3 

I Z3^ * 


.. (III. 2 ) 

3 

Following Cantor and Waugh, 

we write 




= “p I p3 + 

W 5 

p = X, Y, Z 


.. (III. 3 ) 

where 

a) 

“ 4 

(n +3) ; 

a _ 

X 4 

(77 - 1 ) 

.. (III. 4a) 


© qQ. 
toy “ 4 

( n - 3) ; 

</ 

W Y “ 4 

(r? + 1 ) 

.. (III. 4b) 


2 

-e aQ 

“Z = 

n i 

</ _ efaii 

Z 2 


.. (III. 4c) 


The energy level diagrams along with, energies, transition frequen- 
cies and wave functions for a quadrupole nucleus of spin 1=1 in 
zero magnetic field are given in Figure 111.1 for r)= 0 and T) £ 0 
cases, to x’ toy W Z are t * ie t ^ iree pure quadrupole resonance 

frequencies indicated in Figure III . 1 . In FT NQR, at a time one 

3 

deals with only one transition, we consider here the transition 
with frequency to^ (see Figure III.l) and let to^. = to^ ; to^. = . 

Now Eqn. (ill. 3) becomes 

= to q I X 3 + toq ^X4 * ** (l-£-£*5) 

III. A. 2 EFG Iahomogeueitv Hamiltonian 

Various mechanisms which lead to efg inhomogeneity and 
broadening of NQR lines have been discussed in Chapter I. We 



lift [l+l> + l-l>] 


I0> 

(a) n=o (b) n#o 

Figure III .1 Energy levels and wave functions for a quadru- 
pole nucleus of spin 1=1 in the absence of a 
Zeeman field, (a ) T] = 0 and (b) T) / 0 cases. 

represent 5 q to be the spread in efg and A( 0 p to be the frequency 
distribution or resonance offset parameter. Now, the efg inhomo- 
geneity Hamiltonian, V , in QPAS is given by 

V 1 = ^ w p I X3 + A(0p I x4 5 P = X, Y, Z . .. (ill. 6) 

Substitution of A q for q in W and CO gives expressions for Aw 
jr -P P P 

and Ao) p » respectively. As it has been explained in Chapter II, 

4 

QIP is tlie reference frame in NQR spectroscopy and in QIF, tlie 
transformed efg iahoniogeneity Hamiltonian is 

^1 = V 1 v i = v i 5 r U Q* V J = 0 • •• (hi. 7) 

A more explicit expression for is given in Eqn. II . 12 ). 

III. A . 3 RF Hamiltonian 

Writing the rf Hamiltonian (Eqn. (il.io)) in terms of 
fictitious spin-l/2 operators in QPAS we get, 
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K rf = “ 4t0 1 cos ( (0 p t +(jf>(t)) I 


Pi 


.. (III. 8) 


where (rad/sec), y is the gyromagnetic ratio of the 

nucleus and is the rf field magnitude in units of Gauss. 

In QIF, 


^ rf = U Q ft rf U Q 

= - 4( °1 co » (“p* +^(t)) exp[(i# Q t) I p1 exp ( r i Tlqt) ] 


where exp(-i# Q t) = exp(-iu) p I^t) exp(-io)' I^t) . 


.. (III. 9 ) 


. (III. 10) 


In writing this expression, we have used the property £l ,^,1 ^1 == 0 . 

To evaluate the transformation property of I under 9 L , we use the 

P * '•% 

commutation relationships given in Appendix A. Applying the nested 
commutation relationship (Eqn. (l. 2 l)), becomes 


/V 

\f = ~ 4w i cos (“ p t + #(*))[ I p1 cos (0 p t -I p2 sin Wp t] . .. (Ill.li) 


Selecting only one frequency i . e . , w = q 


% 

rf 


{ 2 sin(2co^t) 

1^1 [c°S G)qt - — ] cos <jt>(t) 


_I X2 ^ 


sin(2o), t ) 

JsL 


+ sxn 


in 2 oi 1 s in <£( *)}. .. (III. 12 ) 


Now, truncating high frequency terms, we obtain 


Cs) 


U rf = -2 Wl [ I X1 cos<j^(t) - I y? sin<^(t)] 


.. (III. 13) 
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XXX. A. 4 RF Field Inhomogeneitv Hamiltonian 

Denoting = f AH^ to be the inhomogeneity in rf field 

t 

strength, the rf field inhomogeneity Hamiltonian, V^, in QPAS 
can be written as 

V 2 = -4 cos C® p * + ^(t)] I p1 r • (HI. 14) 

In QIF, 

V 2 = -2 Aco 1 [ I X1 cos 0( t ) — *^X 2 ^ ^$( ^) 1 * • • (lXX*15) 

The totaX Hamiltonian of the spin system in QXF is (refer 
Appendix C) 

^ = ^ rf + V, + V 2 . .. (III. 16 ) 

After thus representing spin Hamiltonians in terras of fictitious 

2 

spin -l/2 operators in the observation frame (QIF ) , in the next 

section, we shall design composite pulses based on Magnus expansion 

5 6 

approximation for spin 1 NQR spectroscopy. 

III.B Composite Pulses for Spin 1 NQR 

The importance of composite pulses in NQR spectroscopy has 

7-9 

been highlighted in Chapter I. Various routes of designing 
composite pulses for NQR have been analyzed in Chapter II. In 
the area of NMR one already has, at present, composite pulses 
starting with a minimum of two pulses going upto several pulses, 
typically nine pulses. 


But in NQR, the application of such 
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longer time duration composite pulses is limited by shorter • 
Therefore, we confine ourselves to construct two and three pulse 
composite pulse sequences* Ve denote the composite pulse sequence 
as (e,)^ 1 -(e 2 ) 9 i 2 -(e 3 ) ^, where represents nth-rf pulse in 

the sequence with flip angle 8 n and phase 0 n > 

III.B.1 Composite Pulses for the Compensation of efg 
Inhomo genei tv 

Since the inclusion of higher-order Magnus terms usually 

9 

results in a longer time duration composite pulse, we consider 
only the main term, namely, zeroth-order term in the Magnus 
series. In the next subsection, we present zeroth-order compo- 
site pulses for overcoming efg inhomogeneity. 

III.B.l(i) Design of Zeroth-Order Composite Pulses 
In QIF, the net Hamiltonian is given by 

ft = + V 1 . .. (11.17) 

Using the idea of % toggling frame * , 1 ° we arrive at (Chapter II) 

U = U At) U (t) . • . (H.18) 

1 

The rf propagator corresponding to different time intervals for a 
three pulse sequence is given below; 



(t) =< 


^exp (2i(i)^ t ,| ) » y • • • 

exp exp(2i(*)^ ) exp ( — ^7 

e xp ( i t • • • 

exp(i^ 3 I X3 ) exp(2i«,t 3 X Xl ) exp^i^I^)' 

exp ( i ^ 2 i X 3 ) ex P ( 1 ^2 Z X1 ) 
exp(-i^I X3 ) exp ( 10^ x xi) 


o < t ^ — — 

' 1" 2 CO, 


06 *2^ fe 


9, 

-r-*- 

3 2«*X| 


where t denotes the nth~rf pulse width. 

In writing the above equation, we have made use of the 
following points: 

(a) phase of the first rf pulse in the composite pulse sequence 

is assumed to be a reference rf phase of the sequence and 
it has been set to zero i • e . , 0 ^ = 0. Alternatively, one 

can assume the phase of any one of the rf pulses to be the 
ref erence -phase . 

(b) In the presence of non-zero phase {0 ^ 0 ) , the rf Haxnil- 

rs» 

Ionian, ^rf (Eqn.lH.13 ) contains two operators, namely, 1^ 

and Ij£ 2 * Since [ 1^ , I^J £ 0, we make a transformation to 

a new frame defined by the unitary operator, U =exp(i{$ I Y )♦ 

' n aj 

/V/ 

In the new frame, the transformed rf Hamiltonian is *}{ _ = 

rf 

-2co^I x ^. During each rf pulse, we take the spin system to 
this new frame which depend on the phase of the correspond- 
ing rf pulse and we bring back the spin system to QIF after 
the effect of % T f interaction. This procedure simplf ies 
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the evolution of the spin dynamics under the action of U _(t). 

rf 

Note that when the rf phase (jf> = 0, the new frame anH qif 
coincide . 

(c) The flip angle of an rf pulse is stated as 

0 n = effective rf field strength multiplied by the pulse width 
= 2W 1 t n .. (III. 20) 

is commonly given in radians. 

The net evolution of the spin system in the toggling frame 
is given by 



L 

exp[-i J dt' V ( t ’ ) J 


.. (III. 21) 


V^(t) is the complete Hamiltonian or the transformed efg inhomo— 

geneity Hamiltonian in the toggling frame and clearly, it is time 
£» , , 

dependent. V^(t) has been derived using the rf propagator in 
Eqn. (ill. 19) to be^ 


V 1 (t) = Aaiq [a(t)l Xl + b(t)l X2 + C(t)l X3 ] + Ao^xi,. 

.. (III. 22) 

Since [ l^t ^(o)l = 0* the last term in the above equation does 

not effect the spin system. Therefore, we write 


A 

v -,(t) = Aou g [(a(t)l Xl + b(t)l X2 + c(t)l X3 j] .. (III. 23) 

where the coefficients a( t ) , b(t) and c(t) are time dependent and 
they can be given in terms of and 0 n at different time inter- 


vals of the pulse sequence as 
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f, 


0 , 


. . oft £ 


ft 


|-sin(2co 1 t 2 )sin9i 2 , 

a(t) =^ sin (26) 1 t 3 ) [| sin(2^ 2 ) cos0 3 (l-cos0 2 ) 


1 ~ 20 > t 

&> 

.. t< - 2 - 
2 2 0 )^ 


-sin$ 3 ( cos 0 2 + sin 2 ^ 2 cos^)]^. . 0 <£. 

- 00 ( 2 ^ t 3)sin0 2 sin$ 2 


7 


^■-3111(2^ t 1 ) , 

-cos(2co 1 t 2 )sin9-, - sin^O).! t 2 ) cos$ 2 cos^, 


.. ( 111 . 24 ) 


6 . 

0 £t. 

1 2 ^ 

9 2 

0 -* 2 - ST 


°OS0 1 [-cos *2 sin & 2 cos(2o> 1 t 3 ) 

O 

0 2 cos ^2 sin(2o> 1 t 3 ) cos^ 




-COS 


b (t) =< 

- 2 sin ( 2 0?) cos $ 2 sin(2co 1 1 3 ) sin0 3 


'2 

- sin 0 2 sin(2< J ) 1 t 3 ) cos0 3 

+ ~ sin <4 2 ) sin(2oo 1 t 3 )sin^] 

^-^sin9^[ cos0 2 cos (2U).j t 3 ) - sin& 2 sin( 2(o ^ t 3 ) . 

•cos(0 2 +0 3 )] 


6 ' 


•• °^ 3 - 20 > 


(III. 25) 



and 
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f 


cos (2GX| t 1 ) 


9 


cos ( 2 V 2 > cose r sin(2 Ul t 2 ) cos (j) z sin^, 



O^t 
0 £t 



1“ 2Q. 


2 2 ( 0 . 


1 


o 



.. (111.26) 


As it lias been discussed in Chapter XI, we apply the Magnus 
expansion approximation-" to get an effective time independent efg 

$5 

Hamiltonian in the toggling frame. V^(t) can be expanded into an 

infinite series 




+ V. 


0) 


+ V. 


( 2 ) 


.. (III. 27) 


i Jl ) 

where is the nth-order term in the Magnus expansion. The 

explicit expression for the zeroth- and first-order terms are 

given in Chapter I (Eqns. (1.25) and ( 1 . 26 )) and higher order 

5 11 

terms are al so available in literature. ’ Now, according to the 
Magnus expansion approach^ (Chapter II) , a zeroth-order composite 
pulse is one which has = 0. In this zeroth-order approxi- 

mation, we truncate the higher order terms whose contribution 
would be negligibly small when compared with the main term V^ 0 ^. 
for a three pulse composite pulse sequence can be obtained 
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by integrating equations (111.24), (ill. 25 ) and (ill. 26 ) with 
respective time intervals of the sequence. on integration 

or, alternatively, ' time— averaged* in the language of AHT, ^ ^ turns 

6 

out to be 


(0) 


_& 
2 W 1 


.a I 


XI 


+ b I 


X2 


+ c I X3 ] 


. . (III. 28) 


Now the coefficients are time independent constants for a given 
rf pulse sequence; otherwise, they are dependent on 0^ and for 
an arbitrary pulse sequence. Explicit expressions for a, b and c 
have been derived and are presented below. 


& 


sinffx^ ( cos 0 2 -sin0 o sin 0^-1 ) 

+ ( cos0^~ 1 ) [cos + — cos© 2 sin2$ 2 

+ cos0 2 sin^sin 2 ^> 2 - — cos0^ sin2$> 2 ] 


cos 


<h 


b = cos 0^-1 -sin 0^ (sin 0 2 + cos© 2 sin6 3 ) 

+ cos co s(^ 2 ( cos © 2 -1-sin© 2 sin^) 

+ (cosD^-l ) [-sin^ sin0 2 cos (M»> 

2 2 

+ cos0^cos^)j{cosB 2 cos 0 2 + sin <p^) 
+ t; cos© 1 sin2^> 2 sin^ 3 (cos© 2 +l)j 

and 

c = cos0 1 (sin0 2 +cos0 2 sin© 3 ) + sin© 1 
+ sin0^ cos0 2 (cos0 2 -sin0 2 sin© 3 -1 ) 

+ (cosO^-l )[sin0 1 sin^ 2 sin(0 2 -0 3 ) 

+ sin0j cos0 2 cos^ 2 cos(0 2 -<^) 

+ cos ©jSinBgCos (0 2 -$ 3 )] 


(III. 29) 


(III. 30) 


(III. 31) 



72 

For a zero th-order composite pulse sequence with = 0, a, b, 

and c should vanish identically . We have optimized 0^ and for 
a three-pulse sequence to get = o and we have arrived at a 

composite inversion jr pulse^ namely, (90)^- ( 270) (90)q* However, 

it is a difficult task to get sequences with ^ ^ = 0 by varying 
0n and 'fo reduce this complexity, we have relaxed the condi- 
tion from = 0 to 

[v 1 (0) ,e(o)] =o .. (hi. 32) 


According to tliis new condition^ ^ ^ does not necessarily have 

to be zero but if V ^ or U (t) commutes with the thermal equi- 

1 

librium density matrix, P(o), then the spin system will not evolve 

■t 

under the efg inhomogeneity Hamiltonian. Hence, a uniform broad- 
band excitation would be achieved. This condition demands only a 
and b to be zero independently. However, it is felt, even at this 
stage, it is not easy to derive many better sequences for our 
purpose. So, we have obtained sequences with even smaller terms 
of a and b as composite pulses.^ Thus, we could arrive at the 
following sequences! (90)q- ( 225) ^ gQ- (31 5 )q has been found to be 
a composite inversion 71 pulse. (90) 0 -(300)^ 0 and (385) Q -(320) 1 g Q - 
(25) q are derived for composite 7l/ 2 pulses. ^ A summary of zeroth- 
order terms in the Magnus expansion associated with these pulse 
sequences is presented in Table III.1. In what follows, we test 
the behaviour of the aboveraentioned composite pulses. 
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Table III.1 Zero th-Order Magnus Expansion term for different NQR 
Composite Pulse Sequences that Compensate efg Inhomo- 
geneity in the Spin I = 1 Case 


Composite pulse Zero th-order Magnus expansion term, 

(a) Composite Tt / 2 pulses ; 

i) (90) 0 ~(300) 90 Ag^ [ -0.0734 - 0.0197 I X2 +0.1469 I X3 ] 

ii) (3®5 ) q- ( 320 ) ^ qq— ( 25) q AcOqC — 0 . 0026 I X2 + 0 . 0026 1^^ 

(b) Composite K pulses : 

i) (90) 0 -(270) 90 -(90) 0 0 

ii) (90) 0 -(225) l80 -(3l3) 0 Awq [ -0.0533 I X2 +0.3105 I X3 ] 
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III* B. 1 (ii ) Evaluation o f the Performance of Composite Pulses 
•£& — — fun ction of EFG Inhomo gene i tv 

It is necessary to know the extent of compensation against 
efg inhomo gene i t y for a given composite pulse sequence* It can 
be obtained by evaluating the time evolution of the thermal 
equilibrium density matrix ^(o ) , under the simultaneous effect of 
the rf pulse in the sequence and the efg inhomogeneity # ^ The 
complete Hamiltonian of the spin ensemble in QIF can be written as 

*= oos ( 4 n- I X2 sin ^J + X X 3 + I X4 •• ( II:t -33) 

since [l xZt , I x± ] = [l x4 , ?(o)] = 0 ; i = 1, 2, 3 . .. (111.34) 

We drop the 1^ term which is not effective in the spin evolution* 

Now , H is given by 

2s *"2co^| 1, cos 5^**“ 1 x.2 ain^^ J + Agoq * • • (ill *35) 

Since the two terms in the above equation do not commute, to 

simplify our calculations we diagonalise * This is done in 

6 

two steps. First, a transformation by a unitary operator, 
namely , «p(-i <f> n 1^), fallowed by a transformation with 
exp(-i)31 X9 ) , where p = tan -1 (2^/Aw^) . Now, the diagonalised 

Hamiltonian is 


d 

% 


exp(-i/31 X 2) exp(-i5& n 
[ 20) 1 sinj3 +Aoo ci cos/3]l 


X X 3 ) % ex P(i0 n 

_ °n J X 3 
X3 ~ t n 



) exp( ij3 


~X2 


) 


f 


(III. 36 ) 



where 
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a n = sinB + c °s 8 ) .. (111.37) 

is the effective flip angle of an rf pulse in the presence of efg 
inhomo gene i ty » In deriving rfi , we have used the commutation 
relationships given in Appendix A* Essentially, the Hamiltonian 

fsj 

'Jl is transferred to a diagonal frame and hack to QIF during each 
rf pulse in the sequence. The net propagator given by Eqn. (II. 16 ) 
becomes, for an rf pulse in the sequence, 

U n (-fc) = exp(i56 n I X3 ) exp(l8I X2 ) exp(-ia n I X3 ) exp(-i£l X2 ) expC-i^I^) 

.. (III. 38 ) 

The complete propagator representing a N-pulse composite pulse 
sequence is given by the product of individual rf pulse propaga- 
tors time-ordered by the Dyson operator^ els 


U(t) 


u», u„, 

N N_1 


.U 


n 


* U 3 U 2 U 1 


.. (111.39) 


For a three-pulse 
work of the type 


sequence that has been considered in the present 

(e,) 9 s i -(| ) S 4 2 -( e 3 V 3 ’ “ e write 


U ( t ) ss e xp ( ® -^P ( i jS'^X 2 ^ oxp ( ) exp (—ip 

exp r 1(^3) i X3 ] exp(i0I X2 ) exp(-ia 2 I X3 ) 
exp(-i0I X2 ) exp [ i( 0 1 “^ 2 ) I x 3 1 ex P(^ I X2^ 

exp(i ai I X3 ) exp(-i0I X2 )exp(-i<^ 1 I X3 ) .. (XXI.4o) 


Now using this time evolution propagator, the expression for the 

rJ 

density matrix at the end of the composite pulse sequence ?(t) 
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can be evaluated (Eqn. I1.1 7 ). Therefore, the performance of a 
composite n/2 pulse is given by the intensity of the NQR signal 
as (Appendix D) (Eqn. II. 27 ), 

1 ' Mi/ 2 ] 

” cos ^[*>,»inoj + bjcosaj] 

-sin^j[-b^si n 0 + b^cos0cosa^ - b 2 cos£sina„] . . (XII.4l) 


where I^ 2 is the observable operator in QIF and this has been 
derived and explained in Appendix D. 

Similarly, the performance of a composite inversion 71 pulse 
is given by (Eqn. (11.28)) 


¥ 


2 


-Tr[l X3 ,f(t)] 

Tr ^- I X3^ 


b^cosjS + b.j sinjS cosa^ - bgSinjSsinctj 


.. (11.42) 


where 

b.j = cosjScos cos /3cosC^-a^cosP sinOL^-a^sinP ] 

+ cos/ 3 sin($2“03 )[ a^ sinc^ + a 2 cosc^] 

+ sin/3 [ a^ sin|6 cosa 2 -a 2 sinjS sina 2 +a^cos^ ] .. (II. 43 ) 

b 2 « cos <M,>[ a^sina 2 + a^cosa^] 

+ sin cosj 3 cosa 2 + a 2 cos8sina2 + arsing] 

.. (11.44) 
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sin/3 c °s (^“^ ) [ -a i cosficosa.^ + a 2 cos 3 sina 2 + a^sinjB] 
-5in 0 sin(9i. 2 -9^)[a 1 sina 2 + a 2 cosa 2 ] 

+ cos3[a lS in0cosa 2 - a 2 sin3sina 2 + a^cos^] .. (ill. 45 ) 

a 1 = 2 sin 2 ^[ sin O 1 sin( 0 1 - 0 2 ) + cosg cos( 0 1 -<?^) (cosa.j -1 ) ] 

+ sin 2 |g[cos 2 j 3 + cos^sin 2 ^] . . (m. 46 ) 

sin^^sin (Xj cos(^-^) + cosjS sin (0^ - ) ( 1 - cos a,)] 

.. (III. 47 ) 

-sin 2 |3 [ sina.j sin (0,-4) + CO s$cos ( 0 - 0 ) (cos ^-1 )] 

+ cosj 3 f cos 3 + cos^sinl . . (m. 48 ) 

Itie performance of composite 7T/2 pulses is compared with 
that of a single Tr/2 pulse with respect to their ability to 
compensate for efg inhomogeneity. In Figure III . 2 the signal 
response is plotted as a function of A(Oq which is expressed in 
G ) 1 units. In Figure III. 3 , we compare the inversion efficiency 
of composite inversion jx pulses with that of a single Tl pulse as 
a function of the efg inhomogeneity parameter. 

I11.B.2 Composite Pulses for the Compensation of RF Field 
Inhomo gene i ty 

III.B,2(±) Design of Zeroth-Qrder Composite Pulses 
The net Hamiltonian in QIF is given by 

AJ ~ fJ 

+ V 2 


and 


a. 


3 


(III. 49) 




Figure III. 2 The signal magnitude as a function of Acoq/co i for 
a single tx/ 2 pulse (squares), the zeroth-order 
composite pulses (90) 0 -(300) 9£) (circles) and 
(385) 0 -(320) l80 -(25) 0 (triangles). 

Figure III. 3 The extent of the inversion of spin population, 
as a function of A Tor a single n pulse 
(squares), the zeroth-order composite pulses 
(90) 0 -(270) 90 -(90) Q (triangles) and 
(90) 0 -(225) l80 -(315) 0 (circles). 






19 


In the toggling frame defined by U (t) (Eqn. (ill. 19 )), we have 




%= V p (t) 


.. (III. 50 ) 


The net evolution in QIF is written as 


u(t) = u ,(t) u v (t) 


.. (III. 51 ) 


where U y (t) = T exp[-i J dt* V (t*)] 
2 o - 


.. (III. 52) 


g(t) has been derived using the commutation relations given in 

6 

Appendix A as 


V 2 (t) = -2Aw 1 [a»(t) I X1 + b ' ( t ) I X2 + c'(t) I X3 ] 


. (111.53) 


The coefficients in the above equation at respective time intervals 

are given below: 


l ' ( t ) =<cos ^ 2 

l( A 1 cos^ + A 2 sin^ 2 ) 


b' (t) 


-COS0.J sin^ 2 

cos A 2 cos$ 2 -A^in (J> 2 J-A^sinf^ 


0 4 t £ -57- 

Bo 


oit 3^ 

.. ( 111 . 54 ) 

0 

1 2 G) 1 

®2 

0 s£r 

.. (111.55) 
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c 1 ( t ) =\ -si sin^ 

sin h + A 2 cos^ 2 ]sin0 1 + A^cos&i 

where A^ = cos 

Ag = cos f^sin^g-^ ) 
and A = sin0 2 sin <4 Z -L). 


. . o 6 ^— 

1 2 to. 

9 n ' 

.. O ^t * 

2 2 (o .j 

. . 0 2 q 


.. (III. 56 ) 


The time dependent rf field inhomogeneity Hamiltonian, Vg(t) can 

5 

be expanded using the Magnus expansion approximation as given by 


v 2 (t) = v 2 


( 0 ) + v ( 1 ) + v 2 (2 ) + ... 


.. ( 111 . 57 ) 


To derive zeroth-order composite pulses,^ "Vg^*"*^ (Eqn. (1.25)) has 
been derived by time- aver aging Eqns. (III. 5 M 1 (ill. 55) and (ill. 56 ) 


at respective time intervals to get 


- U'lv, * b'lv, + C'l ] 


.. ( 111 . 58 ) 


Now, V, is time independent, and the time independent coeffi- 


cients are given below: 


= 0.j + B^cos^ + 9^cos(j>2 cos(<^ 2 - ^ ) 
+ 0 3 cos0 2 sin^ 2 sin(^ 2 -^ij) 


.. (III. 59) 



81 


b ' = ~ ^2 COS sin ^2 + ^3 cos ^ 1 c os 0 2 cos0 2 sin(^ 2 -^) 

- 9 3 ain^a±nB 2 sln(<f) 2 -<£ 3 ) - ^cos^ sin^cos ) .. (III.60) 

and 

c ’ = -^ 2 sin ^1 sin ^2 “®3 sin ^ 1 sin ^ 2 cos(^ 2 -^) 

+ 0 ^si n((^ 2 -(j>^) [ cos fijSin 0 2 + sinS^os^cos^] .. (III.61) 

Now to satisfy the condition V 2 ^^ = 0, we have varied the 
unknowns and found that ( 1 8 o)q— ( 1 80 ) ^ 2 q~ ( 1 8o)q is a composite u 
pulse. Similarly, variables have been chosen to satisfy the 

condition 

[V 2 (0 >, P(o) ]= 0 .. (XII. 62) 

This search produced two more sequences, namely, (90) r .-(l80) Qr -(90)_ 
and ( 90 )q~ ( 360 ) ^ 2 q- ( 90 )q . Both these sequences are composite inver- 
sion rt pulses. Zeroth— order Magnus expansion terms are given in 
Table III. 2 for all the composite pulses constructed for compen- 
sating rf field inhomogeneity. 

III.B.2(ii) Evaluation of the Performance of Composite Pulses 
as a Function of RF Field Inhomogeneitv 

Mathematical expressions to analyze the performance of 
composite pulses, viz., (Eqn. (II.27) and (Eqn. (ll.28)), 
have been derived. These analytical expressions are very similar 
to Eqn. (lll.4l) and Eqn. (111.42). However, some minor modifica- 
tions are needed to get V* and V JJ for the present case. They are: 
the effective flip angle a n in Eqns. (lll.4l) and (ill. 42) has to 
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Table III. 2 Zeroth-Order Magnus Expansion Term for Different NQR 
Composite Pulse Sequences that Compensate rf Field 
Inhomogeneity in the Spin I = 1 Case 


Composite Tt pulse 


Zero th- order Magnus expansion term , 


i) (90) 0 -(i80) 90 -(90) q Aw 1 I X3 

ii) (180) 0 -(180) 120 -(180) 0 0 

iii) (90) 0 -(360) 120 -(90) 0 1.1547 Aw I^ 3 


be replaced by 
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a n = 2 * n ( “| + A «|) » “ ( I:II * 6 3) 

£ = 90° and the efg inhomogeneity must be taken as zero, i.e., 

AWq = 0* Using these expression, we have computer simulated the 
performance of composite ji pulses given in Table XXI. 2 as a 
function of rf field inhomogeneity . ^ They are shown in Figure 
III. 4 along with single tt pulse response for comparison. 


I11.B.3 


Discussion 


Computer simulations presented in Figures III. 2 - III. 4 
illustrate that the composite pulse sequences designed through 


the Magnus expansion approach to compensate for inhomogeneities 
of efg and rf field perform significantly better than the corres- 
ponding single rf pulse. From the analytical expressions of 

V and given in Tables III.1 and III. 2, respectively, 

/ o \ 

one would expect that a pulse sequence with V n (n = 1 or 2) 

of smaller size should perform better than a pulse sequence with 

V (^) of larger size.^ However, in our case, we note the follow- 
ing interesting points. 


(i) The performance of the composite tt/2 pulse ( 90 ) q— ( 300 ) 9 q 

with larger V 1 ^ is better than that of (385) Q -(320) l 80 -(25 ) 0 
with smaller (Figure III.2? Table III.l). 

(ii) The composite inversion it pulse sequences (90) Q -(225) i8q - 
(3,5)0 with vj°Vo snows better uniform population inversion than 
(90) q -(270) 90 -(90) 0 sequence with 
Table III.l). 


= 0 (Figure III. 3 ; 
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Figure III . 4 The extent of the inversion of spin popula- 
tions, as a function of Au.j/to.j for a TT pulse 
(circles), the zeroth-order composite pulses 
(90) 0 -(180) 90 -(90 ) 0 (triangles) and 
(180) 0 -(180) 120 -(180) 0 or (90) Q -(36o) 120 “(90) Q 
(squares). 
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^ *^ r A sequence (90) Q -(36o) 20 -(90) Q with 

larger ^ » the degree of compensation for the rf field 

inhomogeneity is larger than that of ( 90 ) -(l80) qo -(90) whose 
V 2 ' is of smaller size (Figure III. 4; Table III. Z) . It may 
be pointed out here that for both these sequences [v 2 (o) ,p(o)] = 0 . 

These results clearly reveal the importance of higher-order 

terms in the Magnus expansion which are not taken into account. 

It may be pointed out that the Magnus approximation^ works well 

for small perturbations, that is, for f » V n .^* 10 Further, 

1 2 Ag)q 

Magnus expansion does not converge for larger values of ^ or 

A( \ 6)1 

L and the limitations of perturbation approach are inherent in 
CO, 

the Magnus approximation. On the other hand, it is also possi- 
ble that the presence of a non-zero V (n = 1,2) term with the 

requirement fv P ( 0 ) ] = 0 reduces the contribution from the 

higher ox'der terms in the Magnus expansion. For example, if V 
is much larger than then only those components of 

that commute with V ^ will significantly affect the performance 
of a pulse sequence. This effect is similar to that of "second 
averaging"^ discussed in multiple pulse line narrowing experi- 
ments in solid state NMR. Apart from the above mentioned results, 
we highlight the following interesting points: 

(a) The results for (90)q-(36o) ^ 2Q -(90) 0 and ( 1 80 ) Q - ( 1 80) ^ 2Q - 
(180) 0 sequences are exactly the same for a larger range of Ao^ 

as depicted in Figure III. 4. 

(b) The pulse sequence ( 90) Q - ( &) 90 ~ ( 90 ) Q works as a broad- 
band inversion 7T pulse against A(<^ and Aco, • It gives a better 
performance , for compensating A(0 q when c is anywhere between 
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150 and 300 and, for overcoming rf field inhomogeneity, may 
take values from 90 to 270 °. It gives a uniform inversion over 
(i) positive when tt = 90° and ( ii ) negative when CC = 270° 

(c) We could not construct composite 2 pulse for compen- 
sating rf field inhomogeneity. But all those composite 71 /2 
pulses which overcome AWq effect yield responses similar to that 
of a single Tl/2 pulse in the case of rf field inhomogeneity. 

(d) The results presented here are independent of the efg 
asymmetry parameter ( rj) , and hence these composite pulses should 
have a wide applicability. 

(e) It is noteworthy that most of the sequences reported 

here for spin 1 = 1 NQR case are identical to those employed in 

7-9 

NMR spectroscopy. But it is not apparent a priori that all 

the NMR composite pulses behave in a similar fashion in the NQR 
case. For example, our (90) Q - (300)^ is not a composite rr/2 

pulse in NMR, and the (90) 0 -(90)^ Q sequence, which is a composite 

, 6 
n/2 pulse used in NMR, is not so in NQR spectroscopy. 

(f) One might think that the composite Tt /2 pulse 
(385) 0 -(320) l80 -(25) 0 is equivalent to ( 25 ) Q -( 320 ) 1 8q - ( 25 ) Q 
(the flip angle 385 may be written as (2Tt + 25 ) which is equiva- 
lent to 25) and that they would perform identically. To our 
surprise, however, the sequence ( 25 )q~ ( 320 ) ^ qq“ ( 25 ) q is neither 
a composite tt/2 nor even a 3 Tt/2 pulse sequence for compensating 
Aceq. This could be because in NQR one cannot treat the effect 
of an rf pulse to be a simple rotation. 
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A general difficulty faced by us during the course of design- 
ing the abovementioned composite pulses is the treatment of all 
the variables ( 0 n and 0 n ) to be entirely arbitrary. In the next 
section, we consider only the construction of phase alternating 
composite pulse sequences. 

III.C Phase Alternating Comnosite Pulse Sequences (PACPS) 

PACPS consists of rf pulses with phase values 0° and 1 80° 
only. The total number of variables defining an N pulse sequence 
is reduced from (2N-l) to N. We consider PACPS of the format 

of TC when compared to other rf pulses. The benefits of PACPS are: 

(i) They can be easily implemented in an ordinary FT NQR instru- 
ment. However, in the case of composite pulses with arbi- 
trarily phased rf pulses, it is often necessary to use 
sophisticated digital phase shifters. 

(ii) With the reduced number of variables, the complexity in 

evaluating the Magnus expansion terms and the density matrix 
at different times would be simplified to a greater extent. 
Construction of new sequences therefore needs less computer 

t ime . 

III.C. 1 Derivation of PACPS thr o ugh Magnus. Expansion 

It is a general property of all phase alternating sequences 10 
that all even— order terms including the zeroth— order term, , 

of the Magnus expansion can be represented as a linear combination 


6L 


e 3 ^ ®2 S 3 , where the overbar denotes a phase shirt 



88 


’ w * lereas odd-order terms are proportional to 

*^Xl * ^ ie zero th-order terra, V.j ^ for a three-pulse PACPS can be 
derived based on Eqns. (ill. 29 ), (ill. 30 ) and (lll.3l). It is a 
general property of all phase alternating pulse sequences that 
the even— order terms including in the Magnus expansion will 

be a linear combination of I X2 and I^vj whereas the odd-order terms 
will be proportional to I Xl . The coefficients of I X1 , I X2 and 
I X3 °P erators in v/ 0 ) for 9 1 0 2 are given below: 

a a 0 . . (ill. 64) 

b = 1-cos8 1 +2 coa^-Sg) -cos( 0 .,- G 2 + 0 ^) . . (III.65) 

c = 2 sin^ 1 -2 sin (0-,-& 2 ) + sin^- © 2 + d^) .. (ill. 66 ) 


A zeroth-order composite pulse can be obtained if we can find 
, 9 0 and 8^ with vanishing b and c. In fact, Eqns. (ill. 65) 
and (ill. 66) are the real and imaginary parts of a single complex 
equation and the following relationship holds. 


tan [■ 


0, + a, - e n 


]. 


.. (hi. 67 ) 


When 6, = e 3 the sequence is called a symmetric phase alter- 

nating composite pulse sequence. For this case, the zeroth- 

order Magnus expansion term is similar to that of ^ 

reported by Shaka and Pines for resonance offset effects in 

the case of spin l/2 NMR spectroscopy. Therefore, the analysis 

1 3 

presented by Sbaka and Pines bolds good for our case also* 
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There is a continuum of solutions as a function of the net 
flip angTe on resonance, 0= 0^— 2 9^ • Given 0, the solution is ^ ^ 

B = arg [ 1-e -19 + (l +i4 e“ i9 + e" 219 )^ 2 ] .. (ill. 68) 

9 Z = 20 1 + 9 . . (hi. 69 ) 

where arg [re ] = D for real positive r. Hence, from Eqns . 

(III.65) and (ill. 66), we can construct composite pulses of any 
net flip angle 0 so that a uniform excitation is obtained as a 
function of Atoq to the zeroth-order in the Magnus expansion. 
Though the sequences reported in reference J are generally valid 
for the pure spin 1 NQR spectroscopy also, we prefer^ among the 
sequences listed in Table 111.3^ the following, namely, 

60 300 60 as a composite tt pulse and 114.3 318.6 114.6 as 
composite tt/ 2 pulse. Because of shorter T^ in NQR, the longer 
time duration pulse sequences given in reference 13 will not be 
of interest to us here. For the sake of completeness, we repro- 
duce the results for spin 1 NQR from literature in Table III . 4 . 
Apart from these results, we have also derived some more compo— 
site tt/ 2 and 7T pulses with smaller values and these are 

summarized in Table HI. 3- have also considered other possi- 
ble PACPS of the form 0^ For such sequences is 

again given by Eqns. (ill. 64), (ill. 65 ) and (til. 66) but with the 
reversed sign for 'b' . 

We have explored the possibility of designing two pulse 
composite pulses which would tackle the problem of shorter 
Interestingly, if ■*" s a com P os ^ te 71 pulse then 
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Table III. 3 Zeroth-Order Magnus Expansion Terms, V^ 0 ^, for NQR 
Composite Pulses that Compensate efg Inhomogeneity 
in the Case of Spin I = 1 


Composite pulse 



( a ) Composite inversion tt pulses : 


1) 

260 

80 


[-0.0585 I X2 - 0.332 I X3 ] A to, 

2) 

22.5 

225 

22.5 

-0.359 I X2 Aw q 

3) 

45 

270 

45 

[ ~°* 132 X x 2 At °Q 

4) 

60 

300 

60 

0 

5) 

90 

180 

270 

[-0.212 I X2 + 0.424 I X3 ] Awq 

6) 

260 

160 

80 

[-0.2 29 I X2 - 0.451 A(0 q 

7) 

315 

225 

90 

[-0.053 I X2 - 0.310 I X3 ] A(«)q 

(b) Composite 

Tt/2 

pulses 


1) 

22.5 

112 . 

5 

[0.359 I X2 + 0.749' I X3 1 A(0 q 

2) 

45 

135 


[0.132 I X2 + 0.768 i X3 ]Ao) q 

3) 

60 

150 


0.745 I X3 Ac^ 

4) 

114.3 

318 

.6 114.3 

0 


Table XXI. 4 Phase Alternated NQR Composite Pulses for Broadband Excitation 
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terms of the dimensionless offset parameter Atoq/2 (Uj ; •t'Total rotation in 
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^ 2 / 2 ) (that is» a pulse lasting for exactly half of the 
total time duration of the earlier one) is a composite 71 / 2 pulse 
for uniform spin excitation over larger Ao) a range. For instance, 
22.5 225 22.5 is a composite 7T pulse while 22.5 112.5 is a 

composite Tf / 2 pulse. Similarly, the composite 71 /2 pulses 45 135 
and 60 1 50 have been derived from 45 270 45 and 60 300 60 compo- 
site 7T pulses. A two pulse composite n pulse also has been 
derived as 260 80' 

Computer simulations showing the function of various phase 
alternating composite pulses have been presented in Figures III. 5- 
III. 9. For comparison, the single pulse response has also been 
given along with those for the composite pulses. 

III.C.2 Results and Discussion 

It is apparent from the computer simulations that the 
compensation of efg inhomogeneity by composite pulses is much 
better than that by a single pulse. Differences in the perform- 
ance of various zeroth— order composite pulses may be due to the 
contributions from higher order Magnus terms. We highlight some 
of the important observations that emerged from the construction 
of PACPS. 

(i) Construction of only PACPS and symmetric PACPS simpli- 
fies the Magnus expansion approach. Design of symmetric tz PACPz> 
leads to the derivation of two-pulse composite 7t/2 pulses. 

(ii) The sequence of the form ( 60 -x) ( 300 - 2 x) ( 60 -x) for 

x varying appro xxm erfce ly from 0 "to 15 is a. composite 71 pulse • 
Similarly, (l15-x) ( 320 - 2 x) (115-*) with x 0 to 10 is a 
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Figure III. 6 


2 0-5 
c 



L 

- 1-0 


Figure III. 5 Spin population inversion as a function of AgOq/ox, 
due to a single xt pulse squares ) , symmetric 
zeroth-order 60 300 60 (circles) and 22.5 225 22.5 
( triangles ) composite pulses . 

Figure III. 6 Spin population inversion as a function of AoOq/w., 
due to a single n pulse (squares), zeroth-order 
X3 270 45 (triangles) and 90 180 270 (circles) 
composite pulses. 




Figure III. 7 Spin population inversion as a function of the 
ratio of efg inhomogeneity or resonance offset 
(AcOq) to the rf field strength ) due to a 
single 71 pulse (squares), zeroth-order 260 80 
(dark triangles), 260 l 60 80 (circles) and 
315 225 90 (open triangles) composite pulses. 




Figure III. 9 
1-00r 


Acoq/oj, 



Figure III. 8 


The signal amplitude as a function of Ag^/cOj due to 
a single tt/2 pulse (squares), zeroth-order 1 ~hS 135 
(triangles) and 114.3 318.6 114.3 (circles) compo- 
site pulses. 


Figure III . 9 The signal amplitude as a function of AcOq/co^ due to 
a single rr/2 pulse (squares) , zeroth-order 
22.5 112.5 (circles) and 60 150 (triangles) 
composite pulses. 
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composite Tt /2 pulse for compensating efg inhomogeneity. This 
performance, which is insensitive to flip angle errors, Shows 
that PACPS can also be used for overcoming rf field inhomo- 
geneity problems. 

(iii) The performance of the pulse sequence 31 5 225 90 is 
identical to that of its time-reversed sequence, namely, 

90 225 315» This, however, is not the case with other non— 
symmetric PACPS. 

(iv) lhe response of all NQR composite pulses which are 

proposed for broadband excitation over is symmetric with 

respect to Aco y as in NMR case. 11 * This could be due to the 

7 

bilinear nature of the in NQR* This can be proved as given 
below. Let, U( + ) be the propagator at some positive offset 
Aco^. U( + ) will be given as 

U( + ) = exp(ij3l X2 ) exp(-ia n ) exp ( — i0 1 ^ 2 ) , . . ( III . TO ) 

and (J (-) , the propagator for the corresponding negative offset 
i.e., -Acc^» is then given by 

U(-) = exp(-i6l^ 2 ) expCiCX^I ) ex P( i ^ I x2^ *• ( II]: -7l) 

The density matrix at the end of the pulse written as 

P(t) = I^fcos^ + sin j3 cosaj + sinjS sina n I X2 

+ sin2j3 (cosc^-l) ^ Xl 


. . (III. 72 ) 



From this, one can easily show that 
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< I z ( + )> = < C- I Z (“)^> for 7 i pulses 


and <I x (+)> 


^^ x (-)^for ji/ 2 pulses. 


. . (III. 73) 


This means that the population inversion by composite 71 pulses 
and spin excitation by composite 71 / 2 pulses are independent of 
the sign of Agj^« This proof can be extended to other pulse 

sequences in NQR. 


(v) In NQR of spin I = 1 quadrupolar nuclei, the in-phase 
(with respect to an applied rf pulse along the x-axis of QPAS) 

NQR signal is given by ^ Appendix D) , whereas the 90° out-of— 
phase NQR signal is given by ^I^'y * Therefore, for an on-reso- 
nance rf pulse, <^1^^ » I^ ^ » and the 90° out-of-phase signal 
is usually detected in an NQR experiment. In the present work, 
we have given the intensity of the 90 ° out-of-phase signal to 
examine composite 7l/2 pulses (Figures III. 2, III. 8 and III. 9 ). 

The in-phase NQR signal has also been evaluated for all the NQR 
composite Tt/2 pulses. Now, the phase of the NQR signal can be 
defined by 



For a single rf pulse, the phase is zero when = 0 while it is 

a linear function of Atoq. But f° r composite Tt/2 pulses the phase 
of the NQR signal is independent of the offset (A<«^) to 311 extent 
which is smaller than that of the response of composite rt/2 



98 


pulses in NMR spectroscopy. ^ Therefore, it may be said that the 
responses of NQR composite pulses designed via the Magnus expan- 
sion approach are free from phase distortion to an appreciable 
extent . 

(vi) There exists a one-to-one correspondence between the 
spin l/2 NMR and spin 1 pure NQR spectroscopies* This is high- 
lighted by our results on symmetric PACTS. Because of this 
similarity, all those PACTS designed for broadband excitation in 

1 3 

spin 1 / 2 NMR high field case appear to be valid for spin 1 NQR 
also. This similarity could be due to the * two— level ' behaviour 
of the spin 1 NQR. This can be explained as follows. If one 
applies an rf radiation along the x-axis of the QPAS then there 
will be a coherence between only two energy level i.e., 

|0)>« ► -T [ | + f> + |-f> ](Figure III.l) without affecting the 

third eigen state. Similarly, an rf pulse considered along the 
y-axis of the QPAS induce only one transition, namely, 

|0> ^ [ |+1^ - |-1^] whereas a z-axis rf phase connects 

-Jr? [ | + l)> + |-1/> ]and [ | + 1^ - |-i)> ] energy levels. Thus an 

applied rf pulse always connects only two levels and hence the 
real three-level spin 1 problem reduces to a 'two-level 1 one. 

In fact, in 1977* Vega has used this property of spin 1 NqR 
to determine efg axis system of a single crystal specimen. 

(vii) It is noteworthy that composite tt/2 pulses for a 

three-level case in NMR (e.g. , spin-1 case) can be easily 
derived from our PACPS. If 9^ 9^ is a composite rr pulse in 

NQR then 6^/2 ^ is a composite tt/ 2 pulse for compen- 

sating quadrupole and dipole-dipole interactions in spin 1 NMR. 



It should be mentioned here that in a similar manner composite 
TT/2 pulse for spin 1 NMR have been obtained from spin l/ 2 NMR 1 7 
and, in general, they can be derived from phase alternating 
composite pulses of a two-level problem. Symmetric rr/2 PACPS 
thus constructed have been found to be giving broadband excita- 
tion and without any phase distortion in quadrupole echoes. Ve 
present the procedure of designing composite TT /2 pulses for 
spin 1 NMR using the Magnus expansion approach with a proof for 
the above mentioned concept in Appendix E. 

In the next chapter, we present the derivation of NQR 
composite pulses using the Magnus expansion approach for spin 
I = 3/2 nuclei. These results are then further generalized for 
an arbitrary spin I nucleus. 

Summ ary 

In this chapter, Magnus expansion approach has been 
employed to construct composite pulses for spin 1=1 NQR 
spectroscopy. The use of fictitious spin l/2 operator formalism 
has reduced the complexity involved in the spin dynamics calcu- 
lations to a significant extent. Composite ft and 71 /2 pulses are 

-for 

found to compensate^efg and rf field inhomogeneities to a larger 
extent compared to a single rf pulse. It has been illustrated 
that consideration of PACPS simplifies the Magnus expansion 
approach, and better results suitable for experiments in NQR 
have been achieved. PACPS offers an easy way of synthesizing 
composite 7T /2 pulses for compensating quadrupole and dipole- 
dipole interactions in spin 1=1 NMR case. It is inferred 
that there exists an one-to-one correspondence between spin l/2 
NMR and spin 1 NQR spectroscopies. 
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CHAPTER IV 


COMPOSITE PULSES FOR NQR OF SPIN 1 = 3/2 QUADRU- 
POLAR NUCLEI USING THE MAGNUS EXPANSION TECHNIQUE 
AND GENERALIZATION OF THE APPLICABILITY OF THE 
RESULTS TO SPIN 1^.1 NUCLEI 


In the previous chapter, we presented the derivation of com- 
posite pulses using the Magnus expansion method for pure NQR spec- 
troscopy of spin I = 1 nuclei. 1 Although these sequences will be 

2 i ^ 

of use in D and N NQR studies, the number of spin 1=1 quadru- 
poler nuclei that are of significance in nature are very few, 
namely ^D, 8 Li and ^N,^’ 3 Most of the systems studied in NQR^ * 3 
consist of nuclei with half integral spins such as I = 3/2 ( 33 C1, 
37 C1, 79 Br, 8l Br), I = 5/2 ( 27 A1, 55 Mn, l85 Re, l87 Re), I = 7/2 
( 51 V, 59 Co) and I = 9/2 ( 73 Ge, 209 Bi) . It would, therefore, be 

valuable to design composite pulses for spin I > 1 systems in NQR. 
One mig.ii t ask whether a strategy different from that for I = 1 ts 
needed for the design of composite pulses for spins I > 1 . To answer 
this question, the following points may be noted: 


i) Spins with I>1 may have integral or half-integral values. 

ii) Unlike spin I = 1 , 7} / 0 case, all half-integral spins will 

2 _ 

have doubly degenerate energy levels. The presence of r) 
does not split the degenerate energy levels of half-integral 
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spins in the absence of a dc magnetic field (for integral 
spins with I>1 all levels are doubly degenerate in t&e 
absence of a magnetic field except the m 1 = 0 level). A 
schematic representation of the energy levels corresponding 

to a quadrupolar nucleus with spin I > 1 is given in Figure 
IV. 1 . 

iii) By changing the rf coil axis direction in the laboratory 

frame it is possible to connect only a given pair of non- 

degenerate energy levels in the spin X = 1 case for a single 

4 

crystal. Such a clean selection is not possible for half 
integral quadrupolar spin systems. 

iv) For half integral spin systems, the flip angle of an applied 
rf pulse depends on the 77 parameter whereas such is not the 
case in a spin 1=1 system. Thus the behaviour ©f half- 
integral quadrupolar spin systems would be different from 
that of spin 1=1 case. 

Keeping these facts in mind, we decided to investigate the 

problem of design of composite pulses for half integral spin 

systems. The spin 1 = 3/2 case is an important one on account of 

its having a single quadrupolar frequency in the absence of a dc 

magnetic field. In order to study the spin dynamics of the spin 

I = 3/2 system we make use of the density matrix formalism and 

expand the density operator in terms of a suitably chosen set of 

base operators. The operators that we employ here are the irre- 

5 

duclble spherical tensor operators. Appendix B of this thesis 
summarizes the important definitions used in the tensor operator 
formalism. 6 ” 8 Tensor operator formalism has the following 



1 04 



Figure IV. 1 Energy level splitting (schematic) in 

the absence of Zeeman field for a half- 
integral spin I quadrupolar nucleus. 
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advantages over the fictitious spin— l/2 operator formalism £ ^ 

6—8 

i) Tensor operators form a complete basis set, whereas the 
fictitious spin- l/2 operators^ do not. 

ii) For spin I = 3/2, there are 18 fictitious spin-l/2 opera- 

9 8 

tor® compared with the 1 5 irreducible tensors. 

iii) Irreducible tensors obey well-known transformation laws, 

K 

following rotations of the coordinate system ♦ 

iv) Irreducible tensors carry information about the rank of the 
multipolar state and the order of multiple-quantum transi- 
tions * 

Earlier tensor operators have been employed in electron 

paramagnetic resonance 1 *"* and for the study of multiple quantum 

jC ft 1 H 

NMR, * to obtain analytical expressions for the eigenvalues 

of the NQR Hamiltonian for spins 1=1, 3/2 and 5/2 in terms of 
quadrupole coupling constant and T) . More recently, the use of 
tensor operators has been demonstrated for multiple quantum NQEt 

13 

studies also. 

In Section IV. A, we present the tensor operator represent- 
ation of the various interaction energies that constitute the 
Hamiltonian of the system. The Magnus expansion method of design- 
ing composite pulses to compensate for efg and rf field inhomoge- 
neities for the spin I = 3/2 case is outlined in Section IV. B. 
Comments on the applicability of the results to an arbitrary spin 
I quadrupolar nuclei are then made in Section IV. C. 
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XV. A Hamiltonians of the Spin System 

We consider an ensemble of non— interacting spin X = 3/2 

nuclei in a single crystal in a zero dc magnetic field. We assume 

o 

that all the spins are physically equivalent, whereby they can 
ail be defined by a common quadrupole principal axis system (QPAS). 
We further assume that the QPAS coincides with the laboratory 
frame. We consider the following interactions to be present in 
our model: quadrupole interaction {Kq) , efg inhomogeneity (V^ ) , 
the interaction of the applied rf pulse with the spin system (^ rf .) 
and the inhomogeneity in the applied rf field (V 2 ) . In what 
follows, we present these Hamiltonians in terms of irreducible 

tensor operators in different frames, namely, (i) QPAS, (ii) a 

"| 

frame where is diagonal and (iii) QIF. 


IV. A. 1 Quadrupole Hamiltonian 

In QPAS, the quadrupole Hamiltonian, a spin I = 3/ 2 

2 

nucleus is given as 


H n = 

Q 


12 


[31" -i 2 


)] 


.. (xv. 1) 


Using the irreducible spherical tensor operator (t£ , where n is 
the rank and q is the order) formalism, we can write 



(IV. 2) 


e^qQ 

zVT 


p" 


e 2 qQ 

gJz 


V 


. . (IV. 3a) 
.. (IV. 3b) 
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Ihe definition of tensor operators, commutation relationships 
and matrix representation of tensor operators for spin I - 3/2 
;ase have been given in Appendix B. 

can be alao written in the following matrix form. 



The basis vectors used to represent the 7( n matrix follow the 
order f +3/2^> , [ + l/2^> , j - 1 / 2^> and j-3/ 2^> from left to right in 

the rows and top to bottom in the columns. In QPAS, as given 
above is not a diagonal matrix and it will be useful in spin 
dynamics calculations to have the diagonalized form of Tfq. 7 {q can 
be diagonalized by an unitary transformation by the operator T» 
defined by 



where 


(i + i/x) l/2 » 


% 2 m (1 - l/X )^ 2 



(iv. 5) 


(IV. 6a) 


. (IV. 6b) 



and X = (l + 7) 2 /3)^ 2 


1 08 

. . (IV. 6c) 


The diagonalieed quadrupole Hamiltonian, ^ can then be written as 


c K 


d 


T -1 ^ q T 


On 

_s 


1 

o 

0 

0 


0 

-1 

0 

0 


0 

0 

-1 

0 


0 

0 

0 

1 


W « 
H 


* 2 qQ , 2 , sl/2 

— 2 (1 + 03 ) 


. . (IV. 7) 


.. (IV. 8) 


where d refers to the new representation defined by T and W is 

H 

the quadrupole frequency corresponding to j+J3/2^ 4 — ► |±l/2^> 
transition (see Figure IV . 1 ) . Using the matrix forms of different 
T** given in Appendix B, we obtain 

''A 


Wfd , Jia. x 2 


. (IV. 9) 


IV. A. 2 EFG Inhomo genei tv Hamiltonian 

In QPAS, the efg inhomogeneity Hamiltonian^ V^ , is given by 


V, aAP 1 T 2 + Ap" T 2 (s) 
1 o ,2 V 7 


.. (IV. 10a) 


with Ap 1 


» 2 (6q)Q 


2jS 

> 2 (<5q)Q 

67? 


.. (IV. 10b) 


and AP" 


V 


. . (IV. 10c) 
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6q is the inhomogeneity in efg. In a frame where % is diagonal , 
V 1 becomes 


V 


d 

1 



** (iv.1l) 


A(i^ is tiie resonance offset term ? or a term representing the 
distribution of quadrupole frequencies given in terns of 5 q as t 



«» 2 (6q)Q 

2 


(1 + r) 2 / 3 ) 1/2 


.. (iv. 12) 


In the observation frame of NQR ±.e., QIF l2+ defined by the 
unitary operator. 


Uq *= exp(-i t ) , 


.. (IV. 13) 


the efg inhomogeneity Hamiltonian is 


Ai 



IT “1 TT s V d 

U Q 1 u <i 1 


.. (IV. i 4 ) 


In writing the above expression, we have used the commutation 
relationship [u^,V^] = 0. 

IV. A. > RF Hamiltonian 

A general form of the rf Hamiltonian, W rf , has been given in 
Chapter II (vide Eqn. (H.io)). Using tensor operators, in 

QPAS is written as 


% rf m 2C0 t cos(w Q t + 0(t)) t] ( a ) 


.. (IV. 15) 
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where Wq is the frequency of transition between j+3/2^ and }+.l/2^ 
energy levels. The antisymmetric tensor operator T^(a) is expli- 
citly given in Appendix B. in QPAS is given by the following 


matrix! 


0 JJ 0 0 


& 0 


2 0 


o 2 o yj 

0 0 Jj G 


.. (IV. 16 ) 


with D * 0) cos (w^t + </>(*)) • 


.. (IV. 17 ) 


In a frame defined by T, where is diagonal, 'H rf is given by 


T' 1 #, 


0 K 1 O 

K 1 0 K 3 
0 k 3 g 

K 2 0 


.. (IV. 18) 


where * -1/^C/3 + W /J 3]* 

k - 2 = 1 AC /3 (x 2 -i ) 1//2 + x -i] 


and K 3 » 1/x ( 1-77 +X ) 


.. (IV. 19a) 


.. (IV. 19b) 


.. (IV. 19c) 


can be written in terms of tensor operators as 

< f -f m K,*K 3 ] Tj(a) + ^X 2 T3(a) 


. . (IV. 20) 



1 1 1 


To elltninato tfc. ti»o dependency .. make a traaafor „ atlon 

to the QIF defined by the unitary operator, U Q (see Eqn. (IV.13)) 
Now, in QIF , the rf Hamiltonian becomes 




. ( XV . 2 1 ) 


Using the transformation properties of various tensor operators 

under given in Table IV. 1 , ve can write the rf Hamiltonian 
QIF as 


JlH 


a ; 

a 


rf 


T ' (a) Bc °»(v) + 5 K 3 D + K ,('/3- 1 ) ] 


♦ 1§U) [f B SI- [?§ D *, T ?( s ) 

+ T ^ a ^3/5 D£ 1 cos ( w Q t )~ D k 3 ] .. (IV. 22) 


After truncating: the high frequency terms such as s±n(noo n t ) and 

Q 

cos(nco^t), with n ^ 1 , integer, we get the time independent static 
rf Hamiltonian in QIF as 

1 i rt » cos^t) T^(a) + ~=,cos^(t) T^(a) 

+ sin^(t) T^(s) J .. (IV. 23) 

IV. A. 4 RF Field Inhomogeneitv Hamiltonian 

Since the transformation property of the rf field inhorao- 
geneity Hamiltonian, Vg, is similar to that of the rf Hamiltonian, 
&rf* we directly give the final form of V^ in QIF as 
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Table IV. l Transformation Troperties of Tensor Operators 
T q( s , a) under the Quadrupole Interaction, 


WQ 2 

^ Q = /g’ T 0 * >or s P^- n I = 2 after Bowden et al .' 


V t J(») D <i “ j *](»)[> +3 C03(u q t)] ♦ jirf(s) «in(“ Q t) 

* Tff t 000 ^ 4 ) -i] 

U-' T*(a) * lf±Tf( s ) ,m( V ) ♦ Tf(a) cos^t) 

A. 4 3 t 3 / \ _ . _ / 


iT^(s) sin(toqt) 


U Q 1 T 2< a > U Q = T 2^ a > C0 «(« Q t) + ^ T |( s ) s in(o> Q t) 

V ,r ? (a > u « - 

♦ ir ,1 s ) »in(u Q t) + -j T^(a)[3+2 cos (u^t)] 

U- x3 t a)U,,.fi^ (s ) sin(u Q t) + T^ ( a ) CO S ( (Oq t ) 


d e^qQ 2 1 /2 

where Uq = exp(-i’)^ t); C0q = — — ( 1 +7} / 3 ) ' . For the 

symmetric combinations, namely T^( s ) , replace a by s and 

12 3 3 

s by a in the above expressions. T q , T q , T^( s) and T^(a) 

d 

all commute with 'ftq . 
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A/ 

V, 


= ~ K 1 Aw^-^cos^t) t| ( a) + ^ co s ^(t) T^(a) 
+ sin(^(t) T*( s ) j, 


.. (XV. 2k) 


Xn the next section, we present the derivation of composite 
pulses for spin I s 3/2 NQR spectroscopy using a method based on 
the Magnus expansion approximation. 


IV * B Design of NQR Composite Pulses for the Case of Spin I =3/2 
The total Hamiltonian for the spin system considered here 

may be written in QPAS as 


^ + 7i rt + v 1 + v. 


.. (IV. 25 ) 


Making use of the operator T and transforming to a frame where 

'Hq becomes diagonal, we get 


H d = T" Vt T = Xq + + V* + V 


.. (IV. 26 ) 


To suppress the main interaction and to get a time independent 

1 4 

3rf field* we again go to the QXF* Now the net Hamiltonian is 

given by 


A J 


d TT ~ 1 ^ tt -f 

n = U Q ^ U Q * V 1 * 2 


.. (IV. 27) 


Henceforth we shall drop the superscript «d' and write % d , % d f . 


V d and V* as H , K rf t V, and V ^ respectively. 


/*»/ A/ 
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To achieve a composite pulse sequence during which the 

rst nj 

interactions V^ and Vg can be suppressed, we have to tr an sfer 
<w 15 

H to the toggling frame J using the propagator (U rf (t)) of the 

pulse sequence considered. We now consider a sequence of the 

type < e t>o-( 6 2 W-( e 3>(L , where 6 n is the flip angle of the nth-rf 

2 3 

pulse of duration t Q and phase The rf pulse sequence propa- 

gator, u r f( t )» is given by 


ry 

/*©xp( — i ) f 

U rf (t) = { ex P(- i ^2 rf t 2 ) exp( - ± rf V» 


. .0^ t 


<4 

1 




expC-i ^ rf t 3 ) exp (-i '^ 2 rf t 2^ exp ^“ i ^Irf*! ^* *° ^*3 ~ 


where a) 


( 3 + 77 ) CO, 
(3+?! 2 ) 1 1 


.. (IV. 28) 
.. (IV. 29a) 


Bj.d e a - a) rf t n 


.. (IV. 29 b) 


91 „ is the Hamiltonian in QIF for the nth-rf pulse, and can 

n rf 

be written as 


n rf 


>1 K 1 (^i cos 9n T ]( a > + ~^ C ° s( k T 1 U) 

+ sin^ n T^(s) j 


. . (iv. 30 ) 


In the following subsections, we consider the effects of V. J and V £ 


separately. 
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IV * B * 1 — tn- Order Composite P ulses to Compensate for EFG 
Inhomo gene i ty 

In the toggling frame, spins evolve under the transformed 
efg inhomogeneity Hamiltonian (^(t)) alone (neglecting the 
effects of Vg) (refer Appendix C) and the propagator is given by 


t ~ 


U v (t) = T exp[-i / Vl (f) df ] 


.. (IV.31) 


To get the effective time independent efg inhomogeneity Hamil- 
tonian, we use the Magnus expansion 1 ** to write 


ft} 


v t (t) = V^ + ... 


.. (IV. 32 ) 


For the spin 1 = 3/2 system, the transformed Hamiltonian, V (t), 
for a three-pulse scheme has been derived in the following way: 
first , the commutators between the rf Hamiltonian and different 

M 

tensor operators, that is , sure evaluated, and these are 

listed in the Table IV. 2. Then, using the nested commutation 

relationship (see Eqn. 1.21 ), the transformation properties of 

different tensor operators under the effect of rf pulses have 

been derived, and these are summarized in the Table IV. 3. Now 

9S 

using the results of Table IV. 3 V^(t) is obtained as, 

+ a 3 (t)T](a) 


V 1 (t) = ( t)T^ + ia 2 (t)T^(s) 

+ a^(t)T^(a)] 


.. (iv. 33 ) 
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Table IV .-2 Commutation Relationships of few Selected Tensor 
Operators vith the Hamiltonian, % 


AJ <\r 

’Hz: % + V + V 

nrf + V 1 + v 2 


<y „2 


[*, '^(s )] = 


[ T J(a)] 


[ % , T^(a)]: 


[ £ , T?(s) 


[ £ , T^(a>] « 


[ # • T?(a)] 


-£<«, + Ao> 1 ) T^( s ) cos^ n 

+ f J? i ^ C0 1 +Au 1 ) T ](a) sin ^ 

+ y== i(co 1 + Aw 1 ) T^( a) sin^ E 
-V3(w i +Aco 1 ) T^ cos^ n + ■“ Au Q T\j(a) 

+ ^ T ? (a) 

“JF i (“l + A V T f sin & " ( “l + AW i> T f( s ) sin ^ n 

f Aid 

(W 1 +Ao> 1 ) T~(a) cos^ n + *^=? T^(s) 

- (o^ +Aw 1 ) T;^ ( a ) cos^ n - i(« 1 + A(o 1 ) T^ sin^ n 

+ i(w 1 +Aw 1 ) T^®) Sin ^n + \/f ACi) Q 

+A ^) T ]t a ) sin ^n +A( °1 ) T 1 (a) S±I Vi 

*.^'i(a> 1 +A(0 1 ) T^( a) sin^ n+ -^Aco Q T^(a) 

- ( «1 + A( 0 1 ) T^(a) cos^ n + (a) 1 + Aco 1 )T^(a)cos^ i 

+ (<0 1 +Acb 1 ) t](&) co.£ + -J Ao)q T^t*) 

-/I* (c^+Aw,) l|(a) 008^+ (w^ Aca,) T|(s)sin^ n 


• • . « contfd.* 
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Table IV. 3 Evolution of Tensor Operators under the Effect of 
the Hamiltonian, % . 


Hi Z ~ ~ 

+ v + V 

nrf 12 


U T 2 U" 1 
o 


Tq[cos j 5 + sin^gcosa^] + 

.2 


+ iT^(s) sinj 3 [ cos^sina^ + sin^cos8(cos a ^-1 ) ] 
”[3 a ) + T i ( a )] sin 3 [cos 3cos^ n (co SOn -1 ) 


-sin^sin^] 


1-1 

U T (a) U = sinjS [cos j9cosd( 1-cos a )-sina sinC 1 ] — 

1 ' 'n n n rn J 

+ — i— T (s)[sin2flS cosa siB^ -sin2$ sin^/3 - cos/3 sina ] 
2 {2 n n r n n 

. 3 r* -» 

- iT ( a) [sinflf cos/3(cosa - 1 ) +cos <& sina ] sin/3 
j *n n rn n 

« n 

+ T^s) sinjS [cos^cos^t co sa Q -1 J-sin^sina^] 

+ ^ T^( a) sin^jg (l-cos^) 

+ — T^j (a)[ 3cos 2 0 n (sin 2 3 + cos 2 ^cosa n ) + 3sin 2 ^ n cos 

o 

+ sin jSO+cosctjj) +2] 


( a) [2cos^ n ( sin^/3 + cos^cos tt^) 
n— % 




vT5 


+ 2sin ^ n cos c^-sin jg ( 1 + cos^)-2] 


* * . contd. 



U T 2 (s) u" 1 


2 

-iT^ ®in3[ sin^cos 3( co s a n - 1 J-cos^sina^] 

+ [f /f iT ?^ a ) + -^iTjCa)]* 

‘ [ cos^sin^sin^ + sin^cos^costt^-cos^sinCl^cosP ) 
-sin^lcos a n c °s^ n + sin^cosgsinc^) ] 

+ T 2 (s)[sin 2 ^ n ( cos 2 ^ cosa n + sin 2 j3 ) + cos a^cos^] 


U T^U) U' 


T q sin|3[cos3 cos^Cl-cos^) - sin^sin^] 

+ T 2 ( 8 ) [ —— sin2<5^(-cos c^+cos 2 |3 cosa^+sin 2 ^ ) 

“sina^cos^] 

iT^Ca) sinj3[ sina^cos^+siny^cos^ ( l-coso^)] 

+ — — T 2 (s) sin^[ siny^sin c^+cosScos^C 1-cos 0*)] 

2 vjio 

+ T^( a) [sin 2 g ( cos a - 1 ) ] 

2 J 5 J 

3 17*1/ . 1 J _2^/ „ \ , , 2 j 


+ 5 V5 T i( a )[ _1 " 2 sin 3( cosa n -1 ) + sin ^n. 00 S °ii 


+ cos^ n (sin^3 + cos^jScos 0,,)] 


+ ^ T^(a)[l + £ sin 2 0(coscx n -l) + ~ sin^eos 
+ ~ cos 2 y$ n (sin 2 0 + cos 2 /3cosa n ) ] 


• • • contd « 
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'able 


contd. 


(»> U 


2± j 1 

s i n 2^ n ( sin 2 $ + cos 2 $ cosa^-cosCt^) 
+ cosj3sina n ] 


+ , ^2^®)[co s 0L n (si n ^ ri + cos^ n cos^ ) +cos^ n sin^3] 

r-*~* 

+ T (a) sin0[ cos^co s3cosa n + s ina^s in co 38 co s^] 

+ if T 1 (a ) sin ^ CO s ^COS0 ( 1~cosa n ) - sina^sin^] 

+ T^( a) sinS [ cos^ i cos/3( 1-cosa^) - sino^sin^] 

U ( a) U" 1 = - ( a) sin/3 [cos^sin^ + sin(^cos<^-cos/3cosa n sin(^ 

+ T^(a) sinjB [cos^sin^ n ( 1-cosa^) + sin^cosj^] 

+ -j=, T^( a) sin0[cos^sin^ n ( 1-cosa n ) + sina^cos^] 

+ T^(a) [ cos^co s 2 ^ n + sin^ n ( sin 2 ]3 + cos Ciscos 2 /3 ) ] 

- ^ iT 2 (s) [sina^cos/B sin2^(-cosa n + sin 2 8 


+ cosa cos 2 0 ) ] 


u t3(.) U- 1 


— 2 - T? ( a) sin 2 ft ( 1 -cosa) 

10 1 11 

+ — ^n,T^(a) sin 2 3 (cosa n -1 ) 

2 v15 

+ T^(a)[l sin 2 3(cosa n -1 )] 

+ T 2 (s) sin 3 [cosgcos^C 1 -cosa^) + sin^sint^] 

2/2 2 

+ -i- T?(a) sinj3 [cos0sin^(-1 + 003 ^) +sina n cos^] 
rr 2 


# • • co q t d * 
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where 

/c°«( w# f f t ) ’ 


h (t): 


COs6) 1 C ° 8 ^efft 2 )-»in (?, sin(0> eff t 2 )co^ 2 , 

cos ^[cos(w eff t 3 ) co s e 2 -sina 2 sin (w eff t 3 > cos^-^)] 
-sin (cos9 2 -1 ) s in(w eff ,t 3 ) [1 sin 0 3 sin 2$ 2 + cos^cos^] 

^-sin [® in ^ 2 C08 ^2 c ° 8 ( a) efft 3 ) + sin ^ w e ff *3^ cos ^j] 


0 - 


0 £ *2“ 


a. 


“rf 

Ja 

Ci> 

rf 


0 £t 3 £ 


“rf 


.. (IV. 34) 


/sin(<o 0ff t 1 ). 


O^rt^ 


ft 


1 


a 2 ( t ) ={ 


sin 6, cos(w #ff t 2 ) + cos P lS in(0) eff t 2 ) cos^, 

-sin© 1 [ sin0 2 cos(^ 2 -^ 3 )sin(a) eff t 3 )-cos0 2 cos(w ef . f> t 3 )] ^ 

+ cos P, (cosP 2 -1 )sin(w eff t 3 ) [cos^cos^ +j sin20 2 sin0 3 J * * 
\^+co S e i [sine 2 cos9* 2 cos(w eff t 3 ) + cos^sin( % ff t )J 


0 -*2- 


“rf 


CO. 


-*3 


rf 


O'- 


co 


rf 


/n 


J 

.. (IV. 35) 

o £t. 4 


e. 


•j(t) ={ 


5/2 

6 


/ 

* sin$zisin( co ff t ) , 


1 W 40 

rf 

^ e, 

o6t 2 *=;r~ 

rf 


[sinP 2 sin^ 2 cos(( i ) eff t 3 ) + sin^ 3 sin(co eff t 3 )] 


V 


(co»P 2 -1 )sin(cO eff t 3 ) [sin20 2 cos^ + 

( 1 -cos20 2 )sin^ 3 J 


o£t 0 £ 


3 " co. 


rf 


. . (IV. 36 ) 
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and 


( O 
-4 


a. 


a k( t ) 


/30 ® in ^ C eff t 2 ) S W 2 * 

[sin® 2 sin^ 2 cos(w eff t^) + sin^sin((O ef>f .t^)] 
y==r (cos0 2 -1 )sin(co ef>f t^)[-cos^sin2^ 2 

+ sin(^(cos 2^ 2 -1 )] 


. . 0 £ t - — 

1 W rf 

_ o£t 

rf 



.. 0 ^t„£ 


( 9 . 


~ 3- u. 


rf 


.. (IV.37) 


The zero til-order Magnus expansion terra, V^®^ (vide Eqn. X.25) can 
be obtained by integrating Eqns. (IV.34), (lV.35) , (IV.36) and 
(IV.37). Thus we get, 


V 1 0) S 7 *- a 1 T o + ia 2 T ? (s) + a 3 T J ( a ) + a 4 T^(a)].. (IV.38) 

where = sin ^ [ 1 + (cos 0 2 -1 )cos^ 2 ] + cos 9 ^ [sin © 2 + cos & 2 sin B^] 

+ (cos 9^-1 )[sin 0 2 cos 9 ^ 003 ( 02 “^) + sin 0 ^ cos^ 3 J 
-sin sinf^sin ( 9 ^cos^ + (cos 6 * 2 “ 1 ) (cos 0 3 -1 )sin 0 1 • 

.[^ sinjZ^sin2^ 2 + cos^cos^] 

.. (iv. 39 ) 

a 2 * cosS^I +cos(25 2 (cose 2 -l) - sin^sin^cos^ + cos^cos^-l )] 
-sin © 1 [sin8 2 +cos0 2 sin & 3 + sinfl^eos&j-l )cos (0 2 - ] 

+cos0 1 (cos9 2 -l)(cos0 3 -l) [j sln2$ z sln(j> 3 +cos ^cos^]-! 

.. (IV. 40 ) 
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a 3 = Tff (c ° 8 V 1 H2 sin^ 2 +cos^ 3S in2^ 2 (cos6 3 -l) 

+ ®± n ^j( 1-cos2^) (cos 0^-1 ) ] 

+ [®ine 2 sin9 38 in0 3 - s in^ 3 (co S e 3 -l)] .. (IV. 4l) 

and 

a 4 = 435 (cos0 2 ~1 )[sin^> 2 + — 3 . si^^cosG^-l) 

“ — (cos 2^-1 ) (cos^-1 ) ] 

+ ^35 [sin0 2 sin© 3 sin ^ 2 - sin^cos&j-l ) ] .. (IV.42) 

The condition that a zeroth— order composite pulse sequence 
should satisfy vj ^ e O,^ demands that all the coefficients a^ , 
a 2 » aii£ ^ ®2j. ^ lav ® to be individually zero. But the condition 

that [v, (o) , P(0)] = o ' 8 needs only a 2 , a^ and to be zero. Using 
these conditions, we have arrived at many sequences for spin 
1=3/2 case. These sequences with their zero th- order terms, 

V ^ , are given in Table IV. 4 . It should be mentioned here that 
all these sequences which compensate for efg inhomogeneity in the 
case of spin I = 3/2 have been discussed in Chapter III in the 
context of spin 1=1 nuclei (see Tables III.1 and III. 3). The 
underlying similarity between these two cases now emerges. 

IV . B . 2 Zeroth-Order Composite Pulses to Compensate rf Field 

Inhomogeneity 

The time dependent rf field inhomogeneity Hamiltonian in 
the toggling frame 10 of a three-pulse sequence has been derived 
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Table IV. 4 Zeroth-Order Magnus Expansion Terms^vj 0 ^ , for NQR 
Composite Pulses that Compensate for efg inhomo— 
geneity Effects in the Case of Spin I = 3/2 
Quadrupole Nuclei 


Composite Pulse 


( a ) Composite tt/2 pulses 
22.5 112.5 

45 135 

60 1 50 

( 90 ) 0 -( 300 ) 90 

114.3 318.6 114.3 
385 320 25 


vj 0) (>/r “ rf / Aw q > 


1.765 if + 0.848 iT^(s) 

2.4l if + 0.4l4 iT^(s) 

2.73 if - 0.149 ilf( s ) 
if -0.134 iT^(s) + 0.424 Tj ( a) 

+ O.365 T^(a) 

0 

0.033 if - 0.033 iT^(s) 


260 

80 


-1.97 if - 

0.35 iTf(s) 

22.5 

225 

225 

1.67 iT^(s 

) 

45 

270 

45 

0.83 ilf(s 

) 

60 

300 

60 

0 


90 

180 

270 

4.0 T 2 - 2 

O 

.0 ilf( s ) 

260 

i 6 o 

80 

-3.939 if 

- 2.0 ilf ( s ) 

90 

225 

315 

3 .4l4 T 2 - 

O 

0.586 iT^ ( s ) 

315 

225 

90 

-3.4l4 T 2 

O 

- 0.586 iff (s) 


(90) 0 -(270) 90 -(90 ) 0 

w rf * (3+r»)<h 1 /(3+77 2 ) 


where 
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using the results of Tables IV . 2 and IV . 3 to be 

V 2 (t) * ACD-Jb^tjT 2 +b 2 (t)T 1 1 (a) +b 3 (t)T^ (s ) +b 4 (t)T^(a>] 


. (IV. 43 ) 


r° » 

- ^ sin 6 ^ sin ^ , 
t >1 ( t) =/ 2 cos ^ sin ftjCos^f cos^-sin^) 


» • 0 <^t . - 

1 CD. 


rf 


G6t„^ 

2 CD. 


^2 


rf 


+ sin ^ (cos 0 2 -1 ) cos0^sin2<^ 2 -cos 2 ^ 2 sin^ l 


#3 

0 ^ t 3 ~ u 

J rf 

.. (IV. 44) 


j sin2<^> 2 sin^ ] 


2 

5 ’ 


e. 


0 £ — 


1 6 ). 


a. 


rf 


b 2 (t) = \5 COS i 4 ’ 

^ cos^ + Yo ( cos ^2 -1 ) [cos^( 1- c °s20 2 )-sin^sin2^2 J- 0 <lt 0 ^— 

V 


c 3" CD. 


rf 

(IV. 45) 
01 


- • o ^t ^ ““ 


b„(t) 


_i 

2 


cos8^sin^ 2 , 

/— sin^sinSgCos^cos^-sin^) 

j 2 


— “ COS0., (C0s0„- 1 ) B COS <j> S±n2± - 5 cos^sin^ 
5^2 1 2 J 


r “rf 
^ $2 

. . O ■£ - — 

2 “rf 


„ e 3 




+ sin2^ 2 cos^ 3 ] 




rf 


. . (IV. 46) 
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and 


r 




2 

>[T 5 * 


. . 0 fr t.b 




1 ~ « 


rf 


Trf °° s (4 > 

j=~ (c°s9 2 -1 )[cos^( 1-c©s 202) -sin^siBL2^ 2 ] i 


. . o £t_£ 


< 3 * 


1:75 


cos 


A 


r .. o £t 


<*>rf 

8 - 


rf 


.. (IV. 47 ) 


Using the Magnus expansion, ^ ^ we write 


V £ ( t) = V^° ^ + V^> + ... 


(IV. 48) 




To get an effective time independent Hamiltonian, Y^, in the 

1 5 /. T (n) 

toggling frame, we truncate all higher order terms (Y^ ; 

n ^ 1 ) in the series. By integrating Eqns. (IV.44) -(IV. 47), we 

get the zero th-order term, V^ 0 ^ as 

V 2 0) * f b 1 T 0 + b 2 T 1 (l> + b 3 T ? U) + V? (a) l •• (lV - 49) 

rf 


where b 


1 


& 0 

2 sinfijSin^ + ^ cos ^ sin8 2 cos^ 2 (cos^-sin<5^) 


/2 

0 , 


sine^cose^l) cos0 3 sin2^ 2 -5cos 2 ^ 2 sin^ 

i 


4* smc ? 1 \cos& 2 ~ 3 ) ^ 2 ww=> ^3 a>J "“ fc ‘ ST2 W ~ / 3 

4-sin dsin 2 ^ l .. (IV. 50) 
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3 

iss 

2 -5 


3 2 


and 


b. * 


* v/15 


*•* 

+ ^ 2 co ®^2 + ® 3 c °s 0 3 

® 3 (cos g-1 ) [ g c os^ 3 (i-cos 2^) - sin(^ 3 sin2 

.. (IV. 51 ) 

OgCos^sin^ + ^sinSjSin&j cos ^ 2 (cos^-sin^) - 
fijCos 6 j (cos 1 ) cos^ 3 sin2^2- c °s^2 s i 11 ^ 

+ ^ sin2^2 c os^ 3 ] ^ .. (IV.52) 

2(6, + 0 2 CO s^ 2 + @ 3 cos^ 3 ) ~] 

+ O 3 (co 80 2 “1 ) [l-cos2^-sin^ sin2^ 2 ] I • • (IV.53) 


Using the conditions V^ 0 ^ = O 17 and [v^°\ f(o)], 18 we 
obtain expressions for some of the zeroth-order composite n 
pulses for comensating the rf field inhomogeneity . These 
sequences are given in Table IV. 5 along with their zeroth-order 
Magnus expansion terms. To our surprise, these sequences are 

1 2 

exactly the same as that for spin 1=1 composite pulses 
presented in Chapter III. 


IV. B. 3 Efficiency of Compo site Pulses 

The degree of compensation of a composite pulse against 
©fg and rf field inhomogeneities can be obtained by following 
the time-evolution of the spin density matrix under the action 
of the applied pulse sequence (see Appendix D) . The 
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Table IV . 5 Zero th-Order Magnus Expansion Terms, Tor NQR 

Composite Pulses that Compensate Tor rT Field 
Inhomogeneity in the Case oT Spin I = 3/2 Quadrupole 
Nuclei 

Composite Tt Pulses w r T^ w 1 

(90) Q -(180) 90 -(90) 0 -1.57 + 0.94 T^Ca) + 0.81 T^(a) 

(l80) 0 -(l80) 120 -(l80) 0 . 0.85 T^(a) - 1.92iT^(s) 

(90) 0 -(360) 120 -(90) 0 -2.72 T* - 1.41 x](a) -1.22 T^(a) 

2 l/2 

where «* ( 3+7} )to .j / ( 3+7) ) ' 
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reduced tUeneal equilibrium density uatrix in QIF is siren by 


0? d 

P(o)*** 


. . (IV.54) 


Tlie state of the spin system at the end of a composite pulse 
sequence can be defined by the density matrix, f ( t ) , as 


P (t) = u(t) e (o) uur 1 


.. (IV.55) 


where U( t) = T exp [-i f K( t') dt' ] 


.. (IV. 56) 


The effective Hamiltonian in QIF during an nth— rf pulse in the 
composite puls© sequence is 


/Vi (TV 

X--H 

11 


rf 


+ v t + V. 

I * 


.. (IV.57) 


P( t) can be evaluated using the nested commutation (see Eqn. 

rJ f** 

(l.2l)) between 'K and P(t). We have evaluated the commutation 
relationships between the Hamiltonian 'H and different irreduci- 
ble tensor operators T^s.a), and the results are given in 

q 

Table IV. 2 . Substituing the results of Table IV . 2 in the nested 
commutation relationship, the transformation of various tensors 

fs/ 

under the effect of H have been derived and they are summarized 
in Table IV. 3. Now, using the results of Table IV. 3 , one can 

f ^ 

easily calculate the spin density matrix P( t ) for any type of 
pulse sequences in NQR spectroscopy. For a three-pulse composite 

r ✓ 

puls© | v© present in the following P(t) at different time 
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interval, of the sequence . a. density Batril sfter th . 

firat pulse , namely, ( ft, ) o in the composite pulse sequence, is 
derived as 


^ ( 1 1 ) ~ T 0 C cos £ + sin $ cos O,] + iT^(s) sin$ sina^ 

1 

- T -j(a) sin 20 (cosc^- 1 ) -T^(a) ~^=z sin 2 £(cosa.,- 1 ) 

.. (IV. 58) 


where /3 = tan 


VViV 

*- nr J * 

V6 Aoi-t 


and 


XV 


j» W . I T 

1 °= / 3 cosX + sinX 


.. (IV. 59 a) 


a i e [Ag^^ cos £3 + K 1 (w 1 + Aw 1 )sinjS]t 1 \ .. ( IV. 59b) 


.. (IV. 59 c) 


Oj is the effective flip angle of the rf pulse when V and are 
not negligible whereas 0, is the flip angle of the same rf pulse 
when the effects of V., and V^ are not considered. Now, after the 
application of a second rf pulse, namely, ($2 )(^ * the density 

matrix in QIF becomes 


P(t 2 ) = + iB^s) + C^TjU) + ^^(a) 


.. (IV. 60) 


where A 1 = cos^ + sin* 1 # cos 0,003 + cos% sin% (0030,+ cosc^) 

- ~ sin 20 sin 0 (cosa.j -1 )[ cos£ cos^ 2 -sina 2 sin ^ 2 

-cosj 3 cos^ 2 cos (X2 ] 

-sin Oj sin^B [ sin^cosg ( 1 -cos a,) + sim 2 cos^] 


.. (IV. 6l) 
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B, - »ing[co s 2 B . «in 2 ecosai ] [cos^sina., + ain^cosB (co.Oj-t)] 

+ i * in2 e (oo^-Df-sia^ cosa 2 «ln 2 S + al a 2 ^sia 2 P 

+2cosjg sin c^] 

+ sin dj sin 3 [sin 0 2 (sin 2 0 +cos 2 3cos«C^) + cosdgCos 2 ^] 

. . (IV. 62 ) 


y m il 

— y] Q / 2 2 * 

5 »in0(cos 0 + sin 3 cos Oj ) [cos 3cos^ 2 (cos a^-1 J-sinC^sini^] 

~ }o sin2 $(cosa,j -1 ) jcos 2 <^$ 2 ( sin 2 # + cos 2 0cos a,)+sin 2 <^ 2 cos ] 

- sin^sin^ [cos^ 2 (sin^ 2 sin 2 3 +sin^ 2 cos 2 j3cosa 2 

-cos^sin agCosg) 

-sin^(sin0 2 cos3sina 2 + cosa 2 cos 4)] •• (IV. 63 ) 


D 1 = - 3 if 


3 Jj [cos 2 g + coso^sin 2 #] s±njg[cos 0 cos^ 2 (cosa 2 -1 ) 

-sin 

/ 2 r 2 # 2 2 2 / 

in20(cosa i -1 ) /— [cos Q^lsin g + cos |3cosa 2 ) + sin ^ 2 co 5a 2 J 

■j y^sinpsina^ [cos^ 2 (sin^sin 2 3 + sin^ 2 cos 2 j3cosd 2 

-cos^ 2 sino2 C °s3) 

-sin^ 2 (sin^ 2 cos.3sina 2 + cosOjCOs^)] 

.. (IV. 64) 


=* [ Aa^ /T cos 3 + + Ato^ ) sin3 ]t. 


.. (IV. 65 ) 


The density matrix after the third pulse, P(t^) is given by 
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?(t 3 ) = A 2 T^ + iB 2 T^( s ) + C 2 T](a) + D^U) .. (IV. 66 ) 

A m A [cos 0 +• sin 2 0 cosn 1 

3 J 

+ [sin^cos^ ( 1-cosa-j) - cos^sina^] 

c 1 

+ 't~2 s ^ n )3[ cos l3 co s^j + sin^sin^ - cos|5coso^cos^j] 

/ Q 

+ D 1 Y10 sin 3[ cos 3 co s^-j + siaa^sin^ - cosjScos^cosoj] 

.. (IV. 67 ) 


B, 


A^ sin0[cos^sin 0^ - sin^cos0(cos a^-l ) ] 

+ iB^ jsin 2 ^ [sin 2 0 + cos^cosa^] + cos^cos( 0^ +^3 )^j> 

c r 

+ “J-, |sin2^sin 0 ( I-C 05 CI 3 ) -2cos0si nCLjJ. 

r~r fsin2^3 2 2 1 

+ Dj | ' 1 (sin 0-COSCI3+COS 0co S a 3 )- S ina3Cos0 V 

.. (IV. 68) 


3C 


[cos cos^cosa^ + sin^0 ) + sin 2 ^> 3 sina 3 ] 

Ci 2 

+ [sin j3 ( — 1 +cosa^) +2] 

- — ^ B 1 [~cos9^(sin^3sin 2 0 + s±n(p^cos Z ^cosa^ 

+ cos^3sina 3 cos0) 

+ sin^(cos Ojcos^-sin^cos^sina^)] 

*» • 

_ [ si 2 ' 2 -^ cos^(cosOj-l) + sina^sinjSsin^ ] 

+ | J| D 1 [-1- (cosoj-l) + sin^cosOLj 

+ cos^ 3 (sin 2 0 + cos 2 0 cosa 3 )] .. (IV. 69 ) 
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and 


D 2 “ 


“ J 1 5 " s: *- n 0 4.^ [cos j5cos^(cos <Xj- 1 ) + sinOjSin^ ] 

+ \j ~^5 B 1 I- cos ^*j(sin^ 3 sin 2 |3 + sin^cos 2 £cos Ci^ 

+ cos^sin^cosg) 

+ sin(^(sin^cosj3sin cy» cos a^cos^)! 

°1 2 

+ [ 2cos ^-j( sin £ + cos^cosa 3 ) +2 sin^^^cos Cf.^ 

+sin 2 £ ( 1-cosa^) -2] 

2 

+ 3 D 1 [l+ ^2^ (cosa^-l) + -• sin 2 ^cosa^ 

+ ^ cos 2 ^ (sin 2 0 + cos 2 0 cosa J] . . (IV.70) 


a 3“ COS0 + ( a> 1 + Aco.j )K^ sinjB] ( XV . 7 1 ) 


Now, the efficiency of a composite 7l/2 pulse can be judged by 
measuring the signal intensity, (refer to Appendix D) as a 
function of A(<^ or Aco^ . V^' is given as (see Eqn. ( 11.27) ) • 


jfe T i(») f P(S) ] 

M (fjf T 2 ( s ) ) 2 ] 


.. ( iv . 72 ) 


where iT 2 (s) is the tensor operator corresponding to the obser- 
vable spin magne tization in QIF, Similarly, the extent of spin 
inversion by a composite ji pulse is given by (see Eqn* (XX * 28 ) ) 


~Tr[T 2 t p(t 3 )] 

Tr[(Tf) 2 ] 


.. ( iv . 73 ) 


where (l£)t is the adjoint of T* expressed as 
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(-l) q T n 

-q 


.. (IV. 7k) 


Analytical expressions for and V * have been derived 
and, interestingly, these expressions are similar to those for 
V 1 111 • l H) and V (Eqn. III. 42). However, ¥ and ¥ 1 

differ only in the definitions of <X n (effective flip angle of 
an nth-rf pulse in the composite pulse sequence) andj9(i.e. 
the angle which defines the direction of the effective rf field). 
It has been proved through computer simulations of V * and 
(as functions of and ) that ¥^ s and W^ 1 behave in a manner 

similar to and W^, respectively. Therefore, all those compo- 
site 71/2 and 71 pulses reported in Tables IV. 4 and IV. 5 for spin 
I = 3/2 case, perform in a manner similar to that of spin 1=1 
composite pulses as shown in Figures III. 2 -III. 9. 

In the next section, we will discuss the applicability of 
these composite pulses to spin I quadrupolar nuclei in general. 


IV. C Comments on the Annlic ability of Comnosite Pulses for 
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effective flip angles (aj , rf pulse phases ( £ } ^ the 

( 3 ) which defines the direction of the rf radiation. Expressions 
for <X n and 0 depend upon the spin value I and the asymmetry para- 
meter (r)). Apart from the difference in the definition of a n 
and 0 , analytical expressions for ^I^ , ^ I ^ and^X^ are very 
similar irrespective of the spin value. This has been proved 
for spin I = 3/2 and 1=1 problems in the last section. By 
mathematical induction, the proof follows for any value of X. 
Therefore the expression for testing composite pulses, namely, 

W.j and (see Chapter III), are valid for a general spin X. 

Hence composite pulses reported in Chapter III will be applica- 
ble to an arbitrary spin I case. The general expressions for 
a n and 0 are given below: 

a n = t n [ A( jDq COS0+K' sin0] .. (lV.75) 

and 0 = tan” ^ [777—] .. (IV.76) 

Explicit expressions of and K’ for different spin I values 

have been summarized in Table IV. 6 . 

In this context, it is worthwhile to understand the two- 
level nature of higher spin systems in the absence of Zeeman 
field. Though higher spin systems (12-3/2) possess multilevel, 
an applied rf radiation connects only two doubly degenerate 
energy levels. For example, in (i) I = 3/2: |+l/2 )*— * |+3/) , 
(ii) I * 5 / 2 : | + l/ 2 >«-* |± 3 / 2 ) or |+3/2)«-» |± 5 / 2 ) and similarly 

the single quantum coherences in I = 7/2 and I = 9/2 case (see 
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Table XV. 6 Expressions for efg Inhomogeneity ( Aco^) and the 
Effective rf field strength (K' ) for Different 
Spin Values 


Spin of the 

Quadrupole Ato K 1 

Nuclei, I ** 


1 e 2 (<5q)Q 

3/2 ii ±rt 2 jW /z e 2 (6q)(5 

5/2 (i) ^ e 2 (6q)Q for | + l/2^> 4-4 |+3/2^> 

T)«= 0 

(ii) e 2 (6q)Q for |+3/2> 44 |+5/2> 

7/2 (i) for |+l/2> 4-4 |+3/2> 

71= 0 

(ii) 9 for |+3/2> 44 |+5/ 2> 

(iii ) e 2 (sq)Q for |+5/2) 44 J+7/ 2^> 


2 (g> + Aco 1 ) 

(3+T7 2 ) 17 
/8 C^+AU ) 


/5 (w^ + ) 

/l5 (“+ A V 
JTz 

Jl (^, + Aw) 
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Figure IV. 1 ) . This could be the reason for the similarities 
noticed in the analytical expressions of 1^ (where p = x, y, 
z) for different spin I cases. However, this may be no longer 
true if one does not considers non— interacting spins in the 
system • We may therefore propose spin 1=1 pure NQR composite 
pulses for higher spin systems in pure NQR spectroscopy* 

Xn the next chapter 9 we discuss the procedure for pulse 
response studies of NQR powder specimen stitI present a powder 
averaging method for general spin I nuclei. ¥e also present a 
numerical procedure for designing composite 7T/2 pulses to enhance 
the intensity of NQR signals from powder specimen. 

Sutnm ary 

Details of construction of composite TC/2 and 71 pulses for 
compensating efg and rf field inhomogeneities have been presented 
for spin I = 3/2 case in pure NQR. The irreducible spherical 
tensor operator formalism and the Magnus expansion approach have 
been used for constructing composite pulses* Transformation 
properties of the various tensor operators involved have been 
tabulated, and these results could be used for constructing 
higher order sequences and also for studying the time evolution 
of the spin systems under different kinds of pulse sequences* 
Comments on the applicability of composite pulses for arbitrary 
spin X quadrupolar nuclei are made on thebasis of the similarity 
in the performance of composite pulses fcr spin 1=1 and I = 3/2 


cases . 
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CHAPTER V 


DESIGN OF COMPOSITE 71 / 2 PULSES USING A 
NUMERICAL APPROACH TO ENHANCE THE NQR 
SIGNAL INTENSITY IN POWDER SPECIMENS 

Large single crystals of most inorganic and organic compounds 
are not easily available. Smaller sized single crystals do not 
fill up the rf coil used in NQR spectrometers and the poor filling 
factor leads to a poor signal intensity.^ Polycrystalline materials 
are therefore commonly used instead of single crystals in NQR. In 
this thesis, we have so far considered composite pulse sequences 
for single crystals and developed the theory accordingly. The 
validity of such sequences for powders can be demonstrated theore- 
tically only if powder averaging is carried out taking into account 
the random orientation of the efg axes of the crystallites consti- 
tuting the powder specimen. 

NQR signal intensities observed for most polycrystalline 
samples by the application of a single rf pulse are generally 
poor 1 and in some cases it may not be even possible to detect the 
signals. This is because the strength of the rf pulse felt by the 
quadrupolar nuclear spins in different crystallites varies and 
only those nuclei whose moments precess perpendicular to the rf 
coil axle produce the maximum signal. It would be therefore worth- 
while to design composite n/2 pulaes to excite all the crystallite. 
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uniformly in order to achieve better signal intensity. In Section 
V.A of this chapter, we present a theoretical model for calculat- 
ing the rf pulse responses. Section V.B introduces a numerical 
approach to the synthesis of composite tt/2 pulse sequences for use 
with powder specimens. 

V. A. 1 Model for the Spin System 

In a poly cry stall ine sample consisting of quadrupolar nuclei, 
the QPAS of each crystallite is randomly oriented. Hence, unlike 
in the case of a single crystal, a common QPAS cannot be employed. 
Instead, we consider an ensemble of isolated non-interacting 
quadrupolar nuclei whose QPAS ' s are randomly oriented with respect 
to the laboratory frame. The angles a and 3 (see Figure V . 1 ) 
define the orientation of V yy in a given crystallite with respect 
to the laboratory frame in which the coil axis is the x-axis. Here 
it is assumed that the laboratory frame coincides with the rf 
principal axis system (RFPAS) . We consider only the following 
two terms in the spin Hamiltonian, namely, the quadrupolar Hamil- 
tonian , and the pulsed rf interaction with the spin system 

(K ). Admittedly, this is not a good model to mimic exactly the 
nature of the powder specimen. A similar model has been used by 
Peterson 2 in his calculations of rf pulse responses in powder NQR 
and by Cantor and Waugh in their work on pulsed spin— locking 
phenomenon in pure NQR spectroscopy of spin 1=1 nuclei. 

In what follows, we present the details of the calculation 
of the response of powder samples to rf pulses using the model 
proposed here in order to design appropriate composite n/2 pulses 





for NQR spectroscopy. 
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V ’ A - 2 e mulation of the Resnonsft to P..1 BF 

V.A.2(i) Spin. I a i Cass 

Assuming that the rf pulse is applied along the x-axis of 
the laboratory frame with a frequency and phase <fi , the rf 
Hamiltonian in this frame is given as 

^rf = cos ( w qt +^)[l z sinacos0 + IySinasin/3 + I^cosa] 

.. (v.i) 

Ihe frequency of the rf pulse corresponds to one of the three 
transitioxis , namely, jo)> £ [l +1 > + !-1>] (see Figure III, 1 ) * 

Using the fictitious spin- l/ 2 operator formalism, ^ Eqn. (V. 1 ) can 

be written as 


% r f = cos(W^t +<f>) Jl^sinCtcosfS + Iy 1 sino0siiy3> + I^^costt] 

.. (V.2) 

The time development of the model system will be described 
by following the time evolution of the reduced thermal equili- 
brium density matrix, P(o) (see Eqn. (i . 10) ) , in QIF.^ In QIF 
which is defined by (see Eqn. (II. 12 ) ) , the net Hamiltonian 
defining the spin system is given by 

% * U Q = S rf (t) 

« -^cos^t +fi) sinacos3[l Zl cos(-|-) -I Z2 sin( 

+ sinasin^[ I Yl cos ( — 2 — )~Iy2 sin ( 2 ^ 

+I X2 sin(-^-)] 


^)] 


•sat 1 ) 


+ cosed cos v 2 


.. (v.3) 



1 46 

After truncating high frequency terms in % rf {t ) , we get the time 
independent rf Hamiltonian as 

«V 

^rf = - 2W -,cosa[l Xl cos^ -I X2sin ^] .. (V.4) 

In deriving the above equation, we have made use of the time 
evolution of fictitious spin— 1/2 operators^* under given in 
Table V. 1 . These results have been obtained from the commutation 
relations of fictitious spin- 1/ 2 operators given in Appendix A. 
Though the rf components along all the three axes (x, y and z) of 
the laboratory frame are present in QIF, only the x-component (see 
JEqn. ( V . 4 ) ) is effective in inducing a transition between j 0 ^ and 
T 1+ 1^ + |- 1 ^]energy levels. It is important to note that the 
rf field felt by the spin system (see Eqn. ( V.4 ) ) is independent 
of the polar coordinate 0 whereas it depends on a* Now, at the end 
of an applied rf pulse with time-duration, t^ , the density matrix, 
P( t. ) , is given by 

f( t 1 ) = U rf P( 0 ) U~l = exp(-i^f rf t 1 ) P( 0 ) exp(i^ rf t 1 ) 

= “ (>i [l X 3 cos( 2 a) 1 cosat 1 ) - I X 2 sin( 2 w 1 cosat 1 )] .. (V. 5 ) 

Now f for a particular value of CL we can get the 90 out— of— phase 
single pulse response (refer to Appendix D) as the expectation 
value of at the end of the rf pulse 



* 60^ sin(2W 1 cosat 1 ) 


. . (V. 6 ) 
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Table V. 1 Time Evolution of a few Selected Fictitious Spin-l/2 
Operators under the Influence of the Qu adrup ole 
Hamiltonian, 

u (i U q' - cos( V ) + ! X2 sin(^t) 

U Q I Y1 U Q 1 = X Y1 cos (^f") - Iy 2 siQ ("^”) 

I Zt U Q 1 = I Z1 COs (~2 — ) - I Z2 s in(~2 _ ) 

where * exp(-i 'H^t) 

and 'H q = ^ I X 3 ~ I Y3 I Z3^ * 
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“ e now perfo ™ the -« raelng . fraction of crystallites 
with tdeix x-axe s poiating at an angla batve.n a and a .da with 
respect to the rf coil axis i s 


27T sing da _ -d cos a 

hn ~ 2 . . (V.7) 

The NQR signal from a polycrystalline sample at the end of the rf 
pulse can be obtained by averaging over all possible a values. 

77 ” Wq a =/ r 

s coil / * “ 2 J cos a sin (2a) cosa) d(cosa) 

o 1 

n r sin ^ e 1 ) cosie,) n 

= Q 1 ~2 " 5 •• (V.8) 

6 1 1 

where 9 is the flip angle of the rf pulse defined as 0, = 2o) t . 
The quantity in square brackets in Eqn. (V.8) is J (0 ), the first- 
order Bessel function. The dependence of ^ 1 . on 0^ is given 

in Figure V.2. The value of 0 corresponding to the global maxi- 
mum is 119° (0.66rt). This means that the application of an rf 
pulse with a flip angle of 119° gives a maximum signal intensity 
in the NQR of powder specimens. This is the equivalent of a 
1 7 T / 2 pulse* for a powder specimen. The * TT pulse* for a powder 
specimen is not simply twice the time duration of7T/2 pulse but 

it is 256 ° ( 1.4371 ), since this is the pulse that gives a zero 

2 

signal intensity. This powder averaging procedure is equivalent 


to calculating the rf pulse response in a single crystal first and 
then averaging over all possible values of cc » the orientation 
of this crystal with respect to the rf coil axis. In the next 



Pulse flip angle (0 m Radi axis) 


Figure V.2 Response to a single resonant rF pulse of spin I 
quadrupole nuclei constituting (a) single crystal 
and (B) powder systems. 

a t : (n/ 2 , 1 ) ; a 2 : (*, o); ay (3 n / 2 , -l) 

b, : (0.66rr, 0.436); b„: (l. 43 rt , 0); b„: (1.8871, - 
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subsection, the response of a po ly cry s tall ine sample containing 
spin I = 3/2 nuclei to an rf pulse is discussed. 


V • A » 2 ( li ) Tli© Case of* Spin X as 3/2 

W© now consider an ensemble of non-interacting quadrupole 
nuclei with spin 1 = 3/2, with 7) / 0 and in the absence of a dc 
magnetic field* The diagonal quadrupole Hamiltonian, and the 

rf Harail tonian, r H T ^t QIF, using the tensor operator formalism 
have been derived in Chapter IV (see Eqns. (IV # 9) and (IV. 23 ))* 

Now the density matrix after the application of an rf pulse with 
time duration t 1 is given by 

P(t^) * exp(-i^ rf t 1 ) f(0) exp(i !% rf t 1 ) .. (V.9) 

rv» 

f (t^) may be evaluated using the transformation properties of 
tensor operators given in Table IV. 3 » and is given by 

P ( t ^ ) = -^=r. [ cos 0^ + iT^( s) sin 0^] «. (V.IO) 


Now the signal intensity at the end of a single rf pulse in a 

single crystal case is 


<**> = Tr 



7>> u 


8 


sin B 


.. (v.n) 


ft 


is the flip angle of an rf pulse in the case of a single 


where 



crystal, and is given by 


fi . l^Vi 

( 3 + 71 2 ) ,/2 ' -• ( v - 12 > 

Now assuming that the axis of the rf ooii is oriented with respect 
to the crystallite's efg x-axis by an angled, we get 



g sin 



cosa) . 


.. (v. 13) 


Averaging over all possible orientations of the x-axis of the 
crystallites as shown for spin X — 1 case, we get the NQR signal- 

intensity to be 


^coil) 


Sin 0 1 


8 


cos 0 . 


% 


ft 


i] . 


.. (V.14) 


Generalizing, we see that when we consider the NQR allowed tran- 
sition between a pair of energy levels, the response to a single 
pulse for a powder system containing quadrupole nuclei of spin I 
can be written as 


si gixal intensity 


sin ft cos ft 
CL £ - n 1 

ef 


c . 


.. (V.15) 


The expressions for C and 0 1 are given in Table V.2 for different 
spin systems. This generalization is expected to be of much use 
in the NQR spectroscopy of powder and amorphous samples. 
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Table V, 2 Expressions for C and 0 in Eqn. (V. 15 ) for 
Various Values of the Spin I 



3/2 

ii 

2 A 

Q- 

(l*„ 2 / 3 )' /2 



16 

( 3 * 7 )^* 1 ’ 

5/2 

77 = 0 

U) 

e^qQ 

ko 

for |± I") \± 3 /^y 

ii* w 1 1 1 


(ii) 

— for |+3/2) |+5/2> 

t, 

7/2 

T) ss 0 

(i) 

56 

® 2 qQ for |±l/2) *-+ j+3/2> il5 w 1 t 1 


' i — 

(ii) 7 ~ e 2 qQ for [+.3/2^> | +.5/ 2> iTJ t 1 


(iii) 1 ^e 2 qQ for [ + 5 / 2 ) £-» \+7/z) J7 
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V • B — Design of C omposite Pulses to Enhenf'p the NQR 

Signal Intensi ty of Powder Samses 

As seen in the last section, the flip angle of an rf pulse 
which gives a maximum NQR signal in the case of powders is 0.6S7T , 
unlike in single crystals where it is 0 . 5TT . It should be possi- 
ble to enhance the signal intensity from powder materials by using 
composite 7t/2 pulses that correspond to the optimum single Tt/2 
pulse for powders. Such pulse sequences have to be synthesized so 
that they can uniformly excite different crystallites and thus 
yield good signal intensity. For this purpose, we shall use a 
general numerical approach for designing composite TC/2 pulses. 

Consider a cluster of N number of rf pulses, without any 
time delay between them and with constant rf amplitude (w^) and 
frequency (co^) of the type ( Assuming that 

0 1 = 0 as the reference phase of the sequence, we can define the 

/ complete pulse sequence by ( 2N-1 ) variables* Now, using the 

density matrix approach, we can find out the signal intensity at 
the end of the pulse sequence for NQR powders by applying the 
powder averaging procedure given in the last section. The analy- 
tical expression thus derived for the NQR signal intensity is a 
function of (2N-1 ) variables. Now one has to find out the optimum 
values of the variables corresponding to the global maximum of the 

function, . A 

^ max 

In the present work, we assume a three-pulse sequence, 

l»i>o-lft 8 U -(&,>?> • ^ rea5 ° ns f ° r limlting ° urselves “ 

r 2 3 

three-pulse sequence are: 
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(i) Sequences with a large number of rf pulses, have a larger 

time duration and hence will not be useful in NQR on account 
of the shorter T^. As the pulse time— duration approaches T^, 
the signal intensity decreases. 

(ii) Constructing a pulse sequence with large number of rf pulses 
needs enormous computer time to maximize the function on 
account of the larger number of variables* Moreover, the 
spin dynamics calculation becomes cumbersome as we increase 
the number of rf pulses in the sequence* 

V.J3.1 Spin 1 a 1 Case 


We first consider an ensemble of spin 1=1 nuclei, in a 

single crystal specimen* Using the fictitious spin-l/2 operator 

# 

formalism, we can study the evolution of the spin system under the 
effect of the three-pulse sequence. At the end of the composite 

pulse we have 


signal intensity = <^I X 

= W Q (-sin^ 3 )|^^ 1 sin(2^ 2 )cos^ 3 (l-cos0 2 ) + cos 0, sin& 2 sin( 

+ sin0 1 sin^ 3 (sin 2 ^ 2 + cos^cos^) 
sin( 2 ^ 2 )sin^ 3 cos 0 3 ( 1-cos0 2 ) 

- COS ^sinOgCos 6 3 cos((^ 1 -^ 2 ) 

- sin 6, co s 6 3 cos 0 3 ( sin 2 ^ 2 +cos^cos 2 ^) 
+ sine 3 (sin© lS in0 2 cos^ 2 -cosG^osk,) J 


f sin 

+ cosf^j 2 


.. (V. 16 ) 


where 6 n * 2o)j t n . 
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The signal from the powder sample can be obtained in terms of the 
single crystal results using the powder averaging procedure. 


< X coii> <i x > 


powder 


d(cosC) . . . (V. 17) 


is similar to that for 


The analytical expression for /l \ 

^ x / powder 

^ ^ see ® < l n * (^. 16 )) but with B instead of & , where &' - 


g a cos 


a • V- co ii^ is derived 


n 


to be 


<Xoil> “ V X + *] 


.. (V.18) 


, J - -sin & cos & 


where X = 2 sin^sin^cosy^-r^ 

+ sin2^sin^ 2 [cos^ 2 + sin^ 2 J 

„+ sin^^cos^ 2 ( 1 + cos(^ 2 ) 


^ 0 1 


( f cos(R) sin(R) 1 
' L R “ 0 


R 


and 


Y 3! 


sin^(cos^ 2 ~1 )sin(0«j-0 2 ) [ co ® - SisL§X ] 


1 f cosG sinG . cosH sinH n 
r* t ~ - o + h ~ 9 J 


(V.19) 


2 G 


H* 


/ . \ T cosA sinA ") 

+ ^ e 2 +e 3 “ 2 ®V L A .2 J 


A 


D 


, v r cosB sinB cosD sinD 1 

+ ( e 3 -® Z H ~b“ - b 2 + D ■ 

/ \ f cosC sinC "I 

4* ^ + ©3 + 2©^/ L Q ^ 


(v.20) 
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with 


e , * ® 2 * v 


A « 

B - e, - e 2 + By 

c - 0, + e 2 - e 3 , 

D ' ®, - P 2 - e 3 . 

E - e, * e 2 . 


F = 

6 , - 

S 2’ 


G =r 

e i + 

<V 


H a 

e , * 

6 r 



= cos{ 

f>j sin^ 2 sin( (f>^- (j. 

,)■ 

6 „ 

2 

1 

= ~2 co s 

(0 3 cos(^-^ 2 ) , 


e 3 

= 2 co 

is^cos^cos (^- 

4>> 


and 


2 cos r 3 




The problem of maximizing the function is a non- 

linear one, and it is highly time-consuming. We have developed 
a computer programme (in FORTRAN language) and numerically 
evaluated this function by varying the flip angles and phases 
in steps of 3°. If we use a two-pulse sequence and maximize 


/l . ,\ using Eqns. (V. 18 ) - (V.20) , we get the sequence: 

^ coil / 


(l 26 ) 0 -( 300) g0 with relative intensity 0.44. 

For a three-pulse sequence, we get 

( 90 ) 0 _( 70 ) 45 -( 20 ) 25 vith relative intensity 0.64. 

It is clear from these results that as the number of rf pulses 
in the composite n/2 pulse sequence increases the signal intensity 
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from the powder specimen increases. This could be due to the 
fact that the flexibility available with a two-pulse sequence 
is less than with a three-pulse sequence. However, at the same 
time the duration of the total pulse sequence should be much less 
than T„ . 


V.B.2 General Spin I 

It would be useful if one constructs composite 71/2 pulses 
for spin X * 3/2 and for other spins also. We have followed a 
similar numerical approach for spin I = 3/2 nuclei using the 
irreducible spherical tensor operators to design a composite n/2 
pulse to maximize the NQR signal intensity from powder specimens 
in the absence of a Zeeman field. We have derived the analytical 


expression for for the spin 1 = 3/2 case and found that 
this function is similar to that of spin 1=1 case. The defini- 
tion of flip angle, however, is altered depending on the spin of 
the nucleus (refer Table V.2). For example , in the spin 1 = 3/2 

(3+7])Vn 


case we have 6 = 0 , fZ » in comparison ^ v - . — 1 „ 

n (3+7) 2 ) V 

I = 1 case. Keeping this in mind, we may conclude that the values 
of the variables defining the pulse sequence are the same. Hence, 
a (90) 0 -(70)^ 5 ~(20 ) 25 composite 7t/2 pulse sequence can also be 
used for spin I = 3/2 nuclei in the absence of a magnetic field. 
Arguments similar to those given in Chapter IV can be used to 
generalize the applicability of this composite n/2 pulse sequence 
to an arbitrary spin I quadrupolar nucleus. 

In the next chapter we present some experimental results 
using cm pool to pul... for broadband excitation in Nu* spectroscopy 


to 8 = 2(0 t in the 
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Summ ary 

A numerical approach has been, presented for constructing 
composite It/ 2 pulses for powder specimens to improve the signal— 
to-noise ratio. Using the fictitious spin-l /2 operator formalism 
for spin X = 1 and the tensor operator formalism for spin X = 3/2 » 
the expectation values, namely, <^I \ and thus <^I . have 

been formulated and the duration and phase of rf pulses consti- 
tuting the composite Tt /2 pulses have been optimized. Thus, it 
has been found that the pulse sequence (l26) Q -(300)g 0 and 
( 90 ) o” ( 7° ) 4 ( 20 ) can yield optimum signals and could serve as 

good Tt/2 pulses for powder materials. If the pulse flip angles 
are defined taking into account its variation with respect to the 
spin than these composite Tt / 2 pulses can be employed, in general, 
for an allowed NQR transition. 
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CHAPTER VI 


EXPERIMENTAL DEMONSTRATION OF IHE PERFORMANCE OF 
COMPOSITE PULSES IN PURE NQR SPECTROSCOPY 

In Chapters III and IV, we have derived composite pulses 
for use in pure NQR spectroscopy via the Magnus expansion approach. 
This chapter deals with the experimental verification of the 
functioning of composite pulses in pure NQR spectroscopy. The 
composite pulses we have earlier derived pertain to single crys- 
tal specimens consisting of physically equivalent nuclei. Never- 
theless, as general tests of performance, we have chosen to 
35 

observe the Cl pure nuclear quadrupole resonance of mercuric 
chloride powder and sodium chlorate (both powder and single 
crystal), at room temperature. In Section A of this chapter, 
the experimental set up used for testing composite pulses is 
described. Section B outlines the procedure for evaluating the 
performance of composite pulses, while Section C presents the 
experimental results showing the performance of composite pulses 
in mercuric chloride and sodium chlorate samples. 

VI. A Experime ntal Set-up 

This section briefly highlights the instrumental set up 
used for the experimental verification of NQR composite pulses. 
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For the present experimental work, we need to have an 
instrument with the following facilities. Pulse programmer with 
atleastn /2 phase shifts, a high-power transmitter, a receiver 
with large bandwidth and short dead time, fast digitizer and the 
spectrometer frequency going down to a few MHz. The non-availa- 
bility of a pulsed FT NQR spectrometer with the abovemen tioned 
facilities in our laboratory necessitated the use of a commer- 
cial Bruker CXP 90 FT NMR spectrometer elsewhere.* The following 
features of the spectrometer have been used in our experimental, 
s tudy . 

(1) High power, single solenoid coil f solid state NMR probes. 

( 2 ) High power transmitter including broadband driver ampli- 
fier and rf power amplifier. 

( 3 ) Pulse programmer, which is capable of producing rf pulses 
of shortest duration equal to 0.5 JJs with a minimum dura- 
tion step of 10.0 ns and with phases selectable in Tt /2 

steps . 

(4) Broadband receiver with direct detection (no intermediate 
frequency stages). 

( 5 ) Quadrature phase detection is used which allows one to 
place the transmitter frequency in the middle of the multi- 
ple line NQR frequency spectrum; in this manner, an effi- 
cient use of transmitter power is achieved. 

( 6 ) BC-130 fast digitizer for large NQR spectral widths having 
a minimum dwell time of 0.1 ps in steps of 0.1 p s. 


♦Chemical Physics Laboratory, Central Leather Research 
Institute, Madras, India 
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(7) ASPECT 2000 computer to control the spectrometer. 

VI * B Verifica t ion of Composite Pulses 

In any FT NQR experiment, the very first thing one has to 
do is to establish the relation between the flip angle and the 
time duration of the applied rf pulse. As has been explained in 
the last chapter, this relation is not a general one but varies 
from sample to sample. In the present work, we first performed 
a *7l/2 w pulse width calibration by pulsing close to resonance 
and determining the pulse width resulting in maximum signal 
intensity. For this, the time duration of the rf pulse width 
was varied and the signal intensities compared for different 
pulse width values . Similarly, the "tt" pulse width, which gives 
zero signal on resonance, has also been calibrated. From these 
Tt / 2 and tt pulse widths, we have derived the time durations of rf 
pulses that constitute the composite pulse. 

Now to judge the compensation of efg inhomogeneity or 
resonance offset effects by composite pulses, the response of 
the system to a single pulse was compared with its response to 
composite pulses at different transmitter offsets. The quality 
of a composite Tt /2 pulse may then be assessed by the extent to 
which it overcomes resonance offset and efficiently produces 
transverse magnetization resulting in improved signal intensity, 
as compared to the response to a single pulse excitation. 
Similarly, a better composite n pulse is one which gives minimum 
NQR signal by inverting spin population over a wider offset range 
vhen compared to the single TT poise. Important reeult, proving 
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the compensation of ef g inhomogeneity by composite pulses are 
given in the next section. 

VI * C ]£x P erim enta3 - Results Shoving the Performance of Composite 
Pulses 

VI.C.1 Mercuric Chloride (HgCl,,) 

35 

The two Cl NQR lines at 22.23 MHz and 22.05 MHz were 
first observed in HgCl^ powder sample at room temperature in sx > 

unpublished work by Buyle-Bodin and Monfils, which has been 

1 P 

quoted by Bassompiere in 1953; later Buyle-Bodin and Monfils 

studied the temperature dependence of these two lines. In 1954, 

3 

Dehinel t et al. assigned these two different NQR lines to the 

chemically inequivalent chlorines in the HgClg molecule in the 

4 

solid state. Barnes and Hultsch measured the line widths in a 

powder sample of HgCl^ and reported a value of 1 .2 kHz. The 

5 

asymmetry in efg has been measured in the single crystal and 

6 

polycrystalline material of HgCl 2 through the Zeeman perturbed 
NQR spectroscopy. The asymmetry parameter (rj) values for the 
two different 33 C1 sites and their efg values are summarized in 
Table VI. 1 . The crystal structure of HgCl 2 is given in detail 
in well-known book by Wyckoff, 7 and the relevant data on the 
crystal structure of this compound are presented in Table VI. 2. 
Dehmelt measured the frequency of 37 C1 quadrupole nuclei at 303 & 
to be 17.5197+0.9 MHz. In the present work, we focus attention 
on the 35 C1 NQR spectroscopy. This two-line spectrum (line 
separation 180 kHz) would be an important system for testing the 
compensation of efg inhomogeneity by composite pulses. 
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Table VI. 2 Crystal Structure and Molecular Dimension Data 
on HgCl 2 


Orthorhombic Unit Cell 

% = 4.325 X 

b Q = 12.735 X 
c o = 5.963 X 

Space group is V^ 1 ^ . 

There are four molecules in the unit cell. 

Each unit cell has two planes, which contain two mole- 
cules each, at a distance of (l/4) c q and (3/4) c q . The 
nearest distance between the two chlorine atoms in the 
adjacent layers is 3*35 X. 

HgCl 2 is a linear molecule, and the two molecules in a plane 
are such that their molecular axes make an angle of 77°. The 
molecular axis ( Cl-Hg-Cl ) makes an angle of 38 30’ from b Q 
axis. The Hg-Cl bond distances are 2.23 X and 2.27 X. 

The ionic radii of Hg and Cl are 1.1 X and 1.86 X, 
respectively . 
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A well -powdered specimen of pure HgCl 2 was taken in a lO mm 
flat-bottomed sample tube. First, the CAP 90 Bruker spectrometer 
was tuned (all parts of the spectrometer, namely, probe, preampli- 
fier, and transmitter, were tuned and receiver quadrature balance 
adjusted) to a frequency 22.053 MHz which is very close by one of 
the two NQR lines at room temperature (297 K) . Fixing the trans- 
mitter gain corresponds to an rf pulse power of approximately 
500 Watts; such an rf pulse was applied and the intensity of the 
low frequency line alone measured as a function of pulse width. 

Thus the powder *7t / 2 1 pulse was determined to have a width of 
7.5 microseconds (jUs) and the powder 'fl» pulse a width of 15.0 h s. 

The same measurement was repeated for the transmitter frequency 

35 

re-tuned to 22.233 MHz, which is near the higher frequency Cl 
NQR line, and the duration times of tt /2 andTT pulses have been 
confirmed. For different frequencies of irradiation (from one 
end of the spectrum to the other) by single pulse and by composite 
rr/2 pulses, NQR signals have been recorded after Fourier transform- 
ing the FID and then phase correcting the frequency domain signal. 

A few of them for single n/ 2 pulse and for 4 5 135 composite n/Z 
pulse are shown in Figure VI.1-VI.10. From these figures, it is 
apparent that the 45 135 composite TT/2 pulse excites both the NQR 
lines significantly more uniformly when compared to the single n/Z 
pulse. Therefore, it may be stated that 45 135 pulse is a broad- 
band composite W/2 pulse in NQR spectroscopy. In a similar manner, 
the response to all other composite 7l/2 pulses at different 
resonance offset frequencies have been found to be better than 
to the single */2 pulse. One typical NQR spectrum 


the response 



167 


TltmrBs VI. 1 — VI»10 

35 

Cl pur© quadrupole resonance spectrum of HgCl^ powder specimen 
obtained by the application of single and composite rf pulses at 
various offset ( Aw q) values at room temperature (297 &) . 


Common spectral parameters 

Synthesizer frequency * 

Frequency of irradiation = 
Spectral width 
Number of scans 
Pulse time duration 

Normalization constant 
(NO) 

Zero Hertz in the x-axis 
irradiation frequency. 


22.14 MHz 

22.14 + offset 
455.0 kHz 
20 

7T/2 pulse corresponds to 7 »5 P s 
and 7t pulse corresponds to 15»0 ju s 

-2 unless otherwise mentioned 
in the figure 

scale corresponds to the 





Figure VI. 2 ( a ) 
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Figure VI. 8(a) Single n/2 
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of each composite 7T/2 pulse for the trangffiitter frequency at 22 . l4 

MHz (i.e. at the centre of the spectrum) has been given (Figure 
VI.1t). It is clear that with increase of overall pulse train 
duration its performance deteriorates beyond a point, i.e. the 
NQR signal intensity goes down (refer Fig. Vl.ll(g)). This is 
due to the fact that, when the excitation time period approaches 
Tg , the spin magnetization dephases in the transverse plane and 
hence less signal intensity is observed. Therefore , a composite 
pulse duration time must be much shorter than of the sample but 
it should be loxig enough to nullify the errors of individual rf 
pulses to give a broadband excitation. 

In a separate experiment, the time duration of the 7t/ 2 pulse 
was deliberately misset , that is, it was decreased in value by 
30 percent. With this error, both single 7l/2 and composite n /2 
pulses have been applied. The observed results have shown that 
composite TT/2 pulse gives a uniform broadband excitation even when 
fl ip angle errors of this magnitude are involved. But the overall 
signal intensity in this case, when compared to the response by a 
composite Tt/2 pulse without error, is less. We conclude that -our 
composite pulse is tolerant of mis-settings in the flip angle©". 

The response of single n and composite n pulse 
( 90 ) 0 -( 270 ) 90 -( 90) 0 at different frequencies have been given in 
Figure. VI. 12 end VI. 13. **••• signals have been phase 

corrected vitb respect to a single v/2 pulse NQR spectrum, 
less intense end almost absorption free dispersion signals result- 
ing with (90) 0 -(270) 90 -(9°) 0 composite n pulse, as compsred to 
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single 7T pulse response, confirms th«+ 

» umirms that the composite 7t pulse gives 

“ * ffiCi * nt br0 “ dt * nd spin population inversion. Similarly, 

all other composite * p u l ae s have been tested and one typical NQR 

spectrum of each composite n pulse at 22.14 MHz is given in 

Figure VI. 1 3 . 


VI, C. 2 Sodium Chlorate (NaClOj 

1 :u l 1 1 


This is a standard system and has often been studied by 

jp 

various ^Cl pure quadrupole resonance . 8 The pure NQR spectrum 
of NaClO^ consists of a single NQR line at 29.9125 MHz at room 
temperature. The line width is of the order of a kHz. The unit 
cell of this system is cubic, with length of each side 6.57 & and 
consists of four molecules. The CIO^ ion forms a pyramid with 
the Cl atom at the apex and its three-fold axis of symmetry 
renders the system axially symmetric, withT)= 0 . Each molecular 
axis is parallel to the body diagonal of the unit cell. In the 
present work, we have analyzed the performance of composite pulses 
in both single crystalline and powder specimens of NaClO^. The 
single crystal was grown by the slow evaporation of a saturated 
solution of NaClO^. 


The powder pulse width has been measured to be 10 jU s 

and powder rr pulse to be 20 /js. Similarly, the single crystal 
n/2 pulse has been identified as 9 U* and n pulse as 18^5. Now 
in both the samples, the responses to composite n /2 pulse and 
single Tt / 2 pulse have been compared at different transmitter fre- 
quence s. Surprisingly, responses are identical irrespective of 
the kind of pulse applied i.e., single or composite. This 
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suggests that the ce.posite n/z pul .s 8 are not compensating the 
efg intnniogeneity in NaC10 3 .y.t.., Hov.ver, composite „ pulses 
produced much better population inversion for a wider offset 
range than a single ji pulse in single crystalline and powder 
samples. These results of spin inversion after phase correct- 
ing with respect to a single 7l/2 pulse NQR spectrum are shown in 
Figures VI . 1 4 — VI a 24 for both single crystal and powder samples 
of 1 aClO^. 


VI. C, 3 


The theoretically predicted band width of uniform excita- 
tion by the 45 135 composite ft/2 pulse is approximately 60 kHz 


but the experimental results in HgC^ show that it gives a uniform 
excitation over 180 kHz. This could be due to the excitation by 
the side lobes of the 'Magnus’ response (which is the extent of 
excitation against the offset values, see Chapter III). 


Hxe failure of composite rr/ 2 pulses in the NaClO^ system 
as against its remarkable efficiency in HgCl 2 is a phenomenon 
which we do not understand at present. We surmise that the 
presence of more than one physically inequivalent site and/powder 
specimens leads to behaviour that is not adequately modelled by our 


single site/single crystal theory. It is all the more intriguing 
that, nevertheless, our composite rr/2 pulses work well for HgCl 2 , 
and the composite It pulses work for all the samples we have 
investigated. Finally, because of the practical difficulties, 
including availability of suitable samples, growing single crys- 
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Figure VI. 14- VI. 24 
35 

Cl pure NQR spectrum of NaClC> 3 sample obtained by 
tbe application of single and composite * pulses at 
various offset (Aw q ) values at room temperature (297 *) . 


I ..Spectral Parameters 
Synthesizer frequency 
Frequency of irradiation 
Spectral width 


29.9130 MHz 
29.9130 + offset 
4 00 kHz 


Puls© time durations! 

i) for NaClO^ powder case 

it/ 2 pulse width is 10.0 ps 
and jt pulse width is 20.0 Pa 

ii) for NaClO^ single crystal case 

■n/2 pulse corresponds to 9*0 ps 
and ft pulse corresponds to 18.0 ps 

Normalization constant = -2. 



20 ? 


Figure VI. 1 4 
AO)q * 20 kHz 



Figure VI. 15 
Attq * 15 kHz 



Figure VI. 1 6 
AWq * 10 kHz 



NaCl0 3 


Powder Sample 
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Figure 
ACi)^ = 


Figure 
AC Oq * 


VI. 17 
5 kHz 



(a) (90) 0 -(270) 9 q -(90) 0 \ 

\ 



NaClO^ Powder Sample 
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Figure VI. 19 

A(»)q *■* “20 kHz 




NaClO^ Powder Sample 




Figure VI. 22 
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Figure VI . 23 
* 5 kHz 
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(a) 90 180 270 



(b) 45 270 45 




A (c) 315 225 90 


U , 




f/ 

/I (d) 260 160 80 
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(e) 60 300 60 
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Figure VI. 24 NaC10 3 Single Crystal 


Au^ = 10 kHz 
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a more elaborate studv could not , . 

y couia not be pursued in the present work. 

In this chapter, the experimental demonstration of the 
behaviour of composite pulses has been presented. It is inferred 
that composite ft pulses work well in all the systems we have 
tested, whereas composite rr/2 pulses do not. The excellent 
performance of composite ft/2 pulses in the case of HgCl^ is worth 
pointing out. From our experimental studies, we suggest that 
¥5 135 for composite ft/2 pulse and (90) Q -(270) 90 -(90) 0 for compo- 
site ft pulse ore superior excitation trains in pure NQR spectro- 
scopy. 
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CHAPTER VII 


SUMMARY AND SCOPE FOR FURTHER WORK 

VI 1 . 1 Summary 

A general Magnus expansion approach has been developed and 
e r re < tl v«ly employed in the problem of designing NQR composite 
pulses. Ihia approach is shown to be successful in the context 
of broadband excitation problem for any quadrupole spin system, 
i i respec tive of the spin quantum number. Composite pulses are 
obtained to compensate for efg and rf field inhomogeneities via 
the Magnus expansion approach. These sequences are shown to 
behave as a single rf pulse, namely, (0)0 with a net flip angle 
of 0 and phase 

Our main efforts at constructing NQR composite pulses have 
brought out many interesting results related to the fundamental 
properties of the spin system. Particularly, our design of PACPS 
has led us to infer that the two separate cases, namely, spin 1=1 
pure NqK and spin I * l /2 NMR, are closely related. This has also 
led to the construction of composite TT /2 pulses to compensate for 
quadrupolar and dipole-dipole interactions in the case of spin 
1 « 1 NMR spectroscopy (refer Appendix E). Furthermore, NQR 

composite pulses that have been designed for spin I 1 and 3 / 
cases, are also applicable for any general quadrupole sp' 
nucleus. Those results predict that the study of the system of 

nuclei in the absence of a Zeeman field 


non-interacting quadrupole 
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is basic ill 1 y a two*»l,0v0X problem » 

CoiRjtoetl te pulses designed using the Magnus expansion method 
hav© been examined experimentally by studying the 35 Cl NQR of 
HgCl 2 and NaClO^ samples. Our experimental results support rather 
well the applicability of composite pulses for the uniform spin 
excitation and inversion in pure NQR spectroscopy. 

Apart from this, composite n /2 pulses have also been synthe- 
sized using ** numerical method to enhance the NQR signal intensity 

in powder samples (Chapter V). 

In the next section, we point out the future scope of this 

I*# ii|# 

Will A « 

Vil.2 Jicopa for Further Work 

To improve the performance of composite pulses, one may 
consider high* r-order terms in the Magnus expansion (Chapters II - 
IV) * However, this would entail increasing the number of rf 
pulses constituting a given composite pulse. We should keep in 
mind that the total time duration of a composite pulse should be 

much less than T 0 , Construction of composite pulses which compen- 

* 

sate for two or more parameters (such as efg, rf field inhomogenei 
ties and dipole-dipole interaction) would be welcome from the 
point of view of several important experimental studies in NQR. 
Hies© goals could be attained using computer programs for solvi g 
complicated analytical expressions in the framework of th agn 
expansion technique. It might be possible to eva u g 

order terms in the Magnus expansion following 
approaches. 1 In particular convergence problems with the 
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Mncm.s expansion may also be overcome by Maricq’s approach. 2 Hxe 
evaluation of the rf interaction frame Hamiltonian could be 
simplified by using recursive evaluation procedures. 3 

Tho transformation properties of irreducible spherical 
tensor operators given in Chapter IV would be useful in evalu- 
ating tho performance of any kind of pulse sequences in spin 
1 * i/ *■ NtjH spectroscopy. As an additional exercise to simplify 
the complexity involved in spin dynamics studies in magnetic 
resonance, one can develop a computer program for an easy evalu- 
ation of unitary transformations of base operators constituting 
the spin density matrix. 

In this work, we have designed NQR composite pulses for 
single crystals* 'these composite pulses, when actually produced 
as shown in Chapter VI, perform well in powder HgCl 2 also, even- 
though one cannot assume a priori that pulse sequence designed 
for single crystals will work for powders also. Therefore, it 
may he important, and more worthwhile, to design NQR composite 
pulses via the Magnus expansion approach for powders separately. 
Also, It would be desirable to test the composite pulses designed 
explicitly for powders in the present work (Chapter V). Experi- 
mental work with such composite pulses would point the way 
improve their applicability further. It would be thus necessary 
to employ pulses with phases other than 0 , 90 or gr 

multiples of 90° al®°* 

An exciting area of further research is the use of compo- 
site pulses for investigating NQR of high-T c super conducti ty 
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sample# whose spectra are very broad (~MHz) Composite tt/2 
pulses designed using the numerical approach in Chapter V may 
be employed in order to get better NQR signal intensity in these 
samples. Application of these NQR composite pulses might also 
benefit zero field NQR and NMR spectroscopy."* Relaxation times , 
said Tg f can be measured with better accuracy using composite 
n and n/Z pulses in spin-inversion recovery or saturation 
recovery and spin echo methods. In this context, it would be 
worthwhile to study the problem of designing composite n pulses 
specifically for the purpose of refocussing in NQR spin echo 
experiments. 
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APPENDIX A 


FICTITIOUS SPIN- 1/2 
OPERATOR FORMALISM 

In the case of spin 1^.1, the density matrix of the spin 

system cannot be completely described by only three spin angular 

momentum operators and for this it is necessary to have a set of 

(21+1 ) *”- 1 base operators. For spin 1=1 case, Vega and Pines 

have developed an operator formalism * J which is very compact and 

convenient in studying the time evolution of the density matrix for 

spin I « 1 quadrupolar system. ^ They are traceless, independent, 

Hermitian operators and they form the generators of the group SO ( 3 ) 

These operators in terms of three cartesian linear angular momentum 

2 

operators are given by 


I 


XI 





I 

X2 


_ 1 
" 2 


(I I +1 I ) 
v y z z y' 



(i i +i i ) 

v Z X X z / 



(I I +1 I ) 
x y y x 


9 


Ji / 2 2 

x;3 “ 2 v z y ■ 


Y3 


i (i 2 -i 2 ) 


i 

X Z3 


i (I 


2 2v 

2 '~y x ' 


The two indices p, i in I , denote that for each p, we have a 

Jr **“ 

suhspueo i, = 1 , 2, 3 with spin l/2 transformation properties and 

1 

hence the name fictitious spin- 1 / 2 operators. I , I ? and I 

PH P 3 



A2 

behave 1 ike the cartesian angular momentum operators I , I and 

• x y 

■ l z Tor all three possible p's; p = X, Y, Z. Only 8 of these 
operators are linearly independent and this dependence is express- 
ed by I + I Y3 + I Z3 = 0. 

To simplify the spin dynamics calculation in NQR spectro- 

4 

a copy , some more operators have been defined as given below: 


1 


X4 


1 V „-I 

13 1 } 


» **4 3 I 23“ 1 X3 


J *z4 ■ I x 3 - I r 3 


3 5 

It is necessary to know various commutation relations among 

these operators for evaluating the time evolution of the spin 

systems and these have been summarized in Table A. 1 . We have 

employed this operator formalism for studying the spin dynamics 

6 

of spin 1 NQR spectroscopy in Chapters III and V of this thesis. 
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APPENDIX B 


TENSOR OPERATOR FORMALISM 


For spin 1=3/2, a basis set of 1 5 operators is required to 

define the donsity matrix of the spin system. Here, we have used 

1 — 3 

the irreducible tensor operators for studying the time evolution 

of the donsity matrix in the case of pure NQR of spin I = 3/2 

nuclei (vide Chapters IV and V). Irreducible tensor operators T n 

of rank n and order q have been defined by Buclcmaster et al . ^ via 

4 

the Wlgner-Eckart theorem as 


<Im|T n |lm-> = (-1 * n 1 )<I |1 T n |) I^> 

H \-rri q m'/ 

= (21+1 )" 1 / 2 <^Im‘nq|lm^></l | j T n jj I^> 

) and <^Im ' nq | Irr/)> are the well-documented 3d and 


I II I 


-m q m 


Clebseli-Gordan coefficients f respectively, and the reduced matrix 

4 

element is given by 


<i|i"ii i> 


n 


1 n I ( 21 -Mi 4* 1 ) 


[_ 2 11 ( 2n ) I (2I-n)l 


1/2 


Irreducible tensor operators T n thus defined are given in the 

2 

Table B .1 for n 45 , Later Bowden and Hutchison, introduced the 
aymmfs trie and antisymmetric combinations of T ± “ explicitly as 
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*o 

A «? 


T 1 M 
+ 1 


Tensor Operators for n^5 After Buckmaster et al . 


I - (|) 1/2 [3 X 2 - 1(1 + 1)] i 


? 2 t 1 ,. 1 + 1 - 




T 2 -ll ! 
+ 2 2 + 


T 3 * 


'3 

i ^ a 


(^>' /2 [5 1 1 - {3 1(1 + l)-lj Ij 

* i ^> ,/z it (5^ - ^ - ’> -ij. y. 

^)' /2 2 t 1 *- X l\ 


* i ‘i > ,/2 ^ 


2 ( TO ),/2 [35 -{yo 1(1 + 1) - 25} X 2 + 3 x 2 (l + 1 > 2 - 61(1+1 )] 
* • . -,, f 

+ ("^) ,/2 [ {7 - [3 X (I + 1 ) + 1] I z ], Ij + 

i (7) ,/2 [ {7 I 2 - 1(1 + 1 ) -5}, I 2 ] + 


_4 

T +3 - 


1 <i > ,/2 it V 

4 + 


6 ( 74 ) 1/2 ^ 63 I z + 35 [3 -51(1 +1)} I 3 + {12-50 1 (I + 1) 

+ 15 I 2 (I +l) 2 j I 2 ] 

+ (if) l/2 1 + C 21 I 2 ± ^2 I 3 +7 {9-2 I( I + 1)}I 2 ± 

1 4 {3-1 (I + l)}l z +{12-8 1 ( 1 + 1 ) +I 2 (I + O^f] 

1 (|) l/2 l 2 f3 I 3 ±9 I 2 + {12-XU +l)jl z + {6-1(1 +l)j] 

+ (|) 1/2 X 3 [9 1 2 ±27 I z +{24 -1(1 +1)}] 

i(f ) l/2 x + 4 [x z ±a] 
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T q "(a) - Jyfz (T q n +T_^ ) 

and T q a (a) = JTf 2 (T q n -T_ q ) with q / O. 

These T n are called non-unit irreducible tensors and they 
H 

are orthogonal but not orthonormal . But the unit tensors T n 

q 

2 

are orthonormal and are given by 
T n ^ (an+l) 1 / 2 T n 

q <1 || ^ II I> q 

and Tr f ( T n ) + T “ ' 1 = 6 , 5 

1 v q ’ q J 0 nn ' u qq ' 

where the adjoint is 

(T n )* a (-1 ) q T n . 

v q ' v ' -q 

Since the unit tensors form a complete orthonormal basis set, 
any operator can be expressed in terms of the T^ 1 . For example , 
the density matrix or may be written as 

a = E p n T n 

q q 
11, q ^ ^ 

xi 

where the coefficients P are referred to as Fano statistical 

tensors, given by 


p<r - m 
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Thn non-unit tensors are more convenient to use because they 
hold .for all spin X, We have used the non— unit irreducible tensor 
operators in this thesis. Because T n form a basis set, the 
product of any two tensor operators can be written as follows i 


T n T «= E a (n, n' , N, i) <NQ (n, q ,n \ q '\ T Q N 

q q N q / q 

whore a (n, n' , JM, I) *= (2N+ (_l) N+21 ’ 

fn n* N | (l || T* 1 || i) (i 1 I*' || I> 

\l I IJ <1 || t n || 


{ n n' N 1 

_ _ [ is a Wigner 6 j coefficient. 

1 1 I J 

The commutation relations can be given as 


Tt n , T V'l = T n T - T T n 

'• q ’ q' J q q’ q q 

= 2 a (n, n' f N, i) [l-(-1 ) n+n '- N ] <NQ|n,q,n* , q T N 

NQ x 

3 

In Table B.2, we present the commutation relationships between 
different tensor operators which have been used in Chapters IV 
and V in designing composite pulses for 1=3/2 case. In Table 
li. 3, we present the matrix representation of a few selected 
tensor operators which have been useful in Chapters IV and V for 
representing various spin Hamiltonians in different interaction 


frames « 
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Tab 1 a B.2 Commutation Relationships for I = 3/2, with Symmetric 
and Anti -symmetric Combinations 3 


K 

t 


= Tj(a) 


, t3(.)] . 

=-K 3 - 

(s) 

[T 1 

L O 

f 

Tj(a)l 

= Tj(s) 

fT^s), 

. t 3 ( s )] - 


l>0 

1 

T?(.)] 

» T 2 ( a) 

[tJ(s) ( 

■ T^(a)] = 


p’ 

* 

T?(.)] 

* T 2 (s) 

[ T T (») 1 

1 t 3( s )] = 

-- /^2< a > 



» 

*£(•>] 

* 2T 2 (a) 

!>]<»). 

. T 3(a)] . 



[T 1 

1 O 

» 

TjU)] 

» 2T^ ( s ) 

[Tj(a), 

T 2 ] = 

O J 

■ -^Tf(a) 


[*’ 

f 


* T 3 (a) 

[xj (a) , 

, T 2 ( s ) ] = 

= -^ T o- T 2< 
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Table H . 2 (c on td . ) 
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Tab l o B.3 Matrix Representation of few Selected Tensor 

Operators (s , a) 
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APPENDIX C 


IMPORTANCE OF 'INTERACTION FRAMES' 
IN NQR SPECTROSCOPY 

This Appendix explains how an uninteresting interaction can 

be largely removed by properly introducing the interaction repre- 

1 

sentation. We consider a general case with the complete spin 
Hamil tonian, % , in the laboratory frame as given by 


K 



H 


small 


(C.l) 


where the 'big' and 'small* suffixes denote the strength of the 
respective interactions present in the spin system. 


d P(t) 

dt 


= -i [K, f>(t)] . 


.. (C.2) 


Generally, the time evolution of the spin system due to the action 
of the large interaction ( ^H^g) is well known and may not give 


any 


new information, but investigations on the effects of 


small 


would produce much structural information about the spin system . *- 
So , one is interested in suppressing ^^g* This can be achieved by 
the appropriate incorporation of an interaction representation 
which is defined by if hig * Nov » tiie equation of motion is written 


as 


d_£lii. 

dt 


to (O 


.. (C.3) 



C2 


where tho fcildo (#) represents the new frame , 


nt ) « u big 


.. (C.4) 


'ft * u" 1 *H u, 

big big 


.. (C.5) 


U, , * exp ( -i r H.. t ) 

big * v big ’ 


.. (C.6) 


Eqn. (C,4) can be rearranged as 


f(t) - U blg f(t) U' ig 


.. (C.7) 


This leads to 


df(t) 

dt 


) P <*> * u bi g < 


> U big 


♦ V«*| < 


.. (C. 8 ) 


Using Eqn. (C. 7 ) » we write 


TT 1 = (- itH big> P<*> + u bi g < 


) <■ 


+ Ht) (±*( big ) 


.. (c.9) 


Equating fiqns. (C.2) and (C. 9 ) , we get after cancelling out the 
commutators with 


. u bie ( 


) 


.. (c.10) 



C 3 


Multiplying ifiqn, 
right with 


(C. 10) from the left with U 
gives 


-1 

big 


and from the 


* small * 



rJ 


d P(t) 

dt 


. . (C. 11 ) 


It is apparent from this equation that the total spin populations 

rvJ 

evolve only under the transformed Hamiltonian ( ) in an 

SHI ©JL JL 

interaction ropreaentation. I. effect of * has therefore been 
removed by this procedure. However, it should be remembered that 
the *K is partly transferred to the small Hamiltonian. This 

3 

procedure is widely used in magnetic resonance, for example, QIF^ 

2 2 

in NQR, rotating frame in NMR and toggling frame in AHT. In the 
present thesis, we have used the idea of QIF and toggling frame in 
designing composite pulses. 
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APPENDIX D 


CALCULATION OF THE RESPONSE TO AN ON- 
RESONANCE RF PULSE IN NQR SPECTROSCOPY 

At thermal equilibrium, the spin system is defined by the 
reduced density matrix P(o) (vide Eqn. (i. 10)). Effects of 
applied rf pulses can be analyzed by solving the von Neumann 

i 

equation. In QIF, during 

an rf pulse whose time duration 

is t . the spin system evolves 
tj r-' 

only under c H cf ( *.* Ti rf » ^ ±nt ) ♦ 

In the absence of rf pulse, only 
Internal interactions are opera- 
tive in bringing out changes on the state of the spin system. 
Therefore, at time t, the density matrix v(t) in QIF is given 
by 



IV 

f(t) = exp (-i 7< int t) exp ( 


rf w 


) f(o) exp(i^ rf t w ) exp (; 


+ t) 

int ' 


.. (D.l) 

This can be explicitly evaluated for different spin values using 
the nested commutation relationship with an appropriate basis set 
of traceless independent Herrnitian operators. Now, the phase 
sensitively detected NQR signal (FID) in an experiment can be 
mathematically expressed by the expectation values of the obser- 
vable spin angular momentum operator, i.e.. 



D3 


always out-of-phase NQR signal is detected in an NQR experiment. 

in the present work, while testing the performance of composite 

- 

pulses, we have considered only M^q component. Similarly, 
responses of any kind of rf pulse sequences can be analyzed in 

NQR spectroscopy. 
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APPENDIX E 


PHASE ALTERNATED COMPOSITE n/2 PULSES 
TO COMPENSATE FOR QUADRUPOLE INTERACTION 
IN THE CASE OF SPIN I = 1 NMR SPECTROSCOPY 


In solid state NMR spectroscopy, spectral lines are gener- 
ally broad because of the presence of nuclear dipole-dipole (co^) 

and nuclear quadrupole (to^) interactions in the case of spin 1=1 
1 

nuclei. In addition to this, the influence of the quadrupole 

interaction during rf pulses of finite duration and limited power 

2 

leads to severe distortions in the powder line shape. This 
problem may be overcome by the application of short length rf 
pulses with either high rf power or smaller flip angle. Increas- 
ing the rf power indefinitely is not possible because of tech- 
nical problems, affecting the pulse transmitter and rf probe, 
whereas small flip angles lead to poor signal intensity. There- 
fore, the task of achieving uniform excitation over the range of 
line broadening interactions, quadrupole interaction in particu- 
lar, is difficult using a single rf pulse. To some extent, these 

distortions have been overcome by replacing single rf pulses with 

3 4 5 

composite pulses. First, Tycko and coworkers * J constructed 

' ^ 
composite n pulses using the Magnus expansion approximation for 

broadband population inversion over larger range of g)q and tOj-j. 

It should be mentioned here that composite n/2 pulses are important 

8 9 

in quadrupole echo experiments. * The purpose of using quadrupole 



E2 


©cho technique is to circumvent the dead— time problem asso- 
ciated with single pulse experiments. Furthermore, quadrupole 
echo technique especially in H-NMR, has become a powerful 

means for probing dynamic processes in polymers'* - ^ ^ and biolo- 

i 1 5 

gical materials. It would therefore be very useful to 

construct composite tt/ 2 pulses for achieving quadrupole echo 

1 6 

without phase distortions. For this purpose, Siminovitch et al . 
examined several composite pulses experimentally and they 
concluded that these sequences also introduce additional distor- 
tions with the reduction of finite pulse width effects. Recently, 

1 7 

Dongsheng et al* ' constructed higher-order composite 7l/2 pulses 

for quadrupole echo without phase distortions but only at the 

expense of increasing the excitation period. These sequences 

will not be of much use because effects of chemical exchange 

during composite pulses introduce additional complications in the 

1 6 1 8 

interpretation and analysis of lineshapes in solid state NMR. ’ 
The purpose of this work is to construct short duration composite 

/ 4 

TT/2 pulse sequences by the Magnus expansion approach for quadru- 
pole echo without phase distortions in spin 1=1 NMR case. Here 

1 9 

the irreducible spherical tensor operators are employed to 
simplify the derivation of Magnus expansion terms. In the 
following section, Hamiltonians of different interactions pre- 
sent in the spin system have been given. 

Hamiltonian of the System of Spin 1=1 Nuclei 

Consider an ensemble of spin 1=1 nuclei in the presence 
of magnetic field. The total spin system in the rotating frame 



E3 


of NMR i® given by the net Hamiltonian ft as 

* - *rf + \ * •• <*-l) 

where represents the interaction of rf radiation of 

strength ( 0 ^ (rad/sec) and phase with the spin system. It 

can be written as 

^rf “ U 1 cos ^( t ) + I y ain(^(t)] . .. (E. 2 ) 

1 9 

Using the tensor operator formalism, r fl r f is given by 

*^rf * “ w i ( a ) cos *<*>- IT] (s) sin(j>(t)] .. (E. 3 ) 

The term D denotes the dipole-dipole interaction, which can 
be expressed as 

= d ij ^Zi I Zj " 3 I i ,I j ^ ’• 

where d, . is the dipole-dipole coupling constant between spins 

* J 

i and j . The quadrupole coupling seen by a given spin can in 
first order be replaced by its part that commutes with the 
Zeeman Hamil toni an ^ co q is the L armor frequency. 

\ - w ti U z 2 “ i2 /3) •• ( £ -5) 

where ov. is the quadrupole coupling constant and it is 

expressed as 
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* 41(21-1 ) ^2 ^ cos 2 a- 1 ) + T) COS2J3 sin 2 a] . . (E.6) 

a and $ are the spherical polar angles which define the direct- 
ion of the Zeeman field with respect to the crystal frame. 

5? describes the asymmetry present in the electric field 
gradient surrounding the quadrupole nuclei. The transformation 

properties of ^ and ^ are similar to those of the irreducible 

2 

tensor operator T q with rank 2 and order zero. Therefore, one 
of them here ( i . e . , i® considered, but the result will be 

the same for both the cases. In the next section, a brief 
discussion of the Magnus expansion technique to design composite 

pulses is given. Detailed procedure of designing composite 

4 

pulses have been given elsewhere. 


Derivation of Zeroth-Order Composite Pulses Using the Magnus 
Expansion Approach 

The net evolution operator in the rotating frame can be 
given aa* 


U(t) - u rf (t) o Q (t) .. (e. 7 ) 

U y (t) = T exp[-i / .. (E.8) 


where T is the Dyson time-ordering operator and /^,(t) is the 

20 

total Hamiltonian in the rf interaction frame 
defined by u r f(*)» i.e.. 



which is 



Here , only the phase alternating pulse sequences of the type 
^1 S a **® considered which can be implemented in an experi- 

ment easily without the need for a sophisticated digital phase 
shifter. The overbar in our notation indicates that the corres- 
ponding rf pulse in 180° phase shifted with respect to other 
pulses in the sequence. 0 n is the flip angle of an nth rf pulse 
with time duration t n as given by 0 n * ^ t n . 

The rf propagator u r f ( t) can be written for an n-pulse 

sequence as 


D rfh) - u « rf 

n n- 1 


U rf U rf 
2 1 


U rf " exp [ -ixj (a) cos0 n t n ] 

XX 


.. (e. to) 
.. (e.ii) 


^(t) can be obtained from the transformation properties of 
tensor operators under the effect of U rf (t) which are given in 
the reference 19. The time dependent 7^(t) can be written using 
the Magnus expansion^ as 

# Q (t) = + ^ 2) + ... .. (E.12) 


where is the ith order term in the Magnus expansion. 

Here, only the zeroth-order term in the Magnus expansion is 
considered and it is given as, 

*<°> « \ y df *L(f) .. (E* 13 ) 

^ o 

For a pulse sequence of the type ^ one can 

derive 7^°^ 


to be 



Mb 


*( 0 ) 

<4 


pmrn 

/Jw-jt 



+ ib T*(s) + c T*(s)] 


E6 

(£. 14) 


whe r© 

3r (VW 

a * [ 2- + sin2 ^ + sin( 2 & 2 -2& ) + ^ sin( 2 ^-2 ^+2 ^ )] 


#w f ' 

b * T? [ 1-2 cos2^ + 2003(2 0^2 © 2 ) - cos (2^ + 20 3 - 2© 2 )] 

and 


c 


S8S 


[ -( 0-,+ &2 + ^3) +2sin2 0 1 


+ 231x1(202-20^ 


.. (E.15) 
.. (E.16) 


+ s in ( 2 0 1 + 2 © 3 - 26 2 )] .. (E.17) 


According to the Magnus expansion technique^* ^ a zeroth- 
order composite pulse can be obtained by selecting flip angle 
values in such a way that the sequence should satisfy the condi- 
tion = 0 . For this sequence, = JI ; thus, one can 

achieve uniform excitation over larger range of Qq* From Eqns. 
(E. 15 )-(E. 17 ) one can say that it is not possible to satisfy 
this condition with a three pulse phase alternated sequence. So, 
the condition is relaxed to 

=0 .. (E.18) 


where f(0) = is the thermal equilibrium density matrix. 

This condition demands only’b given by Eqn. (E.l6) to be equal 
to zero because in and Tg( s) commute with f(o) but 

[ ( s ) , P(0)] £ 0. A set of zeroth-order composite H /2 pulses 

derived with b = 0 and with smaller size ' b* have been summ uMK 1 
in Table E.1 along with the corresponding zero 


Table E. 1 Composite tt/ 2 Pulses with Their Zero th-Order Magnus Expansion Terras 
Their Bandwidth of Uniform Excitation Over or 
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terms. Interestingly, when the duration of individual rf pulses 

in the sequences given in Table E.1 , are doubled, one gets a set 

o 1 

of composite 7t pulse for NQR spectroscopy of spin 1=1 nuclei 

which compensate for the efg inhomogeneity . Xt is noteworthy 

that any pulse sequence involving phase shifts only in multiples 

of n radians and which generates a broadband population inversion 

in two-level systems will provide a broadband excitation in the 

7 

enharmonic three-level system if the pulse duration is halved.' 
Since halving the individual pulse durations in an composite ji 
pulse of spin 1 NQR case leads to a composite 71/ 2 pulse in spin 1 
NMR case, the spin 1 NQR system must be a two-level problem. In 
the next section the performance of composite 71 / 2 pulses have 
been illustrated. 

Performance of Composite Pulses 

To test the performance of composite n / 2. pulses, the analy- 
tical expressions for the amplitude and phase of the quadrupole 

22 

echo reported by Barbara have been used for a three pulse 
sequence of the type 9^ ^0^ . They can be given as follows: 

2 2-3/2 

Echo amplitude, A = (A^ + A^) •• (e 19) 

Phase ^1= tan -1 (A^A^) .. (E.20) 

where 

A^ m sin K[cos ftjCos O^sin©^ ~ c°s0^sin( 8^“ )] 

- »ine,»in6 2 sinB 3 si„ 3K 


.. (E . 2 1 ) 
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a 2 * *in 2 K sin & 2 sin ( 6 ^ ^ ) . . (e.22) 

K. « tan* 1 (2^/w^) . .. (E.23) 

When > that is ( for a symmetric phase alternating sequence 

of the form 0 ^ ^ ^ = 0 and hence 0= 0, irrespective of the 

strength of the quadrupole interaction, u>q. Hence, symmetric 
phase alternating composite n/2 pulses can be used for getting 
quadrupole echo without phase distortions. The performance of 
different composite n/2 pulses has been tested by evaluating the 
amplitude of quadrupole echo with respect to m and their band— 

/ Q 

width of uniform excitation is summarized in Table E.i. Admit- 
tedly, the degree of compensation of Qq by these sequences is 
not very good. But one can construct symmetric sequences with a 
larger number of rf pulses so as to enhance the degree of compen- 
sation against 6 ) q as well as to suppress phase distortions. Some 

of the composite 77 / 2 sequences obtained for this purpose from 

21 

two— level NMR broadband inversion sequences are given in 
Table E.2. 

Summ ary 

Construction of phase alternated composite 77 / 2 pulses for 
use in solid state NMR spectroscopy using the Magnus expansion 
procedure is presented. It is important to note that the appli- 
cation of symmetric phase alternated composite n/2 pulses gives 
quadrupole echo with no phase distortions. Furthermore, these 
sequences behave like single rf pulses with a constant net flip 
angle of Jl/2 radians over larger range of toq. The irreducible 



Table JB.2 Symmetric Phase Alternating Composite h/2 Pulses Tor Compensating Quadra pole 
Interactions in the case of Spin 1=1 NMR Spectroscopy 
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spherical tensor operator formalism employed here simplifies the 
derivation of Magnus expansion terms and hence the design of 
composite pulses to a greater extent. It has been shown that 
phase alternating composite n/2 pulses for a three-level case 
can be derived from NQR phase alternating composite inversion 
pulses which combat the effects of efg inhomogeneity in the 
case of spin 1*1 quadrupolar nuclei. Hence, the design of 
composite pulses for spin 1 nuclei in NMR and in pure NQR 
spectroscopies need not to be considered separately. 
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